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ABSTPACT 

In the free wake analysis of a rotor, velocities are computed everywhere 
in the wake and integrated to determine the wake displacement. Newly de- 
veloped techniques of free wake analysis of a hovering rotor are presented. 
These techniques are designed to predict better wake structures and, thus, 
blade circulation distributions. Previous analyses have been 
codified to include the roll up process according to the Betz criteria. 

The near wake is modelled in four different ways, resulting in four different 
vortex sheet models. The results show that the azimuthal location of the 
starting point of the roll up process is important. The computational re- 
sults are compared to the most recent available experimental data. 
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jl. iirrnoDUCTioH 

There exists a need for an efficient computational or analytical tool 
for designing hovering helicopter rotors, aa well as predicting the aero- 
dynamics of these rotors. The need for a free wake model is well recognized 
(references 1, 2, 4 and 5). In reference 1, the wake is modelled by a 
finite number of diacrcte vortex filaments; by .contrast,. the model in 
reference 2 assumes a continuous vortex sheet. 

The model in reference 2 and its associated computer code was originally 
developed to analyze the flow through a wind turbine. More recently, this 
model and the resulting modified computer codes have been applied to the 
study and understanding of rotor hovering aerodynamics (3) . The working air 
is assumed to be inviscid and incompressible. Each rotor blade is represented 
by a liftiug line on which the bound vorticity is assumed to concentrate. 
However, the effect of "spilling" of air around the blade tip (tip loss) has 
been neglected. 

The entire spiral wake system is divided into 3 sections. These are: 

(I) the near wake, (II) the intermediate wake, (III) the far wake, a3 shown 
in figure A-l. All the wakes contain vortex sheets, tip and root vortices. 

The near wake usually extends from Che blade to an azimuth position which is 
a parameter of the model, but i3 usually around 90 degrees. There is a para- 
meter K which determines the distribution of vorticity qt the most outboard 
station between the sheet and tip vortex, and also at the most inboard sheet 
and the root vortex. The mesh size of vortex sheets in the near wake depends 
on the spacing of the blade stations (or nodes). These stations further sub- 
divide the blade into sections. A second set of points on the blade, called 


i'tif t«iVfr 








centers. Is defined at the middle of each section. The intoraadiata wake lias 
a coarser mesh sire and extends to the start of the far vaka. Therefore, cn 
overlapping region called the transition vaka, betvaen tho near and inter- 
mediate vaka , is necessary. The far vaka is noddled by constant radii 
vortex cylinders .in which the radii are determined by the final spen3wioe 
spacing of the intermediate mesh. The vertical position of this vortox 
cylinder is located at one vortex spacing below the last intermediate 
wake vortex. All the induced velocities are computed by integration of tho 
appropriate fora of the Biot- Savert Law. Tho far wake contributions to 
the induced velocities in tho near and intermediate wakes are obtained 
by series solution techniques. Self induced effects are included by the 
method described in Appendix I. 

The induced velocities cannot be computed at the points (nodes) which 
define the elements, such as the tvo end points of a segment element or the 

four corner points of a sheet element. • This is inherent in the usq of 
discrete vortex sheet alcacnts of finite size as described above. Thus, new 
positions on the wake, called centers, need to be defined. At these points 
the induced velocities in the vaka are calculated. Finally, the velocity 
Involved at a node is obtained by linear interpolation of induced velocities 
on ths four centers around the node. 

The computational procedure has to ba an iterative one, aa the deter- 
mination of the induced velocity requires tho knowledge -of the location of 
vorticity, which is not known a priori. Consequently, tho results of the 
semi-rigid wake have been used to initiate the calculation in cha free wake 
analysis. The basic assumption of the 3cmi-rigid wake is that the velocity 
of a vortex element is constant with respect to time and equal to the velocity 
through that point in the rotational plane on which it originated. Another 
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approximatioh .is. naclc by assuming tha angle o£ inclination of tho vortex element 
is constant vith respect to time. 

Some difficulties are encountered in the use of this vortex sheet model; 
for instance, the edge singularities do not cancel between sheets, because 
of their relative displacements as tha wake distorts. In addition, a3 the 
intermediate wake goes bayond 230 degrees, the nodel may not give a con- 
vergent solution, due to the modelling of the roll up process. Finally, there 
is a problem vith the root vortex which will be further elucidated in the 

| 

next section. Thus there is a need for an improvement in the free wake model ! 

| 

to study tha development of powered rotor aerodynamics. j 

The results from these two codels in references 1 and 2 compare favora- 
bly with the available experimental data; however, the computational tech- 
niques used arc time consuming and may, on occasion, not give a convergent 
solution for reasons indicated above. 

♦ j 

However, the simpler model proposed in reference 6 requires only the 
determination of induced velocity at the blades, which are represented by 
lifting lines, and at the infinite vortex lines located below the blade. 

It gives results which compared well with the available test data and also I 

yields useful insight into the aerodynamics of rotating wings. 

In this modified model, it is assumed that a roll up of vortex sheets 
into isolated vortex filaments occurs in the intermediate wake. Because 
these rolled up vortices are discrete vortex filaments, the presence of a 
transition wake is not necessary. A tip loss of .985 as used in reference 8 
has been assumed. A more efficient technique Involving the use of elliptical 
integrals for the induced velocity computations due to the vortex cylinder 
has been adopted. This technique is described in (6). The azimuthal posi- 
tion at which the vortices roll up has been retained a3 a parameter than can 
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II. MODELLING OF THE WAKES AIID DISCUSSION OF RESULTS 


5 


Summary of Frea Wake Models 

MODEL 1 Essentially Che same as ref. 2, in which Che near and 
tnternediata wakes are represenced by a segmented vortex sheet 
connecting a tip and rooC vortex. The far wake is represented by 
vortex cylinders as in ref. 6. 


MODEL 2 The near wake is treated as in Model 1 until on orinuth 
of 4> • 255®, after which it is rolled up into three vortex filaments 
as in rof. 6. The location of these filaments at $ ® 255* 13 determined 
as if the filaments had rolled up immediately at <t> m 0. 


MODEL 3 Same as Model 2, except that the near wake is rolled up 
before first encounter, at <{> ■ 45®, and is assumed to remain in the 
rotor plane until roll up. Tha geometry of the rolled up vortices 
■at 4> ■ AS® is determined as in Model 2. This avoids the lengthy 
computations required to establish the geometry of tha vortex sheets 
and hence is computationally much more efficient than Model 2. 

MODEL 4 Some as Model 3, except that the full vortex sheet model, 
Including tip and root vortices, is used to establish the geometry of 
the vortex filaments at roll up. 

First Vortex Sheot Model 

This model is the original formulation of ref. 2, as explained above and 
in tha Introduction. A complete description is contained in the appendix 
and in ref. 2, The following paragraphs will discuss the exploratory results 
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obtained by varying a vide range of parameters, cueh as viscous core also, 
distribution of vorticiey at tha cost outboard station batveen the vortex 
sheet and ths tip vortex and ths effects of tho root vortex. Ths computed 
vaka geometry and tha rotor thrust coefficients ara compared with ths experi- 
nsntal data of rafarcnca 4. A brief suasary of results can ba explained as 

follows : 


1) Ho important offsets aro caused by increasing the coro depth of 
tho vortex sheets, but tha coefficient of thrust (C^) la increased and tha 
vortical location of tha tip vortex (2^) is decreased by enlarging tha coco 
radius of the vortex filaments. (See Figure 2.) 

2) Suppressing tha strength of the root vortex to sere increases 

tha contraction of the tip vortex. is also increased, but is 

decreased - all slightly. (Sea Figure 3.) In general, tho root vortex has 

negligible affects on the solution, as was also shown in ref. 6. 

3) Figure 3 shows the affect of increasing the aaimuth of tha inter- 

mediate wake from 260 to 360 degree. Koto that tha 2^ location for tha 3S0 
degree case starts to diverge from Che experimental solution. However, at 
260 degrees azimuth, tha 360 degree case shows better agreement with tha 
experimental data. 

Further study has been cade to extend the intermediate wake beyond 
360 degrees, but converged results could not be obtained. One of the 
reasons was tho migration of the rcot vortex up through the rotor. As the 
spacing of tho constant radius vortex cylinder approaches aero, tha induced 
velocity approaches infinity, causing convergence problems even though tha 
vertical velocity and circulation near the root approach aero. As discussed 
above, when the strength of the root vortex is suppressed to zero, no 
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inport one effects either in Z^, C^, or conversed solution are obtained and 

convergence is inproved. 

Another effect causing lack of convergence was probably due to the 
roll up of the tip vortex predicted by the sheet coddling of the intermediate 
vake. Decreasing the mesh size of the intermediate wake could improve the 
coddling of the roll up process. Since a finer cash size of the wake would 
result in excessive computational tine, and since the parallel study of 
reference 6 has indicated satisfactory results with a fast roll up while 
avoiding the convergence problems, the free wake program was codified to 
Include the roll up of the internediate wake, according to the Betz criteria 
explained in detail in reference 6. This led to the second vortex sheet 
codcl discussed next. 

Second Vortex Sheet Model 

In this model, the near and the far wakes arc again represented as 
described in reference 2. The near wake is extended to 255 degree azimuth. 

The roll up of the sheet vortices i3 assuced to take place at the internediate 
wake. The roll up schedule as discussed in reference 5 has been adopted. 
Because these rolled up vortices are isolated vortex filaments, the existence 
°f a transition wake is not necessary. These rolled up vortices adjoin the 
near wake and extend to 720 degrees. 

In this modified model there are three rolled up vortices. They are 
the tip, center and root vortices. It Is logical to roll up the vorticity 
from near the tip (allowing for a tip loss) to the peak of circulation to 
produce the tip vortex. Then, the center vortex is obtained by rolling up the 
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vorticity from the peak circulation to about 80S span where dT/dn is minimum. 
The root vortex contains the vorticity froa about 80Z span to the la3t 
otation of the blado and is kept for this case of three vortex filaments. 

Tha determination of the induced velocity due to the presence of these 
vortices requires the knowledge of their positions, especially the starting 
positions, in the intermediate wake region. However, their starting positions 
are not known a priori; consequently a procedure has to be devf -pad to com- 
pute thea. He have noted in the above that the vortex sheet associated with 
a particular segoent of bound vorticity (there are three such segments) on 
the blade is subdivided into several vortex sheet elements of finite size for 
coaputational purposes. The centroid of the rolled up vortices are computed 
at each <J>, according to the Betz roll up criteria starting froa $ ■ 0, 
and the velocity computed at these points due to the entire wake, consisting 
of the bound circulation, the vortex filaments including self induced effects, 
the rolled up intermediate wake vortices and the vortex cylinders representing 
the far wake. It must be emphasized that in the calculation of the induced 
velocity on the blades, the near wake region is taken to be a continuous 
vortex sheet rather than a few vortex filaments as described above. 

The computational results are shown in figures 4, 5 and 5, compared to 
the test data of reference 8 . Figure 4 shows the effect of the parameter K 
on the bound circulation and location of the tip vortex. * IThcn K equals zero 
the tip and root vortices have zero strength; thus the near wake is represented 
only by vortex sheet elements. However, when K equals 1, the strength of tip 
vortex i3 the same as the amount of circulation at the last center of the 
blade, whereas the strength of the root vortex i3 equal to the value of the 
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bound circulation at the first center of the blada and the cost outboard 
and inboard vortex sheets have zero velocity. K equal to .5 means that tho 
distribution of vorticity in divided equally batuaen the tip vortex and 
the east outboard vortex 3heat, and the root vortex and the moot inboard 
vortex sheet. Koto that the vortex locations shown at tho first encounter 
with the blade ore those of the rolled up vortices although, in t.w=putlng 
blade loads, the full vortex sheet syotea is used since roll up does not 
occur until 255*. 

Figure 5 describes the behavior of tho roll up process. Tho result 
shown is for K equal to 1. The interaction of tha tip vortex with the most 
outboard sheet vortices are studied for the first SO degree azimuth. The 
most outboard sheet vortices and tho intermediate wake tip rolled up vortices 
induce upward and Inward velocity components on the tip vortex, resulting in 
its migration up through the rotor. In turn, the tip vortex affects the 
sheet vortices in the opposite manner. Because of the upward displacement 
of the tip vortex in the starting region of the near wake, its vertical lo- 
cation at the first blade/vortex encounter la slightly above the test value. 
Therefore, the velocity induced at the blade (30Z span), especially tha 
vertical component (down wash), yields lower bound circulation than tho pre- 
dicted value. This roll up behavior causes slow convergence of tha computer 
program. Furthermore, because of tho coarse grid size needed to keep 
computational times manageable, tha roll up procedure is probably only on 
approximate model of the actual roll up. This point is discussed further in 
reference 11. 

Tho results shown are in a reasonable agreement with the test data. 

From tho discussion above, this model suggesto that the roll up of the sheet 
vortices should have been assumed earlier than 255 degree azimuth as indicated 
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by f If, urn 6, whore the near uako la extended only Co 40 degree azimuth, 
resulting in a better agreensnt between the experimental end computed circu- 
lation distribution. This conclusion leads to the third vortex sheet nodal. 


Third Vortex Sheet Kodel 

The near uako representation of the second vortex shoot nodal la 
codified to reduco the computational tlso. Tha near vsko is extended only 
to 40 degrees. Because of its negligible distortion over this distance 
the vertical displacements of the vorticos in this uako are asouaad to 
ba zero. This Implies that tha naar uako is positioned on the rotational 
plane of tha rotor. The results of the seal-rigid uako have boon uccd 
to dotaraino tha spansulse positions of thosa vortices. Thorcfora, the 
daterainaclcn of induced velocities at this near uako is not necessary. 

The starting locations of tha intermediate wake rolled up vortices aro 
treated in tha esaa manner ns in the second vortex oheot nodal. Xt 
oust bo emphasized that tha continuous vortex shoot representation of 
the near wake la now used only to detemaine the induced velocity at tha 
rotor blade from which it springs. • 

The results using this model are shown in figures 7, 8 and 9. The 
number of rolled, up vortices is three and tha K parameter is. equal to 1. 

Figures 7 and 8 compare the results for a two bladed and a four bladed rotor 
with tha experimental data of rafercnca 4. Lack of test data on the blade 
circulation distribution and the uako geometry for a four bladed rotor does 
not enable one to compare results, except for tha values of C^/cr and C^/o. 

As tha number of rotor blades increases, tho affect of first blade/vortcx 
interaction becomes core important. An improvement of tho blade representation 
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Is necessary Co ylald bettor results* Roforenco 7 shows cha results using 
lifting surface theory, end these indeed are better than the lifting line 
solution for a four bladcd rotor. 

.'■lother result for a tvo blcdcd rotor is compared with tho teat data 
of reference C in Figure 9. The predicted value of peak circulation is 

r 

slightly low (•— — ■ .0209). This is caused by tho relative vertical 
fin* 

location predicted for the rolled up tip vortex at ita first encounter with 
tho blade compared to teat data. Thus, an improvement of this vortex sheet 
model is necessary to predict tho tip vortex location core accurately, 
loading to Che fourth nodal. 


Fourth Vortox Sheet Modal 

In this model the displacement- of the rolled up vortices at tine of 
roll up is determined, as in the previous codols, by computing tho velocities 
at the location of a hypothetical vortox filament assumed to rqll up starting 
st ^ * 0 in the near wake. However, in the two previous codols, tho effect 
of tho tip vortex in the near wake on tho tip roll up is essentially neglected 
This is because the induced velocities on cha rolled up vortex linos due to 
the continuous vortox sheet in tho near wake has boon excluded in determining 
their geometry, and tho induced effect of the vortex filaments on each 
other and thcnsalvos used instead. In this codified modol the full vortex 
sheer, including tip and root vortices, is used and not the vortex filament. 

It is also assuacd that the vertical position of the continuous vortex 
sheet representation of the near woke is equal to tho vertical position of 
the tip vortex lino. Consequently, this vortex sheet induces only vortical 
velocity on the vortex lines. Tho spansuiso positions of the near wake 
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vortices are aacuned Co be constant radii vortox filancnts and shoots, of 
uhich ths radii ara the scram lines originating froa tha blade. Therefore, 

Che docoraination of tha induced velocities in tha near wako la not 
necessary, except at tha vortex filesaats. Tha vortox coro also becccas 
critical baesusa of tho close opensuiso positions between tha tip vertex 
end tha tip vortox lino. The ccnpu to clonal results of this codel era 
shown in figures 10 to 16. 

In figure 10, tha nuaber of rolled up vortices ia four. Tho roll up 
schadulo is ss follows: froa the pock of circulation to the near tip region, 

the vorticity is rolled up into tho tip vortex; froa tha peak of circulation 
to 79* span, tha roll up process yields tha first center vortex; then froa 
tho 70S to tho 352 span, the roll up producas tho second center vortex; finally, 
froa 352 span to the last center of tha blade, tha vorticity is rolled up 
into the root vortex. Tho strength of the root vortex has been suppressed, 
for this caso only, to zero in order to avoid tha convergence) problca noted 
above. Fig. 10 shows chat tho effect of varying the asiauthal siza of tho 
near uaka indicates that this is not a highly sensitive paraaetor. 

Figure 11 shows the effect of different roll up schedules for the 
center rolled up vortices only. Tho dashed curve has the new roll up schedule 
as follows: tha first center vortex roll3 up the vorticity froa the peak of 

circulation to the 82.52 span; tho second center vortex contains ths vorticity 

m 

froa 82.52 to 352 span; the root and tip roll up vortices have the same 
roll up schedule as described above. 

Figures 12 to 15 cocpare Che results with test data of reference 4. Tho 
new coll up schadulo (above) is used. Finally, figure 16 shows the effect 
of varying the nuaber of rolled up vortices froa four to the original three 
used in previous codeia and of varying the slope of the lift curve froa 


2tt to 6.15. 
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Summary of Free ,yakc Analysi3 

Tha frea'wako analysis nay ba summarized as follows: 

1) The rotor blade is treated as a lifting lino. 

2) The fluid is assumed to be inviscid and inconpressible. A 
vortex core is used to avoid tho singularities. 

3) The wake is divided into three sections, a near wake attached to 
tho blado, an intermediate wake, end a far wake. 

4) Depending on tho representation of the near wake, different vortex 
sheet nodols can be proposed. 

3) Except in tho first vortex sheet model, the intermediate wake is 
assumed to have rolled up. In the second vortex sheet model, the vortex 
sheets are rolled up after the first blade/vortex interaction. In the third 
and fourth sheet models, the roll up process takes place before the first 
blada encounter. The roll up process is based on Betz theory of conservation of 
momentum. 


6) A peak bound circulation i3 expected near the tip due to the 
close encounter between the blade and the tip vortex froa the proceeding 
blade and this determines the strength of the tip vortex. 

7) The intermediate wake is rolled up into either three or four vortex 

filaments. Tha roll up schedule is explained in section 4. The root vortex 
contribution to the circulation distribution can be. neglected as shown in 
reference 6. . 

8) The far wake i3 represented by semi- inf inite vortex cylinders, 
starting at a vertical distance from the rotor one vortex spacing below the 
last intermediate wake vortex. The induced velocities from these vortex 
cylinders are obtained from the solutions in ref. 6. 
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9) Tho' slops Ilfs curve used ta 2 tt and the convergence criteria ia 
3!! for circulation end 62 for vako geometry. A tip loss of .985 is used. 

10) The intcracdiata wake is extended to 720* for a tuo bladed rotor 
and 360* for a four blade 4 rotor. 
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III. CONCLUSION'S . AIID RECOMMENDATIONS 


This report presents newly developed analytical techniques of free 
wake analysis in a hovering rotor. Results are shown for tho locutions of 
the tip vortex and the bound circulation distribution. 

These results indicate the sensitivity of the blade loading diotri- 
butlon to the location of the tip vortex at the first encounter. Investi- 
gation of different vortex shoot models led to two important conclusions. 
These are: 


(I) the roll up of the vortex sheet has to be included; 

(II) the knowledge of the starting point of this roll up process 
is important. 

These have been studied and yield different wake geometries (especially 
the tip vortex location at first encounter) for various vortex sheet models. 

These methods have demonstrated the ability to predict vortex wake 
geometry and performance of a two and four bladed rotor satisfactorily. 
Recommendations for further investigations include: 

1) Determine the storting point of the roll up process. 

2) Develop more efficient computational techniques. 

3) Use a lifting surface representation of the rotor blade. 

A) Obtain experimental data on wake geometry and on blade circulation 
distribution for four or core bladed rotors. 
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Ft*. 6: Ccsparltcn of exparlaencal and theoretical bound circulation 

distribution dad gecaetry. N«ar wait* ex ter. da Co *0 a arisucti 
and .Integrated. 

Cj * 0.00116, 0 teat " 0.C046, V computed geooatry 
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ri|. 10i Effect of varyini tha aalsuthxl alto of neat v»ko tn 
the fourth vortex eheet tsxlel. 
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Fig. 16: Effect of the slope life curve 

2v Cj - 0.00471 o experinentel Cp “ 0.0046 

6.15 Cp « 0.00465 7 coapuced geoastry (close for both esses) 

C ■ 0.05 , e ■ 0.015 
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Sslf-Elenant Induced Velocity Computation 


Tho Biot-Savart Law shown that a3 the point whore the dctarnination of 


Induced velocity la required approaches ito corresponding line vortex, a 
singularity la expected. Consoquently, a special technique is developed to 
calculate the vertical Induced velocity at Che two end points (or nodes) of 
a curved vortex segment. A basic assumption of this approach ia to consi- 
der the helical vortex filament as a vortex ring. Fron reference 6, the 
vertical component of velocity induced at n due to a ring of radius r Is 


Hal. f iTr rfo— iLA 0 .* 

Jo <n 2 + r 2 - Zrr 


- n cos 3) dd 
■ 2 - 2m cos <p) : 



The limit of integration has been changed to a and 2(ir - a) where the value 

of a is half the size of tho vortex segment. This will avoid the singularity 

that occurs with <? equal to zero and r equal to q. The integral interval 

is the same os the sIzq of a vortex segment. After integrating this 

equation numerically, the result is the sum of vertical induced velocity due 
to that portion of the ring from a to 2(x - a) The self-induced velocitv 

is then obtained by subtracting the above result from the vortex ring 
solution given in reference 8, which derives a self-induced velocity on 
tho vortex ring expressed as 

W a " ( t*(8r/e) - I/ 4 ) 
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AFP EM) IS II 


Tho computer coda i3 written in Fortran language. Tha program takes 
between 10 and IS oinutoa to run on a Honeywell 68/dpo for the third end 
fourth vortex sheet models. Listing of tho subroutines and tho main program 
are found in Table 1. Table 2 defines tha variables in alphabetical order 
and tha output variables are listed in Table 3. A typical solution is 
shown in Table 4 and the overall flow chart i3 given in Table 5. An input 
file is required to run tho program and Tabic 6 gives a typical input data 
file for the case in figure 9. 

Note that in the main program, the number of rolled up vertices and 
the type of vortex sheet model must be specified (tha values of nr!2, nr 11, 
irll, irl2 and imodol arc required). The roll up process is coded in 
the subroutines loop22, loopO and compa2. The convergent test is 3% for 
the bound circulation and 6% for the wake geometry. 
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TABLE 1 


QISXmQL FAGS t5> 
0? PC0.1 QUALIFY 


Main Program (fwmaln2) 

The eain program doesn't carry out any calculation. It controla tha 
execution of tha program by calling tha subroutines. It elco defines the 
parameters such as what type of vortex sheet nodal and how many rolled 
up vortices are in the case. 

List of Subroutines 

input2: This is concerned with the input of the case and verifies 

tha input data. It also computes tha constants used by tha 
other subroutines. 

• 

aenri2: Tha semi-rigid wake of rotor hovering aerodynamics is analysed. 

The results are used to initialise the free wake analysis. 

compa2: The function is to modify the results of the semi-rigid wake 

to include the initial roll up process. 

inivnd2: Computes the coefficients which are required to determine 

the strength of the vortex elements (due to unit circulation 
distribution) 


C0FN2: 


Computes the coordinates of the centers in the near and 
intermediate wakes using the coordinates of the nodes 



loopO: 

i 


loop!2i 


loop22: 


vind2: 
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* 

Dct&raines the induces velocity at the blade centers due to the 
ihtcrred late and far wake vortices 

Cosputes the induced velocity at the blcdo centers iron the 
continuous vortex cheat in tha near wake. It also loops on 
the bound circulation end tests the convergence 

Cccputea the strength of the vortex olesents end ales controls 
the ccaputation of induced velocities everywhere in tho vako 

loop on the nodes where the induced velocities aro computed. 

The entries aro nased as follows: 

Entry vtndb2; ‘ * 

Induced velocities on the blade are detemined. 

Entry vindn2 : 

Induced velocities on tha near wake centers aro detemined 
Entry vindnl2; 

Induced velocities on the near wake rolled up nodes are 
coeputcd 

Entry vind!2; 

Induced velocities on the intenaediata wake rolled up nodes 


are- detemined 



indvel: 


prtv: 


wnofc2: 


intgr2: 


wnofc3: 


intcc3: 


39 

v,. Qsjgsaaa, a> 

' ’ •• ; g? qualstv 

Calculates the induced velocities due to a rectangular segment 
or a semi-infinite cylinder element, depending on the value 
of "it" variables. It has 2 entries: wxy 2 and coord. 

Initialize the induced velocity everywhere in the vak®. The 
entry used is clv2. 

Interpolates the induced velocities at the centers to the nodes 
in the continuous near wake vortices. Computes the total 
induced velocities at the anticipated rolled up near wake 
vortices and at the rolled up Intermediate wake vortices. 

AIgo transforms the x,y,z components o e velocities into the 
radial, tangential end axial components. Tills subroutine i3 
not used for the 2nd vortex sheet model 

Determination of the new vortex element geometry by integrating 
the induced velocities along the stream lines. This subroutine 
is not used for the 2nd vortex sheet model 

Same as la vnofc2, but with the addition of continuous near 
wake vortices. (This is used for 2nd vortex sheet model 
only.) 

Sane as in incgr2, but with the addition of continuous near 
wake vortices. (This is used for 2nd vortex sheet model only.) 
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out2: 


I'H.i 


* J 


Computes Cha performance of cha holicopter rotor and prints 
tha output of cha results 


Input data: Input filo of paraaeters needed to run tha programs 



<; A artAAftftftAftftftftrtftnftft 

C 4 ft • ft 

C ft ft PROGRAM FWCHAIH2 * CBtGSNW* 

c ft ft ft 0? POQ>» 

C ft ftftftftftftftftftftftftftftftftft 

c ft ;j - -- . • 

C ft Ffl'tE-V/AKE ANALYSIS OF ROTOR HOVER I MG AERODYNAMICS 

C ft — 

C ft 

c * 

c ft 12/1/80 

c ft 

c ft SUBROUTINES : 

c * 

C ft- UIPUT2 (PUNCH) 

C * SEHRI2 

C ft C0MPA2 

C * INIVN02 

e » C0FH2 

C * L00P22 

C * VI ND2 (VIHDB2 , VIN0M2,VIN0N12,VINDI 2,EXECN, I H I CM) 
c * I NOVEL (COORO.WXYZ) 

c * L00P12 

c <> PR7V (PRTCEH , PRTNOO , PRTVO * CLV , PRTVN , PRTV I ) 
c ft INTGR2 

C ft WNQFC2 (INIWN2) 

C ft WN0FC3 (IMIWM3) 

C * INTGR3 

c * 0UT2 (0UTINT2J 

c * 

C ft 

C ftiMt***eftitrtftrtftftftftftft**Hrtftft*ftftftftftftft**ft*ftft 
C * 

ccnnton/parm/lwr , I rd, 1 trace, Ipscm.nmes, I trcw, I trcl 1, 1 trc!2, 1 trcg, 

St trcnt, i trcf , I plot, Icgcn, i plotv, i trtg, ivcrgr , i test, issaeb, i samat, 
Sliraa, 1 in l, 1 im2,ni tor, i trons, Ipn.jpunch, iplotw, Isw, Icont, iplotg, 
fiipr.iplott 
ccm.T.on/hparm/1 fmu 

common/gcom/nbldsl.nbids, sigma, fmu.eten (25) ,knnvr,etal (6) .knlvr, 
Sltwi3t,thotad(25) .theta (25) ,thotac(24) , thetaO, thotOd, alphas. 
ScdO,cdk,dpsind,dpsin,dpal id.dpsi I ,coeff .ccoff l,c,s, blades, 
Snnvr,nnva,nnvrl,nnval,nncr,nnca,etanv(25) .ctanc (26) * 
Snivr.nivn.nlvr l , nival, nicr ,nica,etaiv (6) ,etalc(7) , 

Sntva.ntval ,ntcn, f psl, fps2,nrl 1 , nr 12, irl 1, lr 12, i mode I 
common /farc/eta (40) ,wx (22,40) ,wy (22,40) ,w: (22,40) 
equivalence (nnvr.nncrl) , (nlvr.nlcr 1) , (nnva.nncal) 
equivalence (niva.nical) , (nnvrl,nncr2) , (nivr l,nicr2) 
namal iat/nparn/iwr, <rd, itraco, lpse:n,nmes, i trew, I trcl 1, 

SI trcl 2, i trcg, i trcnt, 1 trcf , i plot, icgcn, i plotv, 1 trtg, ivergr. 

Si trons , i pn. j punch, iplotw, i sw, i cont, i plotg, i pr . 1 plott, i fmu, 

S I ims, 1 iml , I !m2, i sem 

c * 

C ft 

c ft DEFAULT PARAMETERS 
c * 

nr I !»25 
nr l 2-27 
i r 1 1-2 
irl2-4 


PASS 13 
QUAUTV 
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iroodel-3 

frdU' t - ln0 

i pr«6 ' v *’ ‘ ‘ ‘ 

Jtraca-0 

Ipscm-l 


itrcll-0 

itrcl2®0 

l trcg«=0 

I trcnt-0 

itrcf**0 

ltrtg-0 

ivorgr-0 

itrens-O 

J punch-0 

lcont**0 

icgcn-0 

rtmas-0 

i fcru D 0 


1 8CJ3®0 


coef f«0. 


Sfpoos piw 


C » 

c * OATA BY DEFAULT 

c * 

nb Ids 1-2 
I twist-0 
alphas-!, 
cceff 1-.5 
1 ims-50 
Hnl-50 
l ira2-15 

c ft 

c ft INPUT OF THE OPTIONS 

road (5tnparss) 

If (iscm.sq. 1) wri to (iwr,55) 

55 format (" SEMI-RIGID HAKE ANALYSIS ONLY FOR ALL CASES",//) 

jcont-icont 
limd-liral 
Itrew-i trcw 
ltrc!2-i trcl2 
laos-nsras 

c ft 
c * 

If ({modal .eq.2) wr I te ( 6 , 1003) 
if ((modal .eq. 3) wr I te (6, IDOL) 
if (Icodol .cq.Jj) write( 6 . 1005) 

c * 

1003 format (///," SECOND VORTEX SHEET MODEL ",///) 

1001, format (///." THIRO VORTEX SHEET MOOEL ",///) 

1005 format (///," FOURTH VORTEX SHEET MODEL ",///) 

C * 

c * LOOP OH THE CASES SUMMITEO TO THE FROGRAtt 

c * — — — — 

C ft 

do 999 ncas»l , 100 

c ft 

c ft INPUT OF THE OATA FOR THE CURRENT CASE. 

C * TERMINATION IF AN END OF FILE IS ENCOUNTERED ON (IRQ) 


c * 


C ft 
C ft 

C ft 
C ft 

C ft 


c * 

c * 
c * 
17 


c ft 
c * 

C ft 
C ft 


102 

c * 
10 


13 

c * 

C ft 
C ft 

c 

c 

C ft 
C ft 

c * 
c * 

C ft 
C ft 

106 

c ft 


call input2(ncas) 

I lml»l imd 
ItrcwHtrcw 
itrcl2“ltrcl2 
nr.es°lmcs 

Jtfir-l 

i: If (jcofit.eq. 1) jter«*nlter+l 

INITIALISATIONS 

If (iscn.eq.l) goto 17 

INITIALIZATION FOR VIND2 

If (isamet.eq.O) call inivnd2 
INITIALISATION FOR WNOFC 

If (Isamet.eq.O. and. imodal .eq.2) call Iniwn3 
If (isamet.eq.O. and. Imcdal .no. 2) call Inlwn2 

If (jcont.ne.O) go to 10 

SEMI-RIGID WAKE ANALYSIS 
contlnua 
call ocrart2 
If (iscm.eq. 1) goto 999 

IF ITEST“1 CONVERGENCE IS HOT REACHED IN SEMRI, 

THE CASE IS TERMINATED 
if (itost.cq. 1) goto 999 

DATA INITIALISATION FROM THE SEMI-RIGID WAKE ANALYSIS 
cal 1 ccmpa2 
call cofn2 

if (I trace.cq. 1) wr I te (iwp, 102) 

format ("*ftft MAIN ftftft NODES GENERATED 3Y CCMPA") 

if (I trace.cq. 1) call prtnod2 

continue 

If (lcont.na.3) goto 13 

call compa2 

call cofn2 

continue 

jcont-0 


OrStGiWAS- PASS 
OF P0C3 i 


MAIN LOOP 

MAP OF THE FAR WAKE 
I test“l 


do 1 ni ter=j ter , 1 im2 

VELOCITIES INDUCED AT THE CENTERS AND NODES 
call !oop22 

VELOCITIES AT THE NODES FROM THE VELOCITIES AT THE CENTERS AND TOTAL VELOCITY 
if (itrace.cq.l.ond.ltrcl2.eq.O) wr i te (iwr , 106) 
format ("l*ftft MAIN ***: VELOCITIES AT THE CENTERS") 
if (itrace. eq.l. and. itrcl2.eq.O) call prtv 

I f (imodel .eq.2) go to 1000 







cull wnofc2 
go to 1001 
call wnofc3 

INTEGRATION (NEW POSITIONS OF THE NODES) 

caH“lntgp5.>C>i>«< 

TTl J LI v> y oc--. “o 

NEW POSITIONS OF THE CENTERS 
call cofn2 
go to 1002 

INTEGRATION (HEW POSITIONS OF THE NODES) 
call lntgr2 


°?poof> V *** 1 


NEW OISTRIDUTION OF CIRCULATION 
call lcopi2 

OUTPUT FOR CONTINUATION 

If (J punch. eq. 1) call punch 

INTERNED I ATE OUTPUT 

If (lecan.ga.l) call cutlnt2 
If (lcgcn.eq.2) wrl to(lwr, 107) 
forest (Ihl) 

If (lcgon.cq.2) call prtnod2 

NO CONVERGENCE IN L00P1, THE CASE IS TERMINATED 
If (ltoat.cq.2) goto 999 

CONVERGENCE 

If (Itcat.eq.O) goto 2 

continue 

NO CONVERGENCE FOR THE MAIN LOOP 
wr I to(iwr, 10i«) 1 I m2 

fornatC'l*** MAIN **<■: NO CONVERGENCE WITHIN",!}," ITERATIONS") 
call out2(l) 
goto 3 

OUTPUT OF THE FINAL RESULTS 


continue 

lout«2 

write (Iwr, 105) nltor 

format ("Iftft* MAUI ***» FINAL RESULTS - CONVERGED AFTER", 1 1», 
4" ITERATIONS") 
do 70 lcd«l,l» 
cdk'sO-S^lcd 
If (led. eq. U) lout*! 
call out2(iout) 
continue 
continue 
continue 


stop 

end 


c * 
c * 
c * 

C ft 
C ft 
C ft 

c * 

C ft 
C ft 
C ft 
C ft 
C ft 

c ft 

C ft 

c * 

C ft 

C ft 

c * 


c * 


10 

5 

11 

9 


subroutine Input2(ncas) 

ftftftftftftftftftftftftftftftftftftftftftft 
ft ft 

ft SUBROUTINE INPUT2 * ’ 
ft ' . > ft 

rf PROGRAM FWC * 

ft v ft 

ftftftftftftftftftftftftftftftftftft'ftftftft 
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12/1/80 

INPUT OF THE OATA FOR THE CURRENT CASE. THEN COMPUTES 
THE CONSTANTS USED BY MORE THAN ONE SUBROUTINE 
EXIT WHEN AM ENO OF FILE IS ENCOUNTERED ON (IRD) 


ENTRY PUNCH: THE CURRENT OATA IS WRITTEN ON (IPN) 

canmon/parm/lwr, 1 rd, l traco, 1 pscm, rvnss , I trcw, 1 trcl I « I trcl 2, . trcg, 

H trcnt, 1 trcf , 1 plot, legsn, I plotv, i trtg, Ivergr , 1 test, I sareeb, i street, 
H\ Jos, 11 ml, l In 2 , niter, I trans, ipn,J punch, iplotw, isw, Icont, Iplotg, 
&ipr , iplott 
corereon/hparn/ 1 fnu 

ccremon/geom/nblds 1 , nblds , s i gma, fmu, a tan ( 25 ) .knnvr ,etai ( 6 ) »knl vr , 

Si twi st, thotad (25) .theta (25) ,thetac(21*) , thetaO, thetCd, alphas, 
6 cdO,cdk,dp 3 lnd,dpsin,dpsi id.dpsi I ,coeff ,coeff l,c,s, blades, 
Snnvr,nnva,nnvrl,nnval,nncr,nnca,etanv(25) ..atane (26) , 
Snlvr.niva.nivr 1, nival, nicr,nica,etaiv( 6 ) ,otaic(7) , 

6ntva,ntval,ntC3,fpsi,fps2,nr 1 l,nr 12, ir 1 1, i r 12, 1 model 
common /resul/gamrec (26) » twx (Zk,2L ) , twy (21*. 2L) . twz (2U,2U ) , 

6wxc(21*) ,wyc(2L) ,wzc(2I») ,wxnc(26,19) ,wync(26,I9) ,wznc(26.1S) » 

Cwxic (9.51) ,wyic(9.51) ,wzic(S,51) ,wrnv(23.18) ,wtnv (28. 13) . 
Swznv(28, 18) .wriv(8,50) ,wtiv(8.50) ,wziv(8.50) . 

Sxnc (26. IS) ,ync (26.19) ,znc (26. 19) ,xic (9.51) ,ylc (9.50 »zic (9,50 * 
Rxnv (28. 18) ,ynv (28. 18) ,znv (28, 18) ,xiv (8.50) ,yiv (8.50) ,ziv (8,50) • 

tia.ib 

equivalence (nnvr.nncri) , (nivr.nicrl) , (nnva.nncal) 
equivalence (niva.nicai) , (nnvr 1 ,nncr2) • (n’vr l,nicr2) 
namol ist/case/nb Ids 1, 1 twist, sigma, knnvr ,nnva, dps ind.etan.knivr, 
Sniva.dpsi id.etai .thetOd.fmu, thetod.alphas.cdO.cdk, 
Sfpsi,fps2,cooff l.ntva.i plot, i plotv, i trans, 

&j punch, iplotw, isw, iplotg, iplott, ifmu 
data conv,twopi/.017L53293, 6.283185308/ 


If (icont. eq. 1) goto 1*0 

read (i rd, case, err *399 »end*888) 

nnva-ntva 

dp3 i n«dps i nd*conv 

dpsi i*dpsi idftconv 

I f (1 twi st . eq . 1) goto 5 

thetaO=»thetOdftconv 

do 10 i“I,knnvr 

theta (i) "thetaO 

continue 

goto 9 

do 11 1*1, knnvr 
theta (i) “thetad (i) *conv 
continue 
cont i nue 


SSSSSfR 


1 


t 








JS£^2b“l 

If (ncaa.na.l.and.nblda.eq.nbliisl) goto 77 

I sesab*0 

nbld3=nbldsl 

bladcs°f loot (nblds) 

c°cos (twopi /blades) 

8“sln(twopi/blc$03) 

cant Cn.ua?, • 


is cast* l 

If (cooff.no. coaffl) goto 61 
do 62 i»l,nnvr 

If (otan(l) .no.otanv(I)) goto 61 

continue 

do 63 i^l.nlvr 

If (otal (l) .na.otalv(l)) goto 61 

contlnua 

goto 61* 

loa.taf'O 

coaf f«eeoff 1 

nnvpofennvr 

nivr»knlvr 

nncr“nnvr+l 

nlca°nlva+l 

nicr**nlvr+l 

nlval»nlva-l 

nnvrl«=nnvr-l 

nlvrl®nlvr-l 

do 65 l“l.nnvr 

etsnv (i) ®otan (i) 

“ contlnua 

do 66 l°l,nivr 
otaiv (i) =>ctai (I) 
contlnua 
do 12 l“2,nnvr 

ettne (l)«(ctanv (l)+otanv(I- 2)) *.5 

contlnua 

do 13 l“2.nivr 

etalc (i)»(etalv(l)+otalv(l-l) ) *.5 
contlnua 

otene(l)«*1.5*etanv(l)-.5*atanv(2) 
otal c{l) *1.5*013 1 v(l)-.5*otaiv (2) 
etanc (nncr)“l.5*otanv (nnvr) -.5 ft otanv(nnvr-l) 
etalc (nlcr)*»l.5 ,; Qtaiv(nivr) -.5*otaiv(nivr-l) 

t 

continue 

if (itrans.cq.O) goto 68 
nnva-nnva+dps 1 1 *2 ./dps I n 

if ( (nnva-ntva) *dpsln. I t.dpsi 1*1.9) nnva»nnva+l 
continue 
nnval“nnva-l 
nnca>*nnva+l 
ntval»ntva-l 
- ntpa**ntva+l 
Vt . . do, 32 i-l ,nncr2 

thetac (i) “theta (l) + (theta (i+l) - theta (I)) * 

6 (etanc (i+1) -etanv (i) ) / (otanv (i+l) -ctanv(i) ) 
continue 


VERIFICATION OF THE INPUT 


- ... - • • ■ - ' ' *>/ 

c * ' • 

C ft NUMBER OF BLADES: LIMIT OF 8 BLADES (SEMRI) 
a- if (nblds.lt. 1. or. nblds.gt. 8) goto 60 
C IHI TAT IONS. BECAUSE OF THE ACTUAL DIMENSIONS OF THE ARRAYS 

Jt ' if (nnvr .gt.25*or. nnvr. It. 2) goto 60 
I f (nivr .gt.6>or .nivr . 1 1.2) goto 60 

if (nnva. gt. LL.or.nnva.lt. 2) goto 60 • 

if (niva.gt.JO.or .niva. It. 2) goto 60 cniG«$^* 

c ft VALUES OF ETA RANGE FROM ETA (ROOT) TO 1. 0 ? pCOH QUALITY 

C ft ANO MUST BE INCREASING 

Ifldbarivd) .It.O.) goto 60 

if (etanv'(l) .nc.etalv(l)) goto 60 

if (etanv(nnvr) ,ne.l..or.etaiv(nivr) .ne.l.) goto 60 

do 80 i**2,nnvi* 

if (otanv(i) .le.otanv(l-l)) goto 60 

80 continue 

do 81 i“2,nivr 

if (ctaiv(i) . le.ctaiv(i-l)) goto 60 

81 continue 

c ft ADVANCE RATIO MUST BE POSITIVE 
if (fmu.lt. 0.) go to 60 
if (icont.no. 3) go to 255 

read(ird,LL3) (wxc(l) ,i®l,nncr2) , (wyc(t) , i*l,nncr2) , 

& (w 2 c (i) , J»l,nncr2) , (gammc (i) , I **2, nnvr) 
gammc (l)-0. 
gammc (nner) *0 . 

LL3 format (8 f 10.6) 

wrl tc(iwr,LL3) (wxc(i) ,i«»l,nncr2) , (wyc(l) ,i-l,nncr2) , 

S (wzc (i) , l B l,nncr2) , (gammc (i) , l»»2,nnvr) 

255 continue 

write (iwr, 150 ) ncas.nbldsl .sigma.fmu 

150 • format ("Iftft* INPUT ftftft CASE 1 } :",I3,/, 

6" NUMBER OF BLADES", i2,/," SIGMA", f 12.5./. 

6" MU". f 12.5) 

if (1 twist. eg. 1) write(iwr, 15») 
if (1 twist.eq.O) wr i to (iwr, 152) 

151 formatC BLAOES ARE TWISTED") 

152 format (" BLADES ARE HOT TWISTED") 

if (I twist.eq.O) wr i tc (iwr , 161) thetOd 

161 formatC PITCH ANGLE*", f 10. 5," DEGREES") 

i f (1 twlst.eq. 1) wr 1 te (iwr , 162) (the tad (5) , i»l .knnvr) 

162 formatC PITCH ANGLE DISTRIBUTION (DG) :",5f 10.5,/, 5«.6f 10.5, 
fi/.5x,Lf 10.5) 

wr I te ( l wr , 153) a I phas , cdO , cdk 

153 formatC STALL ANGLE", f 12.5,/," CD-",f5.L,"+".f5.3,"ftALPHA*ft2") 
wr i te (iwr , 15 L) coef f 1 , I ims, I iml , 1 im2, fpsl , f ps2 

15L formatC ROOT ANO TIP VORTICES FACTOR:", f 10. 5. 

&/ ,"&MAX.SHUMBER60FLITERATI 0NS6F0R6SEMR I ,LOOPl,LOOP2:",3i3,/, 

6" CORE SIZE FOR SURFACE ELEMENTS:", f 10. 5," FOR SEGMENT ELEMENTS:" 

s,f 10 . 5 ) 

wr? te (iwr, 155) knnvr .ntva, (etan(i) ,I-l,knnvr) 

155 formatC NEAR WAKE DEFINITION: (", 12,",", 12,") ,",*f 10. 5, 

L/,5x,6f 10.5,/, Lf 10.5) 

wr i te (i wr , 156) kni vr , ni va, (eta i (i) , i**l .knivr) 

156 formatC INT. WAKE DEFINITION: (", i2,",", i2,") ,",Lf 10. 5, 
6/,5x.6fl0.5./.5x,Lf 10.5) 

wr i te ( i wr , 1 70 ) nnva , ntva 
170 formatC* NNVA". i 3 ," NTVA", i 3) 

wr i te (iwr, 180) isameb, isamet 
180 format (" ISAME3",i3," ISAMLT ",i3) 








asaassBa£&aaaa£s 


If {(Icont.eq.O) .or. (leant. eq. 3)) return 


c ^ p* \3 

c ******* SLirt *^****************^*****^^*** 5 r r ' G ^* ? ^ 

c * * • 

j,0 continue 

rewind ipn ' . . . 

- . lt ,o 4/ .. 1 rjMd(lpn.err-999.end-777)tparB.tgeom i tv..nd.tw<»f 

1 VcsdOpn.e rr®999.o™i-777) trcsul 
n,ji -fscaob-1 

If Cjcont.cq.2) rasd (lrd.ca3e,err«*SS9.end«803) 
return 

igi f ornat (" *** INPUT ***: ERROR. IC0HT**1. EOF OH I PH") 

step 

- ^ 

e **************************************** 

c * 

entry punch 

c * 

c ft HOTeT^IGO UHITATION OF 32K BYTES OF OATA/RECORO 
write (Ipn) tparn,t<jecaj,tvind,twnof 

c * NOTE: , TH| ( S P REWIKO U |S HECESSARY. WITHOUT IT THE LA5T RECORO 
c * IS HOT WRITTEM III CASE OF ABEKO 
rewind Ipn 

If (l trace. go. l) write (Iwr, ISO) 

190 f crcatC *** PUHCH ft**") 

return 

c ft 

c ********* ******************************** 

c * 

8g9 wrlte(iwr,201) 

201 format (lhi) 

stop 

ggg write (Iwr, 200) 

200 format (" *** INPUT *** ERROR OH (IRQ) ) 


999 
200 
c * 
99 
102 
122 


fornat(2l2,al0.5) 
format (8f 10.5) 
format (312) 
stop 

continue 

format (" *** INPUT ft** DATA INVALID OR OUT Or RANGE ) 

wr! te (iwr .case) 

stop 





subroutine semr!2 


e 

c 

e 

c 

c 

c 

c 

e 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 
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A 

A 

A 

A 

A 

A 

ft 

ft 


ft 

ft 

ft 

ft 

ft 

ft 

ft 

A 

ft 

A 


ft 

ft 

ft 


•Sftftftftftftftftftftftftftftftftftftftftftftft 


A ft 

ft . PROGRAM FWC ft 

ft . ft 

ft "SUBROUTI ME SEHR 1 2 * 

ft ft 


ftftftftftrtftftftftftftftftftftftftftftftftftft 


OJwSSiAL PASS Jp 
0 ? PC 03 GUfiLtifr 


12/1/80 

THIS SUBROUTINE USES THE SEMI-RIGIO MODEL 

TO COMPUTE A FIRST APPROX I MAT I OH OF THE GEOMETRY Or THE WAKE 

AND OF THE DISTRIBUTION OF CIRCULATION ALONG THE BLADES 

THE SUBROUTINE RETURNS ITEST1-0 
IF CONVERGENCE IS REACHED IN (LIMS) ITERATIONS 
ITESTl-i IF NO CONVERGENCE 

eomson/parm/iwr , Ird, I trace, lpsem.nmes, i trcw, I trcl l, I tr cl 2, 1 trcg, 

•SI trcnt, I trcf , I plot, legen, iplotv, I trtg, Ivcrgr, I test, I ssraeb, Isamet, 
S liras, I ini, I i m2, ni ter, i trsns, I pn, j punch, iplotw, isw, icont, iplotg, 
Slpr.iplott 
common/hparm/i fmu 

common/ gecm/nb ldsl.nblds.sl graa , f mu , etan (25) . knnvr , e t a I (6) , kn i vr , 

SI twist, thetad (25) , theta (25) , thetac (24) , the tsO, tl.etOd, alphas, 
ScdO,cdk,dpsind,dpsin,dpsi id.dpsi I .coaff ,coeff l,c,s, blades, 

Snnvr , nnva , nnvr 1 , nnva 1 , nncr , nnca , etanv (25) , atanc (26) , 
Snivr,niva,nivrl,nival,nier,nica,etaiv(6) .etaic (7) , 
Sntva,ntval,ntca,fpsl,f‘ps2,nrl l,nrl2,irl l,irl2, i model 
common /resul/gammc (26) ,twx (24.24) ,twy (24.24) ,twz (24.24). 

Swxc (24) ,wyc(24) ,wzc (24) .wxnc (26,29) ,wync(26.2S) ,w 2 nc (26,29) , 

Swxic (9.50 ,wy ic(9,51) »wzic (9,51) ,wrnv (28, 13) .wtnv (28, 18) , 

Swznv (2R.18) ,wriv(8,50) .wtiv(8,50) ,wziv(8,50) , 

Sxnc (26,19) ,ync (26,19) ,znc(26,19) ,xic(9,5U ,yic(9.5U ,zic(9,51) , 
Sxnv (28. 18) ,ynv (28, 18) ,2nv(28,l8) ,xiv(8,50) ,ylv(8,50) , 2 iv( 8 , 50 ) , 
Sia,lb 

dimension alpha (21*) ,wyct(2i*) ,wzct(2I*) ,wyv( 25 ) , 

Swzv (25) . sn (36 . 8) . cn (36 , 8) 
dimension f lamda (21*) ,u (21*) 

rca Mi p (21,) . 1 p, 1 1 ,mu, 1 1 i p (24) .flip (24) .dp (24) ,f dp (24) . tdp (24) , 
6tp(24) ,fp{24) 
dimension bigf (25) ,fmp(25) 
dimension savl (24) ,sav2 (24) ,facv(24) 
equivalence (nnvr.nncrl) , (nivr.nicrl) , (nnva.nncal) 
equivalence (niva.nical) , (nnvr l,nncr2) , (nivr l,nicr2) 
data phi 5. phi 10, pi ,fpi , twopi .pi 2/. 087266463. . 174532929, 

63. 14 1592653.. 07957747151. 6. 283185308,9- 869604404/ 
data Iwx, Iwy, lwz /"WX ","WY ","WZ "/ 

data lota, lu. llan-.da/"ETA ","U ","LAM "/ 
data Ithet, lalpha, Igamma/"THET","AIPH","GAJW7 
data 1 1 ip, ltl ip, If 1 ip/"LIP ","TLIP","FL!P"/ 
data ldp, 1 tdp, 1 fdp/"DP ","TDP ","FDP "/ 
data ltp,!fp/"TP ","FP "/ 

data I bf/"F "/ 
data 1 fac/"FACV"/ 


ARRAY OF SiNES AND COSiNES 


do 10 nb»-l,r.blds 

ps l ° (nb- 1} *twop i /b 1 ads* 

do 10 j«l»36 

phls“phi5"P5i+phi 10* (j-1) 


.. pn.Q a nb)“sin(phis) 
■' eh 1M Vn!>) “cos (ph i *) 
5 r cohtinu8 i t .i , 




</ 

PRAMOTl'S CORRECT! OH* FACTORS FOR THE NUMBER OF BLAOES 


f l d“f rcu-s 1 jpa*p i A . 
do 20 I“l,nncr2 

fp.p (I) “.5* (f Id+sqrt (f ld**2+*!gjna*pi *otanc (i+U *thctac (J) ) ) 
fs°t>ladasft.5*sqrt ( 1 .-t-frap ( i ) **2) /fr.p(i) 
saall f“fs* (l .-otane (i+1) ) 

If (snallf .gt.IO.) go to 1000 
sma-exp (-sisal 1 f) 
go to 1001 
saa“0 . 

blgf (I)«2.*aeos (saa) /pi 
conti nus 

If (lpsca.cq.2) wrl te(iwr, 110) Ibf , (bigf (n) ,n“i,nncr2) 

INITIAL APPROXIMATION 
ltest*0 

•f Lp".5* (f l d+sqrt (f 1 d**2+* i graa*p I *thetaO) ) 
wi“abs (fkp-fou) 
do 31 i"i.nncr2 
wyc(i)“0. 

W 2 C (0 “frcp (I ) -fmu 

facv(t)».8 

continue 

do 32 i't.nnvr 

wyv(i}**0. 

wivOJ^fmplO-fmu 

conti nu • 


SPAHWISE MOMENTUM THEORY 


ca».25 
do 16 k“l , 10 
do 15 l“l#nncr2 
if (k.cq.l) goto 650 

f l«=atan ( (wzc (I ) ^-frau) / (etanc (l+l) +wyc (i) ) ) 

wyc (l) *wyc (I ) *f 1/f 1 srada (i) 

continue 

tg I £3?» (w2C ( i ) +f tr.u) / (etanc ( t+1) +wyc ( I ) ) 
f lssda (i)=atan(tg1ain) 
alpha W larsda (i) -thetac (i) 

If (alphal.gt. alphas) alphal“aiphas 
If (alphal . 1 t. -alphas) alpha l“-alphas 
wzct (i) — (pi 25 ) *s i gma* 

&sqrt ((etanc (i+l)+wyc (i) ) *«2+(frr,u+wzc (i) ) **2) * 
Salphal/(tglam*etanc(i+l)) /bigf (i) 
wzc (i) “v»ct (!) *. 25 +wzc (i) *.75 
wzc (i) “(wzct (i) -ca>«wzc (i) ) / (l.-ca) 
wyc (i) “-wzc (i) *tglam 
continue 
continue 


31 


I f (lpscn.no. 2) goto 13 
wr i to (twr, 110) Iwi, (wjct (n) ,n*»l ,nner2) 
wrltallwr.110) I Ismda, (f Umda(n) ,n«l,nncr2) 
wrlto(iwr, UOl lwy t> (wyc (n) ,n»l,nncr2) . 
wrltodwr.UO) lwz,*(wjc(n) ,n«l,nncr2) 

18 contlnua » * 

c * 

do 17 n*>2,nnvr 1 ORIGINAL PAU& h> 

«aa“(etanv(n)-ecanc(n))/(etanc(rr+l)-atenc(n)tep P003 QUAlitY 
wzv (n) «*wzc (n- 1) + (wsc (n) -wxc (n- 1) ) *ana 
wyv(n)“wyc (n-l) + (wyc (n) -wyc (n-1) ) ftaaa 
17 . continue 

-wyv(l)«>wyc(l) 
wyv (nnvr) »wyc (nncr2) 
wzv (l) -wac (1 ) a. 5 
wzv (r.nvr) “wsc (nncr2) a. 5 

c ft 

c * HAIM LOOP 
c ft 

gammc (1) «0. 
garwjc (nncr) *0. 
do £»00 nl tor«l, 1 Ibis 

c * 

o«. 10*wl 
eps«.0005*wl 
gtnin e O. 
gtnax*0. 

c ft 

c ft CALCULATION OF THE ANGLES AMO CIRCULATIONS 
c * 

• do 1*0 t«!,nncr2 
wxc (I) »0. 
wyct(l)**C. 

W2Ct (!) °0. 

c ft 

C ft CORRECTION TO ACCELERATE THE CONVERGENCE 

f l-atan ( (wzc (i) +fsiu) / (atanc (1+1) +wyc (l) ) ) 
ccr«f l/f lamda(i) 

If(ccr.gt.l.l) ccr«i.l 
1 f (ccr . 1 1. .9) ccr«,9 
wyc (I) -wye (I) Acer 

c ft 

c ft INFLOW ANGLE 

f I smda (I ) -a ton ( (wzc (i)+fmu) / (etanc (1+1) +wyc (I) ) ) 
u(l)“(otanc(i+l)+viyc(i))/co3(f lemda(l)) 

c * 

C * ATTACK ANGLE 

alpha (i) »f lamda (1) -thetac (i) 
c ft STALL EFFECT 

a I pha >a 1 pha (i) 

I f (alpha (i) . la. -alohas) alphal«-alphas 
I f (alpha (i) . gt. alphas) alphal M alphos 

c ft 

C ft CIRCULATION ALONG THE BLADE 

garrwc (i + 1) -pi ftpi ft* i gr.aftu (i ) *a l pha 1/bl ados 

c ft 

C ft STRENGH, FOR THE FAR WAKE 

I f (gtmin.ge.gair,"nc (i M) ) gt. , ain*ganunc (i + 1) 

if (olnaK . It. oam.%* ( i + ll 1 lUntonumi- ( 





*0 


c * 

c ft 

,C 

C ft 


51 


50 

c 

e 

c 




continue 

gtt—gtmex 

If {-gtain.0E.gtf5ax)gtt«»-gtnin 


4>i CAL&>LA$0M OF THE VELOCITIES INDUCED BY THE HEAR WAKE 

A 

do JO cst«l,nnvr 

tg I (wsv (fita) r>fnu) / (e tenv (m) -fvfyv (ra) ) 

> WaEtewsv (ronj+frau 
etav-atanv(frsi) 
da=phl lOftatav 

Ocnsnat D g£mnc (ras+1) -gsswc (en) x t V5 

H-d<j*g8n.«*fpi '’XLtf* 

do 50 nn»l,nner2 
etep-otanc (nrn*i) 
wy2®0. 
wx2-0. 

W32“0. 

do 51 nb-l.nblds 

do 51 JJ-1,9 

sp-sn(jj ,nb) 

cp«cn(jj.nb) 

f-otnp-otavftca 

h»nai2tft(phi I0*jj-phl5) 

rO=»sqr t (f ft t* (etavftap) **2+Hfth) 

thc®daft.5/r0 

os®theftthcft(tgltfEfttglfio+'l.) 
bb®thoft (hfttgl aamstapftjp) /rO 

f j“.5 ft ( (aa+bb) /sqrt (1 .+2.ftbb+ae)+(sa-bb) /sqrt (l.-2ftbb+ao) ) / 
6(aa-fcbftbb) 
fu«ffftfj/(rOftrOftrO) 
wy2»w«2+f uft (tg l toftf -hftjp) 
wx2®wx2+f u* (etavftspfttg I sa*>h*ep) 
w 22 «w 22 +fu«* (otav-atapftcp) 
eontinuo 

wxc (nn) «wxc (nn) +wx2 
wyct (nn) *=wyct (nn)-K*y2 
W 2 tt (nn) “wsct (nn) -H«2 
continue 


52 


* 

* 

* 


CALCULATION OF THE VELOCITIES INDUCED BY THE FAR WAKE 

do 60 ssa-l.nnvr .nnvrl 

tg I as- (wzv (ra) +feu) / (otsnv (sr a) +wyv (era) ) 

wn>2t®wzv (era) +fmu 

etav-atanv (rrsa) 

da-phi lOftotav 

f f»gttftda<»fpi 

if (nsa.cq.l) ff—ff 

do 60 nn-l,nncr2 

otap-etanc (nn*l) 

wx2-0. 

wy2-0. 

wz2«0. 

do 61 I-M5 
atest-0. 

hO»twopi ftf loat (l-I) -ph?5 
nnn-l 

if(l.cq.l) nnn-10 


raraaBaB 



do 62 nbM.nblda 
do 62 jj»nnn,36 
«p°sn(jj,nb) 
cp«*cn(jj ,nb) 
f»etnp-etav*cp 
h«(hO+phi lOAj j) :Wn2t 
rO"aqrt (fftf+(etav*sp) rt*2+h*h) 
fu"f f/ (rO*rOftrO) 
otest-atcst+fu* (otav-etapftep) 
wy2 rc wy2+fuMtglanftf-hftsp) 
wx2“wa2+fu* (etavftspAtglam-hftcp) 
continue 
wz2«wz2+atest 

If (1 .gt. l.and.abs (atost) .It. op?) goto 63 
continue^ 

• ^:dS&inue , 
wxc (nn)*^t«tec (nn)+wx2 
wyct (nn5«»wyct (nn)+wy2 
wzct (nn) *wzct (nn) -fwz2 
continue 

CONVERGENCE TEST 
do 70 n«*l,nncr2 

If (abs (wzct (n) -wzc (n) ) .ge.e) goto 71 
If (abs (wyct (n) -wyc (n)) .go.*) goto 71 
continue 
goto 80 

CALCULATION OF THE NEXT APPROXIMATION 

continue 
wl®0. 

do 55 I°l.nncr2 
wl®wi+wzc(l) 
if (niter. le. 3) goto 5k 
xO" (wzct (I ) *savl (i) -sav2 (!) *wzc (I) ) / 

& (sav! (I) -wzc (i) -sav2 (i) +wzct (I) ) 
f acv (I) «.5* (( (x0-sav2 (i) ) / (savi (i) -sav2 (I) ) )+facv (I) ) 
if (faev (i) . 1 1. . 1) facv(i)*.l 
If (f acv (i) .gt. .9) faev (I) - .9 
continue 
savl (i) “wzc (i) 
sav2 (i)«*wzct(i) 
facn»l.-facv(i) 

wyc(i)«wyc (i) *facv (i)+wyct (I) *facn 
wzc (I) -wzc (i) *f acv (I ) +wzct (I ) *f acn 
continue 

wi«abs (wi) /nncr2 
do 7k n«2.nnvrl 

wyv (n) *=wyc (n- 1) + (wyc (n) -wyc (n- 1) ) » . 5 

wzv (n) -wzc (n- 1 ) + (wzc (n) -wzc (n- 1 ) ) * . 5 

continue 

wyv (1) *wyc (1) 

wyv (nnvr) »wyc (nncr2) 

wzv (1) <=wzc (1) A. 5 

wzv (nnvr)“wzc (nncr2) *.5 
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IS 
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If (lpsem.lt. 2) goto 302 



uo 


302": 
c * 
400 
e * 
c * 
c * 
c * 


102 

80 

c * 

C ft 
C ft 
C ft 


500 


5 

109 

103 


wri tc(iwr, 140) niter 
format (//," ftft* SEWil **ft HITER“",i3) 
wri ta(iwr, 110) twx. (wxc (n) ,n«l,nncr2) 
writa(iwr, 110) Iwy, (wyct(n) ,n*l,nncr2) 
wrlte(Jwr,110) iwz, (wzct(n) ,n«l,nncr2) 
wri to(Iwr. 1 10) Igamma. (gairme(n) ,n«2,nnerl) 
wrlto(iwr^iOO) lfac, (facv(n) ,n«l,nncr2) 
wrlte(iwr, IfiQ) lwy, (wyc (n) ,n«l,nncr2) 

-• '• Write (Iv#r ; 1 10) Iwz, (wac (n) ,n®l,nner2) 

J,V; '.CCOt1<UIO > 

, . ' 

. ecntinOd* • 

it 


34 



CONVERGENCE IS HOT REACHED H1THIH LlttS ITERATIONS 
RETURN WITH ITESTl-l 


lt53t»l 

write (iwr, 102) 

format (" ft** SE/UU ft*ft HO COHVERGENCE ") 
conti nua 

If (Ipsem.eq.O) return 

OUTPUT SECTION 

lt«0. 

lp»0. 

do 500 n®l,nncr2 

cd*cdO+cdkfta 1 pha (n) ft* 2 

l f (aba (a l pha (n) ) . ge . a I phaa) cd»cd-«-edQ 

I ip (n) *u (n) ftgsraic (rt+l) 

tl ip(n)=»l ip(n) ftcos (f l£cda(n)) 

flip (n) *-l ip(n) ftsin (f lsmda (n) ) 

dp (n)»u(n) *ft2ft(piftsigma/blades) ftcdft.5 

f dp (n) »dp (n) ftcoa (f I cmda (n) ) 

tdp (n) *dp (n) fts i n (f 1 zmda (n) ) 

tp(n)"tl ip(n) +tdp(n) 

f p (nj “flip (n) +f dp (n) 

1 t"l t+tp (n) * (ctanv (n+1) -etanv (n) ) 

lp*»lpfrfp(n)ft(otonv(rrt‘l) ft*2-etanv (n) ft*2) 

continue 

lt*ltftbl£dcs/pi 

lp"lpftblcde3/twopi 

if (ifnu.cq.O) go to 5 

cp“lp/fmuft*3 

ct"lt/fmu**2 

ctcp«ct/cp 

continue 

writetiwr, 109 ) 

format (lhl) 

wr 1 te (i wr, 103 ) niter 

format (10a. "SEMI-RIGID YttKE ANALYSIS - CONVERGED AFTER", 14, 
4" ITERATIONS ") 
write (iwr , 120) 

write(iwr, 115) leta, (etanc(n) ,n«2,nncrl) 
writc(iwr, 120) 

wr i te (iwr , 1 10) Iwx, (wxc (n) ,n»l . r>n cr 2) 
wr i te (iwr , 1 10) Iwy, (wyct (n) ,n»l ,nncr2) 
write(Iwr, 110) Iwz, (wzct(n) ,n-l,nncr2) 
wri te( iwr, 110) lu, (u(n) ,n«l,nncr2) 



write (Iwr, 120) 

wr I to (iwr , 1 10) 1 Ismda, (f lamda (n) ,n«l ,nncr2) 
writo(!wr, no) lthet, (thctac(n) ,n“l,nncr2) 
write(lwr, 110) lalpha, (alpha(n) ,n«*l,nncr2) 
write (Iwr, 110) l gemma, (gcmmc(n) ,n**2,nncrl) 
f ' ,;*wr l tc(lwr , 120) 

• v<r l tcjiwr 10) 1 1 ip, (1 ip (n) ,n*»l,nncr2) 
wr fte (iwt\«110) 1 tl ip, (tllp (n) ,n a l,nncr2) 
wr l te ( iwnx. 1 10) If 1 Ip, (PI ip(n) ,rt"l,nncr2) 
wrlte(iwr, 120) 

write (iwr, 110) l dp, (dp (n) ,n«l,nncr2) 
wr l te(iwr , 1 10) I tdp, (tdp (n) ,n*»l,nncr2) 
wrl te (iwr , 1 10) 1 fdp, (fdp (n) ,n*l,nncr2) 
write (Iwr, 120) 

wrl te (iwr , 1 10) 1 tp, (tp (n) ,n«l,nncr2) 
wr i te (iwr , 110) lfp, (f p (n) ,n«l,nncr2) 
if (ifmu.ne.O) write (iwr, 130) it, lp,ct,cp,ctcp 
If (ifmu.eq.O) write (iwr, 130) 1 1, ip 
130 format (//.l Ox, "LT- ",f 10.5,5x."LP- ", f 10.5,/, 10x,"CT- ".flO,5,5x, 

S"CP-S",f 10.5, /. 10x,"CT/CP«t",f 10.5) 
roturn 

c & 

115 format (lx, alt, lx,9fl0.6,/,6x,9f 10.6) 

110 format(lx,al»,lx,9fl0.6./.6x,9f 10.6) 

120 format (lx) 

c * 


55 



0 ? 


end 



subroutine caopa2 

* 

ft ftftftftftft*ftftftftftftftftftftftftftfrft 

ft ft ft 

ft ft PROGRAM FWC ft 

ft ft ft 

ft ft SUBROUTINE C0MPA2 ft 
p . ft ^ ft 

p, ( 'ftftftftftftftftftftftftftftftftftftftftftft 
ft )r , 

* OI»/2l»/78 

ft 

ft INITIALISATION OF THE POSITIONS OF THE NODES 

* ANO OF THE IK'OUCEO VELOCITIES, 
ft FROM THE RESULTS OF SEMAI 

ft 


56 


O? * 


ft IF THE NUMBER OF ROLLED UP VORTICES IS 3 THEN NRL22, 1RL22.CG2.Q2. 
ft ARE NOT USED 

* 


cession/ 1 enb/w I s 
ccssson/rol l6/ql.q2,q3.q4. Iw 

cccr-.on/parn/iwr, ird, i traca, lpaea.mes, I trcw, I trcl 1,1 trc!2,ltrcg, 

&i trcnt, I trcf , I plot, Icgcn, iplotv, I trtg, Ivorgr, I teat, loantsb, istciet 
SI l os, I Ini, 1 la2,ni ter, itrsns, Ipn.jpunch, I plotw, Isw, (cent. iplotg, 
Slpr, I plott 

ccsnon/gccn/nbldsl, nb Ids, signs, frnu.etan ( 25 ) .knnvr.etai ( 6 ) .knlvr, 
Sltwist,thctad(25) .theta (25) ,thatsc{2J*) .thotsO.thatOd, alphas, 
tcdO.cdk.dpsInd.dosin.dpsI Id.dpsi I ,coof f ,ccaff I, c,s, blades, 

Cnnvr, nnva. nnvr l.nnval. nncr, nnca.atanv (25) .etanc (26) . 
Snlvr.niva.nivrl, nival, nl.cr,nica,etalv(6) ,etaic(7) . 
tntva,ntval,ntca,fpsl.fps2,nrl l,nr 12. Irl l,lrl2, i nodal 
cession /reaul/garoc (26) . twx (24.21*) . twy (24.24) . twz (24.24) , 

Cwxc (24) ,wyc (2U) ,wzc (21*) ,wxnc (26, 19) .wync (26. 19) .wznc (26. 19) , 
Swxic (9,51) ,wyle(9,5l) ,w2ic(5.5l) ,wrnv(23.18) ,wtnv (28,18) , 
fiwznv (23, 18) ,wr Iv (3,50) ,wtlv(8,50) ,wziv(3.50) , 
ftxnc (26.19),ync (26,l9).znc (26. 19) ,xic (9.5D ,y ic (9.5U .ztc (9.51) , 
&xnv(28, 13) ,ynv (28, 18) ,znv(28,13) ,xiv(8,50) ,yiv(8,50) ,ziv(8.50) , 
&la, lb 

dinsnsion wzn (28) ,wzi (9) ,dsnl(25) *dgra2 (25) ,dgn3 (25) ,ds*n4 (25) . 
8dgn5(25) .* 1306 ( 25 ) 

equivalence (nnvr.nncrl) , (nlvr.nicrl) , (nr.va,nncai) 
equivalence (niva.nical) , (nnvr l,nncr2) , (nlvrl,nicr2) 


wzc(l)-.8*wzc (1) 

wzc(2)-.8*wze(2) 

nncr3°nncr2-l 

W 2 c (nncr2) (nncr2) 

wzc (nncr 3) ■ . 8*wzc (nncr 3) 

eten!»(otanv(l)+etsnv (nnvr)) ft. 5 

do 1*0 i“2,nncr 

If (etanc(i) .gt.etan) goto 1*1 

ik«i 

if (ctanc (ik+l) -otan. I t.etam-etanc (ik) ) ik**lk+l 

wt«»wyc (ik+l) 

cps 1 "l ,+vrt/etanc (i k) 

ft 

ft MEAN INDUCED VELOCITY 

* 


wzm=0 . 

do 1 i-»l,nncr2 



1 

C ft 

c ft 

C ft 


100 

• c. * 
c ft 
C ft 


20 


101 


102 

103 


21 


500 


91 
90 
505 
510 
c ft 
c * 
C ft 


W2m«w2m+wzc (1) ft (etanv(i+l) ftft2-etanv(i) ftft2) 
continue 

w2n“W2m/(l.-ctanv(l) ftft2) 

APPROXIMATION OF THE EXPANSION OF THE WAKE 

OSttGlHA- 

r3"sqrt ( (fmu+wzm) / (fmu+wzmft2.) ) 0? 

drd 2 “wxc (nncr2) / (fmu+w 2 c (nncr2) ) 

J. ' ^9drdz/(r3-l.) 

wwz^wzm+fmu 

If {I trtfec.ne.O) wrl to(iwr, 100)w2ia,r3 
format (lhl//" ftftft COMPA ftft*: WZM-",f9.&," R3-".f9.6) 

NODES AND VELOCITIES FOR THE NEAR WAKE 

do 20 !*»2,nnvrl 
wzn (i) »fmu+wzc (i -l) + (wzc(i)-w 2 c(i-l)) ft (atanv(i) -etanc (i))/ 
6 (etanc (i+1) -etanc (i) ) 
continue 

wzn (1) »fr,u+w 2 c (1) ft. 5 

wzn (nnvr) »fmu+wzc (nncr2) ft. 5 

do 21 j®l,nnva 

p$i*f loat (j-1) *dpsinftcpsi 

pslc*f loat (j - 1) *dps i n 

z«psi*wwz/cpsi 

vr«drdz*cxp (-amftz) *wwz 

if {imodel .eq.A) go to 101 

aexp^rS- (r3”l •) ftoxp (-amftz) 

go to 102 

aoxpol.O 

ccc°cos (psic) *1. 

sss“sin (psic) *1. 

go to 103 

cce«cos (ps i ) ftaexp 

S3!k«sin(p3i)ftaexp 

do 21 i**l,nnvr 

wrnv(i , j) **etanv (i) ftvr 

wtnv(i ,j)*wt 

wznv (i , j) *wzn (i) -fmu 

xnv (i , j) “cccftetanv (I ) 

ynv(i .j)-sssftetanv(i) 

znv(i , j)»psiftw 2 n (i) 

continue 

If (icont.ne.3) goto 90 
wr i te (iwr, 500 ) 

formatC//," COMPA: WAKE GEOMETRY (Z ONLY)",//) 

do 91 j“l,nnva 
wr I te (iwr.505) j 

wri te (iwr , 510 ) (znv(i,j), i»l,nnvr) 

continue 

continue 

format (/, i A) 

format (lx,9f 10.6,/, 5x,9f 10.6) 

NODES AND VELOCITIES FOR THE INTERMEDIATE WAKE 

do 30 i=l,nivr 
et“etaiv (i) 
do 31 i2°2,nnvr 
i 1” i 2- 1 


mhL>s. 10 
QU^tV 



31 

32 


30 


V * j 


22 


93 

92 


1000 

2000 


c 

c 

c 

c 

c 

c 

c 

c 3000 
c 
c 
c 






0»- o 

o * 9 


k 000 


If (at. lt.etanv(J2)) goto 32 
contl nuo 
continue 

s® (etonv (1 1+1) -at) / (etanv (I 1+1) -etanv (1 1) ) 
jA V (I. l)°a*2nv(l l,ntva) + (i.-a) *2nv(l 1+1 , n tva) 

WzVjl) ciaftwzn (i l) + (t.-a) *wzn(i 1+1) 
contl nije 

P3 H b»dp3 i nft (ntca-2) 
do 22 j®l,nlva 

P9i°(pal ib+f loot (j-l) ftdpsl I) *cpal 
2®P3|£*fcf2/cp3l 
vp«*dPd2*oxp ( - ra^£w2 
sexp®p3-(p3-l.)fto«fxC*CvA2) • . 
ccc»coa (psl) *aar.p . » 

sss®8 1 n (pal) ftoexp 
do 22 t®l,nlvr 
wrtv(l ,j) e ctaiv(l) *vr 
wtlv(i ,j)®wt 
W 2 lv(i ,J)«*w 2 l (I) -feu 
x I v (I ,j)®cec*etalv(l) 
yIv(i,j)®ass*otalv(l) 

2lv(l,j)®2iv(I,l) + (j-l)*dpsii^2l (!) 
contl nuo 

If (Icont.no. 3) goto 92 
do 93 j®l,niva 
wrl to (iwr ,505) j 

wr!to(Iwr,510) (ziv(t,J), i«l,nlvr) 
contl nuo 
conti nua 
nncr5«nncr2+3 

COMPUTE THE STRENGTH OF TIP, ROOT C rtIDDLLE VORTICES 

do 10C0 j“3,nncr2 

k®nncr5”j 

If (aba (garac(k)) . lt.aba (gaK=ac(k-l))) go to 1000 

go to 2COO 

conti nua 

qij— gerac (k) 

q3®gt!nnic (k) -gc-tsc (k— 3) 

ql®gsranc (k-3) -gawnc (2) 

LOCATION OF SECOHO HI DOLE VORTEX 

k3-k-3 

ki»®k-l» 

p®0. 

q®0. 

kl®k-i 

do 3000 l«/,.k4 

dgml (1+i) « (garrac (1+1) -gasac (1) ) ftetanv (1) ft*2 
p®p+dgral (1+1) 

q®q+(gaanc (1+1) -gaansc (1) ) ftetanv (1) 

continue 

cg2-p/q 

LOCATION OF TIP VORTEX 

pl®0. 

p2®0. 

do LOOO l®k,nncrl 

dgm3 (i + 1) - (garonc (i) -gamme (i+1) ) *etanv(i) ft<v 2 
pl“dgm3 (i+I)+pi 

p2®p2+ (gantac (i) -gamac (i+1) ) *etanv(i) 
continue 


S8 




j 




cgi*»pl/p2 

c5»0.0 


c * 

c * LOCATtOH OF FIRST MIDDLE VORTEX 

k3»k-3 

b B 0.0 

do 7000 J-k3.kl 

a- (gairaic (j+1) -g3nm;c (j) ) *etanv(j) **2 ORIGINAL 

b»b+a • ■■ .... G? P033 QSALt'ifv 

c5“c5+(gam:nc (j+l) -garrcnc (j) ) ftetanv (j). 

7000 continue 

cg3-b/c5 

* 

eft LOCATION OF ROOT VORTEX 

c* 

ppl*>0. 

pp2**0. 

k I°k- 1 
ki*-k-l» 

gannic (1) -0.0 
do OOOO i»2,kl» 

ppl»ppl+ (gcrCTc (i+1) -garsne (I) ) ftetanv(I) ft*2 
pp2« , pp2+(gsmrac(i+l) -gatranc (i) ) *etanv (i) 

8000 continue 

cgl-ppl/pp2 
nrl2i-nrl2-l 
c nri22-nr 12-2 

lrl2l“ir 12-1 
c irl22-irl2-2 

do 5002 i»l,nnvr 

If (etanvd) .gt.cgl) go to 5003 

5002 . I c*» i 

5003 tauS^cgi-etanv (ic) 
tau6**e tanv ( i c+ 1) -etanv ( i c) 

W2n(nr 1 l)=wzn(ic) + (w2n(ic+l) -w2n (ic) ) ft (tau5/tsu6) 
c do 500L i*l,nnvr 

c if (etanv(i) .gt.cg2) go to 5006 

c 500J* id-1 

c 5006 tau7-eg2-etanv(id) 

c tau8*etanv (id+1) -etanv (id) 

c W2n (nr 122) **wzn (id) + (wzn (id+l) -W2n (id) ) ft (tau7/tau8) 

do 7001 i«l,nnvr 
if (etanv(i) .gt.cg3) go to 7002 

7001 ie-i 

7002 tau9-cg3*etanv(ie) 
toulO-atanv (ie+1) -etanv (ie) 

W2n (nr 1 2 1) “wzn ( i e) + (w2n (i e+1) -W2n (i e) ) ft (tau9/taul0) 

do 8001 i-l,nnvr 

if (etanv (i) .gt.cgl,) go to 8002 

8001 iq-i 

8002 taull«cgi4-etanv(iq) 
taul2-etanv(iq+l) -etanv (iq) 

wan (nr 12) =wzn (iq) + (wan ( iq+1) -wzn (i q) ) ft (taui l/tau!2) 

do 5007 j-l.nnva 

ps i -float (j — 1) ftdpsin*cpsi 

2-psiftwwz/cpsi 

vr=>drdz>‘‘exp (-arnftz) ftwwz 

acxp-r3" (r 3- l .) ftexp (-amftz) 

ccc-cos (psi) ftaexp 

sss-sin(psi) ftaexp 



i 


| 


5007 


95 

9U 


c 


e 


c 


5008 

10 


97 

96 

C* 

C* 


do 5007 t*nr 1 l,nr 12 
if (i .eq.nr! I) roil«»egl 
if (j .cq.nr!22) rol1»cg2 
i f (i .cq.nri 2 1) roll**cg3 
if (i .cq.nrIZ) roli»cgl* _ * 

wrnv(i , j) <*rol 1*vr 
wtnv(i,j)«wt 
/W 2 ny*(l,,j)T 5 «n{i) -fnu 
• ../,.apv() , j) “GQCftrol 1 
ynv(l ,j)°535ftrol I 
znv(i tj)*psiftw 2 n(i) 
continue 

if (icont.ne.3) goto 9U 
do 55 J”l.nnva 
wri to (iwr,505) J 

write (iwr, 510) (znv(i,j), I— nr 1 1 .nr 12) 

cont i r.uo 

continue 

wzi (irl i)«*wzn(nrl 1} 

wzi (i r 1 2) «wzn (nr 12) 

wzi (Irl2l)««wan(nrl21) 

wzi (irI22)»wzn(nrl22) 

zlv(Irll t l)«=znv(nrl l.ntva) 

ziv(irl2,l) n znv(nrl2,ntva) 

ziv(lrl 2 I,l)° 2 nv(nr! 21 ,ntva) 

Ziv(lrl22, l)°znv(nrl22,ntva) 
do 10 j*»i,niva 

psi"(pst ib+f Icat (j-I) ftdpsi i)*cpsi 

2°P3 ? ftwwz/cps i 

vr«drdzftexp (-araftz) *wwz 

aoxp®r3~ (r3“l .) *exp (-am*z) 

cec^cos (psi) ftaexp 

sss«5in(psi)ftaexp 

do 5008 i“ir 1 1 , i r 1 2 

if (i .eq.irl 1) rolll^cgl 

if (i .cq. 1 r 1 22) rolll a cg2 

if (i .eq. ir 121) roiil“cg3 

if (i .eq. irl2) rol11«cgl* 

wr iv (i ,j)»rol 1 lftvr 

wtiv(i.j) «wt 

wziv(i ,j) =wzi (i)-frau 

xiv(i ,j)*>cccftrol i 1 

yiv(i ,j) “assftrol 1 1 

ziv(? ,j)nziv(i , 1) + (J-1) *dpsi iftv/zl (i) 

continue 

continue 

If (icont.no. 3) goto 26 
do 97 j*l.niva 
write (iwr, 505) j 

write (iwr,510) (ziv(i.j), l»trl t,ir!2) 

continue 

continue 


o«"£ 0 o» ^ 

0 Tt ' 


*elf induced velocity 

cl=0. 

c2-0. 

r 1-.95 

r2»l .05 

do 1009 kw« 1 , 2 
do 1009 i-1, 180 



61 


1010 

too? 

c* 


p S rnjft0.03U2+0.3U2 

lf(kw.eq.2) go to 1010 
b«i-r l*cos (psr) 

d« (1+r lftrt2-2*r l^cos (p3r) ) .5 

cl-cl+(b/d) *0.0311*2 

go to 1009 

b«l-r2*cos (psr) • 

d- (l+r2**2,-2*r2*C03 (psr) ) .5 

c2-c2+(b/d)**0.031/*2 

continue 

wl»(cl+c2)/2. 

Lamb induced velocity (vortex ring) with d».0l,D=»2 

wi 1-5.52 

wlc-wi 1-wi 

return 

end 



subroutine lnivnd2 


ftftftftftftftftftftftftftftrtftftftftftftftft 
* ft 

ft PROGRAM FWC * 

t{ A ^i. r 't^’j ft 

if -'SUBROUTINE INIVHD * 
ft ft 

ftftft**ftft**ftftft*ftftftftftft*ftftft 


«&$* 


of^' 

O* v • 


02/20/78 

PARAMETERS FOR THE EVALUATION OF THE STRENGH OF THE ELEMENTS 
FROM THE DISTRIBUTION OF C i RCUIAT I ON ALONG THE BLADES 

FOR EACH ELEMENT (I) 

HMl(l.J) IS THE STRENGH OF THE ELEMENT ON THE INNER SIDE 
MULTIPLIED 8Y THE SIDE LENGTH OF THE ELEMENT, FOR A UNIT 
CIRCULATION AT THE INNER CENTER (J) OF THE BLADE 
HH2 (I , J) IS THE STRENGH OF THE ELEMENT ON THE OUTER SIDE 
MULTIPLIED BY THE SIDE LENGTH OF THE ELEMENT, FOR A UNIT 
CIRCULATION AT THE INNER CENTER (J) OF THE BLADE 

C(KEson/cvl nd/ha In {2L.2L) ,hra2n (21*. 2L) .hcrnnjtl (21*) ,hscat2 (21*) ,her--nt3 (21*) 
cemaon/parni/lwr, i rd, i trace, Ipsca, antes, i trcw, i trcl 1, 1 trc!2, 1 trcg, 

SI trcnt, itrcf , I plot, icgen, iplotv, itrtg, ivergr, i test.isamab, isartet. 
Sties, I leal, I ica2,ni ter, i trans, ipn.j punch, iplotw, isw, iccr.t, ipiotg, 
Sipr.Iplott 

eocoon/gecra/nbldsi, nblds, sigma, fnu.ctan (25) .knnvr.etai (6) ,Univr, 
Sltwist,thetad(25) , theta (25) ,thetac(21*) .thetaO.thetOd, alphas, 
ScdO,cdk,dpsind,dpsin,dpsi id.dpsi i ,coeff,coaffi,c,s, blades, 
Snnvr.nnva.nnvr l,nnval,nner,nnca,Qtanv(25) ,ctanc(26) , 
Snivr,niva,nivrl,nival,nicr,nica,etaiv (6) ,etaic (7) » 
Sntvn,ntval,ntca,fpsl, fps2,nr l l,nr 12, irl 1, ir 12, i model . 
d Intension gne (30) ,gic(7) ,hntv(2l») ,hitv(6) ,gntv(25) ,gitv(2i*) 
equivalence (nnvr.nncrl) , (nivr.nicrl) , (nnva.nncal) 
equivalence (niva.nical) , (nnvrl,nncr2) , (nlvr l,nicr2) 
data Ihli , Ih2t/"HHII","HM2I"/ 
data Ihln, ih2n/ ,, HHlN ,, ,"HM2N ,t / 

data 1 ht 1 , « ht2, 1 ht3. i hti*/"HMTl", "HMT2 ,, ,"HMT3" ."HHTi,"/ 

CLEAR ARRAY OF UNIT CIRCULATION ALONG THE BLADES 

do 10 i*»i,nncr 
gne (i) “0. 
continue 

LOOP ON THE UNIT CIRCULATIONS ALONG THE BLADES 

nncr3*>nncr2-l 
do 12 j«l,nncr3 
gnc(j+l)«l. 

LUMPED STRENGH OF THE TRAILING VORTICES 

do 11* i«l,nnvrl 
gntv(t) «gnc (i+1) -gne (i) 
continue 


DISTRIBUTED STRENG OF THE TRAILING VORTICES 


c * 


nnvr2»nnvr 1-1 

do 15 i»2,nnvr2 QSKQiftAl. PAGE tS 

hntv (1) "gntv (i) / (etanc (i + 1) -etanc (i) ) q~ PCOH QUALITY 

.15 continue 

i . tttl-gntv (1) *coeff 
w ''»^hC2"gntv (nnvr 1) >'cocf f 

’hnt\>(l)’ a 2;fc (gntv(l) -htl) / (etanc (2) -etanc (1) ) 

hntv (n'nv?l0'“2.* (gntv(nnvr 1) -ht2) / (etanc (nncr 1) -etanc (nncr2)) 

c ft 

C * COEFFICIENTS FOR THE NEAR WAKE 
c * 

hsnmtl (j)“htl 
hamnt2 (j)«ht2 
do 18 i"l,nnvr2 

hmln(i ,j) n hntv(i)ft(ctanv(i+l)-etanv(i)) 
hn2n(i , j) H hntv (i+1) ft (otanv (i+1) -otanv (i) ) 

18 continue . 

c* 

c * 

gnc(j+l)“0. 

12 continue 

c * 
c* 
c * 

c ******************** 
c ft 

if (i trace. eq.O) return 
write (iwr, 150) 

150 formatC *** I Ml VI HD ft**") 
wr i te (iwr , 151) 

151 • format (/," I NT.COEFF . (NEAR WAKE)") 

do 68 ki“l,nnvrl 

wr i te (iwr , 1 10) Ihln, (hmln (k i , k j) ,kj**l,nncr2) 
write(iwr,110) lh2n, (hm2n(ki ,kj) ,kj«l,nncr2) 

68 continue 

write (iwr, 110) 1 htl, (hammtl (kj) ,kj»l,nncr2) 
wri te (iwr , 1 10) lht2, (hammt2 (kj) ,kj«*l,nncr2) 
write (iwr, 152) 

152 format (/," I NT COEFF. (INTERMEDIATE WAKE)") 
return 

c * 

110 format (lx, al,, lx, 9f 13.8../, 6x,9f 13-8) 

120 format (lx,a4, lx,9i 13./.6x,9i 13) 

a * 



end 


subroutine cofn2 

c a 

-C ftftftftftftftftftftftftftftftftftftftftft 

~ r ,#wn 1 -.»•»*#., i; \ 

C3 “-•* ■* 

c ft * PROGRAM FWC ft . 

c * ft ft r « V» 

c * ft SUBROUTINE COFN ft 

c * * * 0^°%^ ^ 

c * ° f 

e * 02/03/78 

e * 

c ft EVALUATION OF THE COORDINATES OF THE CENTERS 
c * GIVEN THE COORDINATES OF THE NODES 
e * 

cccnon/rol 13/a tar (A) 
concqn/gecoi/nbidsl.nblds.signia.finu.etsn (25) ,knnvr,«tai (6) ,knivr , 

61 twist, thotad (25) • theta (25) , thetec (21*) , thataQ, thetOd, alphas, 
6cd0,cdk,dp3ind,dpsin,dpsi id.dpsi i ,cocff ,coaff l,c,s, blades, 

6nnvr , nnva , nnvr l , nnva l , oner , nnca , stanv (25) , otanc (26) , 
6n!vr,niva,nivrl, nival, nicr,nica,otaiv(6) ,staic( 7 ) , 

6ntva,ntval,ntca, fpsl,fp32,nrl 1 ,nr 12, Ir 1 1, ir 12, irsodcl 
cessaon /resul/gcccc(26) ,twx(2i»,2 U) ,twy (21*. 2L) ,tws(2i»,2i») , 

6wxc(2l») ,wyc (2A) ,wzc (2L) ,wxnc (26,19) ,wync (26,19) .wane (26, 19) . 
$wslc(9,5l) ,wyic (9,50 ,wzic(9,51) ,wrnv(28,18) ,wtnv(28.l8) . 

6wznv(23. 18) ,wriv(8,50) ,wtiv(8,50) ,w2iv(3,50) , 

6xnc(26,19) ,ync(26,19) ,znc(26,l9) ,xic(9,50 ,yic(9.50 , 

6zic(9.50 , xnv (28. 18) ,ynv (28,18) ,znv{28,18) ,xiv(3,50) , 

6yiv(8,50) ,ziv(8.50) ,ia,ib 

. d icons ion rnvl ( 19 ) ,rnvr (19) ,r Jvl (51) ,rivr (51) ,dpnl (18) ,dpnr (18) , 

6dpi 1 ( 50 ) .dpi r ( 50 ) 

equivalence (nnvr.nncrl) , (nivr.nicr l) , (nnva, nnca 1) 
equivalence (niva.nical) , (nnvrl,nncr2) , (nlvrl,nicr2) 

c * 

c ft INNER POINTS 
c 

do 1 I"2,nncr2 
do 1 j“2,nncal 

xnc (i »j) ".25* (xnv(i-l,j-i)+xnv(?-l, j)+xnv(i ,j-l)+xnv(i ,j) ) 
ync (I ,j) ».25ft(ynv(i-l,j-l)+ynv(l-l,j)+ynv(i , j-l)+ynv(i ,j)) 

2 nc (i , j)“.25* ( 2 nv(i-l,j-l)+ 2 nv(l-l,j)+ 2 nv(i ,j-l)+znv(i , j) ) 

1 continue 

c * 

c ft CENTERS OUTSIDE THE WAKE (EXEPT CORNERS) 
c * 

do 3 l"2,nncr2 

xnc(i , 1) »xnv (i , l)+xnv (l-l, 1) -xnc (? ,2) 
ync (i , 1) ®ynv(i , l)+ynv (1*1,1) -ync (i ,2) 

2 nc (i , 1) - 2 nv (i , 1) +znv (i -1 , 1) -znc (i ,2) 

xnc ( i , nnca) «xnv ( i , nnva) +xnv ( i - 1 , nnva) -xnc ( i , nnca I) 

ync (i ,nnca) “ynv (i , nnva)+ynv(i-l ,nnva) -ync (i ,nncal) 

2 nc (i , nnca) -znv (i , nnva) + 2 nv (i -1 ,nnva) -znc (i , nnca 1) 

3 continue 

do 7 j "l, nnca! 

rnvl (j) -sqrt (xnv (1 , j) ft*2+ynv (l, j) **2) 
r nvr (j ) °sq r t (xnv (nnvr 1 , j ) **2+y nv (nnvr 1 , j ) **2) 
if (j .cq.nncal) go to 7 

dpnl (j) »atan2 ( (ynv (1 , j)+ynv (1 , j+1) ) , (xnv (1 , j)+xnv (1 , j+1) ) ) 

dpnr (j) -atan2 { (ynv (nnvr 1 , j) +ynv (nnvr 1 , j+l) ) , (xnv (nnvr 1 , j) +xnv (nnvr 1 , j-*- 



. / 0 





6 

c ft 
c * 
c ft 


c ft 

C ft 
C ft 



do 6 j “2, nnca l gSG&AL PA 6S & 

xnc ( 1 . j) ■ (rnvl (j)-trnvl (j-1) ) ftcos (dpni (j-1) ) -xnc (2, j) a^&lTV 

ync(l.j)*(rnvl(j)+rnvl(j-l))*sin(dpnl{j-l))-ync(2,j) C,r 
2 nc(l, J)“ 2 nv(l.j)+ 2 nv(l,j-l) -znc (2,j) 

xnc (nncr 1. j) ■ (rnvr (J)+rnvr (j-1) ) *cos (dpnr (J-1) ) -xnc (nncr2. j) 
ync (nncr l , j) - (rnvr (j)+rnvr (j-1) ) ft 3 in (dpnr (j-1) ) -ync (nncr 2, j) 
rnc (nncr 1 , j) »znv (nnvr l , j) +znv (nnvr l , j- 1) -znc (nncr2 , j) 


continue* .. ‘ ... 




CORNERS 


xnc (1, l)«* 4 .*xnv(l, t) -xnc (1.2) -xnc (2,1) -xnc (2.2) 

ync (1, l)»l*.ftynv(l, l) -ync (1,2) -ync (2,1) -ync (2,2) 

znc (1. l)"l».* 2 nv(l,l)-snc(l, 2 )- 2 nc (2.1) -me (2,2) 

xnc (nncr 1, 1) <*4*xnv (nnvr l, 1) -xnc (nncr 1,2) -xnc (nncr2, 1) -xnc (nncr 2, 2) 

ync (nncr 1 , l)* 4 ft ynv (nnvr l, 1) -ync (nncr 1,2) -ync (nncr 2, 1) -ync (nncr 2, 2) 

2 nc (nncr 1 , l) ® 4 * 2 nv (nnvr 1 , 1) -sne (nncr 1, 2) -znc (nncr 2, 1} - 2 nc (nncr 2, 2) 

xnc (l.nnca) -4.*xnv (l.nnvo) -xnc (2,nnca) -xnc (l.nncal) -xnc (2,nncal) 

ync (l.nnca) «4.ftynv (l.nnva) -ync (2, nnca) -ync (l.nncal) -ync (2, nnca l) 

me ( l.nnca) *4. *mv (1, nnva) - 2 nc (2,nnea) - 2 nc (l.nncal) - 2 nc (2,nneal) 

xnc (nncr l.nnca) «4.*xnv (nnvr l.nnva) -xnc (nncr2,r.nca) -xnc (nncr l.nncal) 

fi-xnc (nncr2.nncal) 

ync (nncr 1 . nnca) "4 . *ynv (nnvr 1 , nnva) -ync (nncr 2 , nnca) -ync (nncr 1 , nnca l) 

5- ync (nr.cr2,nnccl) 

2 nc (nncr f, nnca) -4 . * 2 nv (nnvr l , nnva) -znc (nncr 2 , nnca) -znc (nncr 1 , nnca 1) 

6- znc (nncr2,nncal) 


return 

end 





subroutine !oop22 



c 

c 


ft 

ft ftftttftrtftftftftftftftftftftftftftftftftft 

. n w ,n . v ft 

ft ' 'a ‘' PRQSRKftfyc ft 
ft ft * y ’ ft 

ft ft SUBROUTINE LOOP22 ft 
* ft . o 

ft ftftftftftftfetftftftftftftftftftftftftftft 


of&o** 

o* 


ft ,12/1/EO 
ft 

ft LOOP OH THE ELEMENTS, TO COMPUTE THE VELOCITIES 
ft AT THE CENTERS, OY VI 1103,11,1 

ft 
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ft 

ft 

* 

ft 

ft 

ft 

ft 


IF THE NUMBER OF 

ARE HOT USED 

BLADES 

HEAR WAKE 

IHT WAKE 

FAR WAKE 


ROLLED UP VORTICES IS 3 THEN IRL22.HRL22.CG2.Q2 
EL 

/ A / B / / / 

/ / O' / C / / 

/ / / / F / 


ccc23n/ron6/ql,<q2,q3,qJ*. !w 

cccsaon/al l/cgl,cg2,cg3,cg2*,Ir 

cosscn/ro) 1 10/wxctr2(25) ,wyctr2 (25) ,wzctr2 (25) 


cocson/partn/ivfr f Ird, I trace, Ipsea.reaas, I trcw, i trcll, I trc!2, t trcg. 

Si trcnt, i trcf . 1 plot, l cgen, I plotv, l trtg, Ivargr , 1 teat, I seneb, I sasict, 
Slir.s, l Ini, 1 icl.ni ter, I trsns, i pn, j punch, iplotw, i sw, icent, iplotg, 
Slpr, iplott 

ccicnen/gccjn/nbldsl.nblds.sigina.fmu.ctan (25) .knr.vr.etai (6) ,knivr, 
Cltwist.thctad (25) .theta (25) ,thetoc(24) .thetaO.thatOd, alphas, 
CcdO,cdlc,dpslRd,dp3in,dp3i id, dps i t ,cce? f ,cesf f l.c.s.bl edes, 
6nnvr,nnva,nnvrl,nnval,nncr,nnca,etanv (25) ,ctanc(26) , 

Snivr, niva.nivrl, nival, nier.nlca.ctalv (6) ,ctaic(7) , 
Cf5tva,ntvsl,n:c3,fpsl,fps2.nr I l,nrl2, ir 1 1, ir!2, Iccdal 
coszan /resul/ga^c (26) ,twx (24,24) ,twy(24.24) ,twz (24.24) , 

Swxc (24) ,wyc(24) ,wzc{24) ,wxnc(26,19) ,wync(2&,19) ,wznc (26.19) . 

Gwxic (9,51) ,wy ic (9,51) ,wzic(9.51) .wrnv(23.1S) ,wtnv(28.18) . 
6w2nv(23.i8) ,wr iv (3,50) ,wtiv(S.50) ,wziv(3.50) , 

Gxnc ( 26 , 19) ,ync(2&.l9) ,znc(26.19) ,xic (9.51) .yic(9.5D .xic(9.51) . 
Sxnv(28.l8) ,ynv(28,18) ,2nv(28. 18) .xlv(8.50> ,y iv (0.50) ,2iv(3.50) , 
Gia.ib 

ccraen/cvind/hnln (24,24) ,ha2n (24,24) ,hacantl(24) ,hamt2(2r») .ha.'iatS (21,) 
C02r»n/rol I l/hcrcat5 (26) .hc-tstG (26) 

ccs=on /val/x.y.z,ux,uy,u2,xi.yl,zi,x2,v2,z2.x3,y3.23.x4.y4.24. 

5gn,gal.gn2,dga. i t.rho.zpf ,tl ,epsl,cps2.stre,jca 
ccsson /vort/graln(24) ,ga2n(24) .grali (5) ,gra2i (5) ,” 

Sgascnt 1 , garat2 , garcnt3 , gasrantU 
cacion /savtw/twyt (24,24) , twit (24.24) 

CCEison/tipco/wxntl.wyntl.wzntl ,wxnt2,wynt2,wznt2,wxntp,wyntp,wzntp 

coaTiOn/sejni /de 1 2 

equivalence (nnvr.nncrl) , (nivr.nicrl) , (nnva.nncal) 
equivalence (ni va.nical) , (nnvr l,nncr2) , (ni vr l ,ni cr2) 
dirsension dgrai (25) ,dgm2(25) .dgm3(25) ,dgm4(25} . dg.*n5 05) ,dgm6(15) 

A 

* S2-SIH(15) .C2-C0SO5) ,TT2“2.*TG (15/2) 

data s2.c2.tt2/. 253319045.. 965S25326..26320L995/ 
data twopi/6. 283185306/ 

data Igln, lg2n, Igli ,lg2i/"GmN","GM2N","Grtll","CM2l"/ 



data lgora/"GAMC"/ 
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c . 

e Afti5?TAe’ftAT:ftftftftftftftftftftftftftftftftftftftftftAAAAAAftAAAftA 
C A ^ *• » . ' '• 

c * STRENGH OF THE ELEMENTS 

C A 

C TIP LOSS 

c 

gammc (nncrl)-O.O 
' , nnvr2-nnvr I-l 
' nncr3*nncr2-l 
d6 600 l-l,nnvrl 
grain (i)»0. 
gra2n(i)“0. 
do 600 j»l,nncr2 

grain (i) -grain (i)+hmln(i,j)ftgammc (j+1) 
gra2n (1) -gm2n(i)+hm2n (i , j) ftgammc (j+1) 
600 concinuo 

C A 
C A 



gcmmtl-O. 

• gsmmt2»0. 
do 604 j«l,nncr2 
gammtl-gammtl+hammtl (j) Aganrnc (j+1) 
gamrat2-gammt2+hamrat2 (j) Agammc (j+l) 

601* continue 

jgra=0 

if (nmes.ne.l) go to 723 

gtrain-0. 

gtraax-O. 

1 f (gtmiri.gt.gamnc (j+1)} gtmin-gammc (J+1) 

I f (gtmax. ! t. gammc (j+1) ) gtmax-gsmmc (j+l) 

If (gtrain.gt.gammc (j+1) ) jrain-J 
i f (gtnax. 1 1. gammc (j+1) ) jrnax-j . 

725 continue 

gtt-gtmax 

I f (-gtmin.gt.gtmax) gtt-gtmin 
jtrn-jraax 

If (-gtmin.gt.gtmax) jtm-jrain 
723 continue 

C A 

If (ltrcg.eq.0) goto 610 
wri te (iwr , 153) niter 

153 format (///, "aaa 100P2 aaa (ITERATION:", 13,")") 

wr I to (iwr , 1 10) I gam, (ganrnc (ki) ,k?-2,nncr 1) 
if (i trcl2.eq.0) goto 610 
write (Iwr, 110) I gin, (gmln(ki) ,Ui«l.nnvr 1) 
wr i te (iwr , 1 10) lg2n, (gm2n(ki) ,ki-l,nnvrl) 
write (iwr, 110) 1 g 1 i , (gml i (k i) ,k: -1 ,ni vr 1) 
wr I te (iwr , 110) 1 g2 i , (gra2 i (k i ) , k i -I , n i vr l) 
wr i te (iwr ,615) gamratl , gcmmt2,gammt3.gamratl* 

615 formatC GArtMT 1 TO l» ",Ael3.1») 

110 format (lx. 3k, lx.gf 10.6,/,6x,9f 10. 6) 

610 continue 

c A 

C AAAAAAAAAAAftAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
C A 

c A SAVE INFLUENCE COEFFICIENTS OF PRECEDENT ITERATION 

C A 


do 60 i-l,nncr2 



do 60 j«l,nncr2 
twyt(l ,j)»twy (i .J) 
tw2t(l.j) a :wzd.j) 
continue 
ft 

It *• ^ 

’T, ^7ppRl<»Sj^t , i. t » 

' /OpaS^piSV ? 

Itt2<>itrcl2 

If (ltt2.eq.M ItrcI2“2 
If (Itrcf.cq.l) Itrcl2®0 
continue 
call clv2 

I dcb«*5 

If (I trans.oq.O.or. I tref .cq.l) ldcb«l 
If (Itref .eq.J) goto 2g0 
If (Itrcnt.cq.l) call inicn2 
If (ltrcl2.no. 0) wrl tc(lwr,2C0) 
forcwtP'lftftA IOOP2 ftftft") 

* 

If (ltscdol.cq.3) Co to 111 

lf(lcodQl.cq.l«} go to 222 

0° to 333 

do kkk l«l,nnvrl 

do l)l»li jol.rnva 

2 nv(l ,J)iK).0 

go to 323 

do 2111 {■l.nnvrl 

do 2111 l*»l,nnva 

1 znv(i * j) °>znv(nr 12, J) 
do 2222 5»nrl t,nrl2 
do 2222 lM ( nnva 

2 znv(l ,J)«znv(nrl2,j) 
ft 

continue 




■^v 


ftftftftftftftftftftftftftftJVftftfc** 

* 

ft VELOCITIES INDUCED DY THE CIRCULATION OH THE BLADES 

* 

TIP LOSS FACTOR 

gcccc (nncrl)==0.0 
1 1-2 
Iw-1 

If (itrcl2.gt.l) wrlta(iwr,301) 

format (" VELOCITIES INDUCED BY THE CIRCULATION ON THE BLADES") 

xl*otanv (1) 

yl-O. 

zl«o. 

y3«0. 

a3“0. 

do 1 I»l,nncr2 
x3 n ctanv (i + i) 
gm"-ga»mc li + 1) 
cal l coord 
call vindnl2 (1 ,nnva) 
call vindi2 (2,niva) 





If (imodel ,eq.2) call vlndn2 (l.nnca) 69 

xl«x3 

I continue 

If (itrcl2.eq.2) call prtv 
If (itrc!2.cq.3) call prtvd2 

c * 

*' wxntl«=wxnc (nnvr 1,2) ORlGUi^l* PASS* (5 

■' ■ .wynt l“wync (nnvr 1 . 2) o* GUAUTV 

wrnt l“wznc (nnvr l , 2) 
wxq£2^j«nc (nncr 1 , 2) 
wy nt2«”wy nc (nncr l , 2) 
wznt2«w2nc (nncr 1 , 2) 

wxntp- (wxnfe<tyivr 1 , l) +wxnc (nncr t , 1) +wxnc (nnvr 1 , 2) +wxnc (nncr 1 , 2) ) A . 
wyntp* (wync (nnvr l, l)+wync (nncr l, l)+wync (nnvr 1,2) -5-wync (nncr 1,2} ) A, 
wzntp- (wznc (nnvr 1,1) +wzne (nncr 1,1) +wznc (nnvr 1 , 2) +wznc (nncr 1 , 2) ) A. 

c * 

c A VELOCITIES INDUCED BY THE NEAR WAKE 
ph I b-twopi /blades 
kapa«lf lx (phib/dpsli) 
nvsp»niva-kapa 

c * 

c A SURFACE ELEMENTS 
I t*»i 

If (itrcl2.gt.l) write (iwr, 302) 

302 format ("1VEL0CITI ES INOUCED BY THE HEAR WAKE (SURFACE ELEMENTS)") 
do 3 I a l,nnvr2 

gtnl«S!nln(i) 
gm2»sra2n (!) 
do 3 j°l,nnval 
xl»xnv(i , j) 
yl»ynv(i ,j) 

2 t a znv(i,j) 
x2®xnv(i+l,j) 
y2*>ynv(i+l,j) 
z2»znv (»+l,j) 
x3“xnv (I ,j+l) 
y3-ynv(i ,j+l) 

23 » 2 nv(i,j+l) 

x4»xnv(i+l,j+l) 

yl»»ynv(i + l,j+l) 

zl»-znv(i+l,j+t) 

call coord 

call vindb2(l,i) 

if (imodel ,cq. A) i t**8 

if (imodel .cq.4) call vindnl2 (2, nnva) 

if (imodel .cq.l«) it**l 

i f (imodel .eq. 2) call vindn2(l,nnca) 

3 continue 

if (i trcl2.cq.2) call prtv 
if (i trcl2.eq.3) call prtvd2 
if (i trcnt.cq. 1) call cxecn2 

c A 

c A SEGMENTS ELEMENTS 

c * 

it-2 

if (i trcl2.gt. 1) wr i tc (iwr , 303) 

303 format ("1VE10C1TIES INDUCED BY THE HEAR WAKE (SEG. ELEMENTS)") 
do 500 i“l , nnvr 1 , nnvr 2 

if ((nmcs.eq.2) .and. (i .eq. 1)) go to 500 




<*■ 


510 
500 
e 
e 
e 
e 
e 
e 


GR-ganantl 

If (I .eq.nnvr l) gm»ga.Ttnt2 
do 510 j»l,nnval 
xl»xnv(l ,j) 
yl»ynv (I , j) 
zl»xnv(l ,j) 

«3"*r>v(l ,j+l) 

•• y3«yrftf(h<J+l) 

•aS««wf(l i j+l) 
call coord 
call vlndb2(2,l) 

If (tRsdol.oq.il) it»9 

If (Ictadol .cq.lt) call vindnl2 (2,nnva) 

If (lcodol.cq.il) lt«2 

If (Inxodol .cq.2) call vlndn2 (l.nnca) 

continue 

continue 


***** ***************** 

* 
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.♦j 

o^ 0 o^ 


ft 

ft 

ft 


COHPUTE THE STRENGTH OF TIP. HIOOLE.AND ROOT VORTICES 
lt *»2 


1000 

2000 


c 

c 


5000 


liOOO 

c* 
c * 


nncr5“nncr2+3 
do 1000 j<°3.nncr2 
k«=nncr5”j 

If (aba (gzunrac (k) ) . I t.sbs (garrac (k-1) ) ) go to 1000 
go to 2000 
continue 
qli«-55T3nc (k) 

q3«SS3Sx (k) -gairac (k-3) 
q l-gcrrsnc (k-3) -garac (2) 

LOCATIOM OF SECOND HI DOLE ROLLED UP VORTEX 

ki»k-l 

k3-k-3 

kli«k-l| 

P“Q. 

pp»0. 

do LOOO loli.kli 

dgral (1+1)» (gacaoc (1+1) -garsnc (1)) *ctanv (1) 
p"p+dgn»l (l + l) 

pp«pp+ (gasrac (1+1) -gaisnc (l) ) 

continue 

cg2«p/pp 

LOCATION OF THE TIP ROLLED UP VORTEX 

pl»0. 

p2"0. 

do LOOO i“k,nncr2 

dgrn3 (i+l) « (gairwic (i) -garamc (1+1) ) *etanv (I) 

pl*pi+dgn3 (i+l) 

p2-p2+ (ga.Tsnc (i) -garrenc (i+l) ) 

continue 

cgli«pl/p2 


LOCATION 

b“0.0 


OF THE FIRST HI DOLE ROLLED UP VORTEX 
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c5*0-0 

do 7000 jolO.kl 

b»b+ (gommc (j+l) -gernmc (j) } *etanv(j) 
e5"c5+ (gsmme (j+l) -genvnc (j) ) 

7000 continue 

cg3 n b/c5 
c* 

eft ' LOCATION OF THE. ROOT ROLLED UP VORTEX of&KftAf. PA6S ‘ S 

ppi-o. n/pao^ Quatrrt 

pp2«0. 

ganssc (1) ,: 0.0 
do 8000 1-2, U.- , 

, ppl=ppl+ (gcrex (i+1) -gsnrnc (I) ) ftetanv(i) 

, pp2“pp2+(gsrrinic (i+1) -gsrenc (i)) , 

8003'/. 'continue 

cgl-ppi/pp2 

c * / 

wr i te (iwr ,1(001) egl ,cg2,cg3,egi»,ql ,q2,q3,qi* 

2*00 1 fornat(" location of rollup vortices - "#8fl0.5) 
if (itrcl2.gt. 1) wri to(iwr,307) 

307 format ("1VEL0C I TIES liJDUCEO BY THE IHT. WAKE (SEG. ELEMENTS)") 

irI2i*ir12-l 
nr!21»nr 12-1 
c nr!22-nr12-2 

c ir!22«irl2-2 

c VELOCITIES INOUCEO BY ANTICIPATED NEAR WAKE ROLLED UP VORTICES 

c 

do 7 i -nr 11, nr 12 
if (i .eq.nrl 1) grn*ql 
c if (i .eq.nrl 22) gw»q2 

if (i .eq.nrl21) gn«q3 
If (i .eq.nrl 2) gm-qi* 
do 6 j=*l»nnvai 
xl-xnv(i,j) 
yt-ynv(i ,j) 
zl-znv(i ,j) 
x3°xnv(i ,j+l) 
y3 B ynv U .j+l) 
z3»znv(i .j+l) 
call coord 
call vindnl2 (2,nnva) 
call vindi2 (2,niva) 

6 continue 

7 continue 
c 

c VELOCITIES INDUCED BY INTERMEDIATE WAKE ROLLED UP VORTICES 

do 18 i“irll,irl2 

if ((nmes.eq.2) .and. (i ,oq. 1)) go to 18 
If (i .eq. iri 1) gn-ql 
c if (i .eq. i r 1 22) gm»q2 

if (i .eq. i r 121) gm*>q3 
if (i .eq. ir 1 2) gm=*qA 
do 17 j-l, nival 
dps2“eps2 
xl-xiv (i ,j) 
yl*yiv(i . j) 
zl”z i v (i ,j) 
x3*»xi v (i .j+l) 
y3“yiv(i .j+l) 

c CORE BURSTING AFTER FIRST BLADE INTERSECTION 



rr 


phll**dp3inft(nnva-l)+dpsllft (j-1) 

If (phll.gt.phib) cps2 a .l 
call coord 
call vlndnl2 (2.nnva) 
call v!ndl2 (2,niva) 

If (imodel .eq.2) call vlrdn2 (l.nnca) 
cp32»dps2 

17 ; continue* 

10 , '^cjUlnuo"‘ 

If (ltrcli.sq.2) call prtv 
If (Itrcl2.cq.3) call prtvd2 
If (Itrcnt.eq.l) call e»ocn2 

ft 

ftftftftftftftftftftft*ftftftftft***ft**ftft*ft 
ft * - 

T VELOCITfE* INDUCED BY THE FAR WAKE 

ft 
ft 

continue 
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c 
c 
c 
c 
e 

c * 

250 

c ft 
c 
c* 
eft 


O* T 


FOR 2 SPIRALS : 720 DEGREES 


308 


c 

c 


205 


220 

c 

201 

c 


do 50 1-Irll,!rl2 
If (1 .oq.lr 121) gm«q3*l. 

If (l.eq.lrl l) gm**qlftl. 

If (l .cq. Irl2) grr.=qi|ftl. 

If (I .eq. I r ( 22) g:n»q2ftl. 

If (i trel2.gt.l) wri ca(iwr,303) l ,gn> 

fornat (" 1V„ L. INDUCED BY THE FAR WAKE (SEC) , l<=", 13." GK«".el3.4) 

dols“2iv (1 ,niva) -siv (I ,nvsp) 

zpf*»2lv (I ,niva) ft2-2iv(l ,nvap) 

rho®.5 ft sqrt ((xlv(t ,niva) - trIv (I .nival)) **2+ 

&(y!v(I ,niva)+ytv(l .nival)) **2) 

If (itref .na.l) goto 201 
ft£S£5S5£S£ &UUteUtetete&te&&te 
ft VERIFICATION OF THE MODEL OF THE FAR WAKE 
nk«!sS0 

do 220 ktt»l,2 
lt®4-ktt 

wrltc(lwr,205) It, I 

fornat ("1** VERIFICATION OF THE FAR WARE a** IT**”, 1 3.” 1 1- 4 '. 1 3) 

if (It.eq.2) goto 210 

call vlndnl (l.nnca) 

call vlndl (2,niva) 

call prtvd 

continue 

conti nuo 


. ! 


c ft 
c * 
210 
50 

c ft 


!t“3 

call vindnl2(l,nnva) 
cal) vlndl 2 (2,niva) 

I f (Icnodel .cq.2) call 


vl ndn2 (l.nnca) 


continue 

continue 

If (i trcl2.cq.2) call prtv 



73 


75 

3io 


838 

889 


If (Itrcl2.cq.3) call prtvdZ 
contl nua 

If (ltrcJ2.eq.l) writB(lwp,310) 

format (" VELOCITIES INDUCED AT THE CENTERS") 

If (itrcU.eq.l) call prtv 

1 f (l trcnt.eq .-fc^and. I trcf .eq .O.and. i trcl 2.ne.3) 

If (Itrcf .eq. Desoto 88S 

If (irrc^.cq.J) goto 889 

,f (I-tccl'2.eq.2) i trcl 2-1 

If (IttS'.'cq.i,) Itrcl2»4 

return 

contlnuo 

I trcf *>0 

Itrcl2“ltt2 • 

goto 999 

continue 

I trcl2*l 

goto 999 

end 


cal 1 cxcen2 


obkssnal pass b 

OF PC 0.1 QQALfFY 


subroutine or tv 
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************ ft* ******** 

* ft 

ft PROGRAM FWC ft 

ft * 

* SUOROUTIME PRTV ft 

ft ft 

**ftftft*ftft*ft*ft**ft******* 






* 02 / 03/78 

* 

*•• PRINt 'OUT ^OR VERIFICATION 
ft ENTRIES: ' ’• 

ft PRTV: VELOCITIES AT THE CENTERS AND INFLUENCE COEFFICIENTS 
ft PRTCEN: POSITIGNS OF THE CENTERS 
ft FRTNOD : POSITIONS OF THE NODES 

* PRTVO VELOCITIES AT THE CEHTERS, 
ft THE CLEAR THE CORRESPONDING ARRAYS 

ft CLV CLEAR THE ARRAYS FOR THE INDUCED VELOCITIES 

ft AND THE INFLUENCE COEFFICIENTS 

ft PP.TVH VELOCITIES AT THE CENTERS OF THE NEAR WAKE 

* PRTV I VELOCITIES AT THE CENTERS OF THE INTERMEDIATE WAKE 

ft 

ccssnon/roll2/v«ivr(8#50) ,wxlvr (8,50) ,wyivr(8,50) .wznvr (28. 18) .wxnvr (28. 18) . 
Swynvr (28,18) 

common/ral 15 /wxctr (25) .wyctr (25) .wzetr (25) 

common/parra/iwr, I rd, I trace, Ipsem.nmas, I trcw, itrcl l, I trcI2, 1 trcg. 

Si trcnt, I trcf , I plot, icgen, iplotv, i trtg, ivergr.f test, isaneb, Isamot, 

• 61 Ies. 1 ini. 1 in2.nl ter, i t.rana, ipn, jpuneh, iplocw, isw, icont, Iplatg, 

Si pr, iplott 

co:maon/goc: 3 /nbld 5 l,nblds, 3 lsma,fmu,etan( 25 ) .innvr.etal (6) .krtivr, 

&1 twist, thctad (25) .theta (25) , thatoc (24) , thcazO.thotOd, alphas, 
ScdO,cdU,dpsind,dpsln,dpsi id.dpsi i ,coeff,coeffl,c, 3, blades, 

Snnvr , nnva, nnvr 1 , nnva 1 , nncr , nnca , e tanv (25) . ctanc (26) » 

Snivr .niva.nivr l.nival.n'er ,nica,ctaiv(6) ,ata»c(7) , 
Sntva,ntval,ntcs,fpsl,fp32,nrll,nrl2, irl l,ir!2, I model 
common /rosul/gammc (26) ,twx(24,24) ,twy (24.24) ,twz (24,24) , 
twxc (24) ,wyc (24) »wzc (24) ,wxnc (26, 19) ,wync (26.19) .wznc (26.19) • 

Swxic(9.5i) ,wyic(9,51) ,wzic(9,5l) ,wrnv (28, 13) .wtnv(28, 18) . 

Swznv(28.13) ,wriv(8,50) ,wtiv(0,5O) ,wziv(8,50J , 

Sxnc (26.19) ,ync(26.19) ,znc (26,19) ,xic(9.5U.yic (9.5D .zlc 0.51) , 

6xnv(28. 18) ,ynv(28.18) ,znv(28.t8) ,xiv(3,50) ,yiv(8,50) ,z!v{R,50) , 

Sia.Ib 

dimension wr (26) 

equivalence (nnvr, nncrl) , (nivr.nlcrl) , (nnva, Racist) 
equivalence (niva.nical) , (nnvr l,nncr2) , (nivr3,nicr2) 
data Iwxnc, Iwync, Iwznc, 1wr/"WXNC","WYNC","W2NC' , ,"WR "/ 
data lwxic.lv/yic, Iwzic /"WX I C" , "WY I C" , "WZ I C"/ 

data ltwx,ltwy,ltwz/"TWX ","TWY ","TWZ "/ 
data Ixnc, lync, !znc/"XNC ","YNC ","ZNC "/ 
data lxic. lyic, lzic/"XIC ","YtC ","ZIC "/ 
data Ixnv, lynv, lznv/"XNV V'YNV ","ZNV '*/ 
data Ixiv, ly iv, lziv/"XIV ","YIV ","ZIV "/ 
ft. 

w (x.y.wx.wy) - (wx*x+wy*y) /sqrt (x*ft2+y**2) 

ft 

format (lx,a4, lx,6f 10.6, /,6x,6f 10.6./, 6x,6f 10. £) 
format (lx) 




lde!°0 

write (iwr , 100) 

100 format (" *** PRTV ***") 

c ft 

60 continue 
c * 

do 1 JM.nnca 
do It i**l,nncr 

wr (i)> s w(xnc(i,j),ync(J,J),wxnc(i,j),wync(i,j)) 
it ' continue 

' write (iwr, 102) 

write (iwr , 101) Iwxnc, (wxnc (i , j) , i«*l ,nncr) 
wr i to (iwr , 101) Iwync, (wync (i , j) , 1®1 ,nncr) 
write (iwr , 10 1) Iwznc, (wzne (i , j) , 1*1 ,nncr) 
write(iwr.lOl) Iwr, (wr(i) ,i-l,nncr) 

1 continue 

if (idol. eq. 2) return 

c * 

61 continue 

do 2 j°l,n?ca 
do 5 i“l.nicr 

wr (i) ®w(xic (i ,j) ,yic(i,J) ,wxie(l,J) ,wyic(l,j)) 
5 continue 

wr i te (iwr, 102) 

wr i te (iwr , 101) Iwxic, (wxic (i ,j) , ?«1 ,nier) 
write (iwr, 101) lwyic, (wyic(5 ,j) , i«i,nicr) 
wr I te (iwr , 101) Iw 2 ic, (w 2 ic(i,j) ,i“l,nicr) 
write (iwr, 101) Iwr, (wr (i) , i = l,nicr) 

2 continue 

1 f (idol .eq.3) return 

c * 

if (idcl .cq. 1) goto 70 

c * 

c ft PRINT OUT OF THE INFLUENCE COEFFICIENTS 
c * 

do 3 J**l.nncr2 
wr i te (iwr , 102) 

wr i tc (iwr , 101) 1 twx, (twx (i , j) , i"l ,nncr2) 
wri te (iwr , 101) Itwy, (twy(i.j) ,i«l,nncr2) 
write (iwr, 101) Itwz, (twz(i,j) ,i-l,nncr2) 

3 continue 
return 

c * 
c * 

entry prtcen2 

c * 

c ft PRINT THE POSITIONS OF THE CENTERS 
c ft 

wr i te (iwr , 103 ) 

103 format ( ! * *** PRTCEN ftftft") 

do 30 j*I,nnca 
wr i te (iwr, 102) 

write (iwr, 101) Ixnc, (xnc (i ,j) , i«l,nncr) 
wr ite(iwr, 101) lync, (ync(i.j) ,i-l,nncr) 
write(iwr, 101) lznc, (znc (i ,j) , i-l.nncr) 
continue 
do 32 j“!,nica 
wr i te (iwr, 102) 

write (iwr ,101) Ixic, (xic (i ,j) , i»l,nicr) 
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wrlto(iwr,101) lyic, (y 1 c (t ,j) , i»l,nicr) 
wr i te (Iwr , 101) l 2 ic, (zic (I , j) , i*l,nicr) 
32 continue 

return 

c ft 

C ft 


entry prtnod2 


c * 

c ft 

e 

y n ♦ . 

102 . ’ 


50 


52 


pniirr 4 ^* positions of the nodes 

% W 

;/<S if/ • 

. ,-wf ite(iwr,l02,) 

forma t(" ftft* PRTHOD ftftft") 
do 50 J«*l,nnva 
wr I to (iwr, 102) 

writo(iwr,101) lxnv, (xnv(l,J) ,i*»l,nnvr) 

write (iwr, 101) lynv, (ynv(i ,J) , t*»i P nnvr) 

write (iwr, 101) Iznv, (2nv(i ,j) , l«*l,nnvr) 

continue 

do 52 J**1 ,rIvo 

writs (iwr, 102) 

write (Iwr, 101) Ixlv, (xlv(i,j) ,i»I,nivr) 
write (iwr, 101) lyiv, (yiv(i.J) ,i«»l,nivr) 
wr i to (iwr, 101) l2iv, (ziv(i.j) ,i«*l,nivr) 
continue 
return 








oS 


eft 

C ftAftftirAftAfr&ftrtAArtTirtAAirAftiiAAAirAftAAiSftrtilAftiiAAAAAAAftAAAftjYAnftft 

C * 


entry prtvd2 

c * 

c * PRINT THE VELOCITIES AND CLEAR THE CORRESPONDING ARRAYS 
eft 

write (iwr, 105) 

105 formatt" ***« PRTVO ft*ft") 

ldal«*i 
goto 60 

c * 
c * 

c * 

entry clv2 

c * 

C * CLEAR THE ARRAY5 OF INDUCED VELOCITIES 
c ft 

70 continue 

do 80 i*»l,nncr 
do 80 j"l,nnca 
.wxnc(l,J)«0. 
wync (i , j) ®0. 
wznc (i , j) ®0. 

80 continue 
C * 

do 81 i"l,nicr 
do 81 j*»l,nica 
wxic(i ,j) “0. 
wyic(i,j)-0. 
wz I c ( i , j ) «0 . 

81 conti nue 
c * 




widft U iit * t«V,. Ji.h. If 
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do 82 i«i.nncr2 
do 82 j»l,nncr2 
twx (I ,j) “0. 
twy(i,j)-0. 
twr(l,j)“0. 

82 continue 

d©i33>,|“l .nncr 2 on.u 

wxcl^ fc. } oQ. QF (^JAUTV 

wyctr (i) «*Q. 
w2Ctr(i)“0. 

83 continue 

do 84 i*ir 1 1, !r 12 
do 85 j“l,niva 
wyivr (i , j) =0. 
wzlvr (I ,j)«0. 
wxlvr {I 

85 continue 

84 continue 

do 5017 il“nrll,nrt2 
do 5017 j l*t.nnva 
wxnvr(il,jl)“0. 
wynvr (i 1 • j 1) “0 . 
vanvr (I l,jl)«0. 

5017 continue 

return 

c * 

C ftftftft*ftftftft*ft*ftAft*ftft*****ftftaftftft*ftftft*ft*ftrtft**ft*ftft*aft 

C * 

entry prtvn2 

c * 

c * PRINT THE VELOCITIES AT THE CENTERS OF THE HEAR WAKE ONLY 
c * 

ldel“2 

writn(iwr, 106) 

106 formatC *** PRTVM ft**") 
goto 60 

c * 

C ftttftftftftftftftftftftftftftaftftftftftftftftftftftftrtftftrtftftftftftftftftftftftftftftftft 

c * 

entry prtvi2 

c ft 

c * PRINT THE VELOCITIES AT THE CENTERS OF THE INTERMEDIATE WAKE ONLY 
c * 

idel«3 

write(iwr, 107) 

107 formatC *** PRTVI ftftft") 
goto 61 

c ft 

end 


a 11 



c 

c 

c 

e 

c 

c 

c 

c 

c 

e 

e 

e 

c 

c 

c 

c 

c 

e 

c 

c 

c 


subroutine vind2 

ft ftftftftftftftftftftftftftftftftftftftftft 

ft ft ft 

ft ft PROGRAM FWC ft 

ft ft ft 

ft * SUBROUTINE VIN02 * 
ft ft ft 

ft_ ftftftftftftftftftftftftftftftftftftftftft 
ft ' ‘• 4 -• u 

ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
* 
ft 
ft 
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07/14/80 

LOOP OH THE POSITIONS 
ARE COMPUTED 
ENTRIES : VIHDB 


WHERE THE IHOUCEO VELOCITIES 


: 


* 

ft 


VIHDM 

VIMDI 

EXECM 

I men 


: 


VELOCITIES INDUCED ON THE BLADES 
(EVALUATION OF THE INFLUENCE COEFFICIENTS) 
VELOCITIES INDUCED OH THE BLADES HEAR WAKE 
VELOCITIES INDUCED ON THE BLADES 
INTERMEDIATE WAKE 

VERI F I CAT I CM OF THE INFLUENCE COEFFICIENTS 
INITIALISATION OF THE ARRAYS FOR 
VERIFICATION OF THE INFLUENCE COEFFICIENTS 


CMEBon/rol16/ql,q2,q3 t qi*, iw 

common/rol 15/wxctr (25) ,wyctr (25) .wzetr (25) 

ccmmon/roi 1 l/hamratS (26) ,h£imit6 (26) 

common/ro II 2/wz i vr (3.50) .wxivr (8.50) ,wyivr(8,50) »wznvr (23.18) * 

6wxnvr (28. 18) .wynvr (23. 18) 

ccmncn/cvi nd/ha In (24.24) ,hss2n (24.24) ,hamati(24) ,hcsat2 (24) ,hsmat3(24) 
common/parm/iwr, ird, i trace, lpsera.nraes, i trew, i trcl 1. • trc!2, i treg, 

61 trent, I tref , iplot, iegen, iplotv, itrtg, i verge, I test, issracb, issraet, 

61 ins, I irat, 1 ici2,ni ter, i trans, ipn.jpuncn, ipiotw, isw, icont, Ipiotg, 
6lpr, Iplott 

cc.TO!on/gcca/nbldsl,nblds,3igaa,fir.u,etan(25) .knnvr.etai (6) .knivr, 

61 twist, thetad (25) , theta (25) ,thetac(24) , thctaO, thetOd, alphas, 
6cd0,cdk,dpsind,dps in.dpsi id.dpsi i .cocff .coeffl ,c,s. blades, 

6nnvr ,nnva,nnvr! .nnval.nncr ,nnca.etanv{25) ,etanc{26) . 
6nivr,niva,nivri,nival,nicr,nica,etaiv(6) , ctaic ( 7 ) , 
6ntva,ntval,ntca,fpsl,fps2.nr 1 l,nr!2,irll,irl2, i medal 
carmen /r osu 1 /gansac (26) . twx (24 . 24) , twy (24 . 24) . twz (24 , 24) . 

6wxc (24) ,wyc (24) .wzc (24) ,wxne (26. 19) .wync (26. 19) ,w2nc (26. 19) . 

6wxic (9,51) ,wyic (9,51) .wzie(9.51) ,wrnv(28.18) ,wtnv(28.18) . 

6wznv (28. 18) ,wriv(8,50) ,wtiv(8,50) ,wziv(3.50) , 

6xnc (26. 19) ,ync(26. 19) .znc (26. 19) ,*ie (9.51) ,yic (9,51) ,zic (9,51) , 
6xnv (28.18) ,ynv(28, 18) ,znv(28, 18) ,xiv (8.50) ,yiv (3.50) ,2lv(8.50) , 

6ia, ib 

common /vol/x,y,2.ux,uy,u2,xl,yl,zl,x2,y2.z2,x3.y3.23.x4,y4.z4. 

6gra,gml,g’n2,dcia, i t, rho,2of . tl ,epsl,cpa2,stre, jgm 
common /vort/galn (24) ,gm2n(24) .gmli (5) ,gm21 (5) , 

6gammt 1 , gammt2 , gsmmt3 , gamatA 
dimension wzr (1,6) 

dimension wxcc (24) ,wycc (24) .W2cc (24) 
dimension wxcea (2A) ,wycca{24) .wzcca (24) 
common/rol lup/cg, icg.p.q 

equivalence (nnvr.nncrl) , (nivr.nicrl) { (nnva.nncalj 
equiva I ence (ni va.nical) , (nnvr l ,nncr2) , (n i vr 1 ,ni cr2) 
data lxc, iyc, lzc/"WXC "."UYC ","WZC "/ 
data Ixca, lyca, lzca/ l ‘WXCA' , , ,, WYCA" ,"WZCA"/ 


entry vindb2 (k, ikk) 


riVV,.^ A L.V.Ia'L., 




y*- v -Vi 



I • 


c * 

C ft 
C ft 
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./ 


CALCULATION OF THE INFLUENCE COEFFICIENTS 



e 

c 

c 


11 


e 


c 

c 


c 

c 


c 


3 

c 

c 

c 

c 

12 


* 

ft 

ft 


ft 


ft 

ft 


ft 

ft 


* 


* 


1 r 1 21-i r 12-1 

lr!22-irl2-2 

savl«gm 

sov2»dgnj 

nbl»nblds 

tf(it.eq.3) nbl*l 

If (k.eq.k) go to 101 
nncr3*nncr2-l 
do 1 j«*l,nncr2 
goto (11. 12.13). k 


ORIGINAL PAQ3 IS 
OF POOR QUALITY 


NEAR WAKE. SURFACE ELEMENTS 
continue 

If (hnln(ikk.j) .eq.0.and.hm2n(ikk,j) .eq.O) goto 1 

g:ss“.5ft (h?ain(ii(k,j)+hm 2 n(lkk,j)) 

dcu»-.5ft(hni2n(ikk,j)-hnln(ikk t j)) 

gn«gaft.5 
dgw»dgnft.5 
do 3 kk»l,nncr2 
do 3 Jj-1.2 
x«xnc (kk+l.jj) 
y-ync(kk-H.jj) 
z«anc (kk+l.jj) 
cc«l. 
ss“0. 

do 3 n-l,nti 
call wxyz 


ROTATION OF THE VELOCITY VECTOR OF' - (N-l) »2.*PI/BLA0ES 
twx (kk , j) “txx (kk , j)+uxftcc+uy*ss 
twy (kk , j) »twy (kk , j ) +uy*cc-ux*S3 
twz (kk.j)-twz(kk, j)+uz 

ROTATION OF THE POINT OF +2 .*PI /BLADES 
sav**x 

X"sav*C“yftj 

y-sftsav+yftc 

INCREMENT CCSS OF 2.*PI /BLADES 
sav**cc 

CC-'»CC*.C-SS*S 

ss»isftc+sav*s 

continue 


goto l 

ft 

* NEAR WAKE, SEGMENT ELEMENTS 

ft 


continue 

if (ikk.eq.l.ar.d.hammtl (j) .ftq.O.) goto 1 
if (ikk.cq.nnvr 1 .and.Ham.nt2 (j) .eq.O.) goto l 
gm*>ha.T.-nt 1 lj) 

i f ( i kk . cq . nnvr 1) gn=«hamrat?. (j ) 
goto 20 


c 

c 


ft 

ft 



00 


\ 


13 





continue 

If (hnll (ikk. j) .eq.O.and. hn2l (ikk, j) .cq.O) goto t 
Cn*.5*(hsti (;kk.j)fhrj2i (ikk.j)) 
dg3!*.5* (h=»2l (ikk,j) -hrali (ikk.j)) 
goto 20 


continua 

do k kk*l,nncr2 

x*M»tanc (kk+l) 

y»0. 

z*0. 

cc*l.' * 

*s°0. 

do k n*l,nbl 
call wxyz 

twx (kk, j) «twx (kU,J)+uxrtcc+uy*8« 
twy (kk,j) *twy (kk, j)+uyftcc-ux*83 
twa (kk, j) «>twz (kk,j)-K)2 




3 * 


*av*x 

x»sav*c-y*s 

y*s*sav+y6c 

sav*cc 

CC*CC*e-8S«S 
S8 0 38*c+sav*8 
k continue 

l continua 

ga^savi 
ds»“8ov2 
go to 1C2 

c * 

C * ROLLED UP VORTICES 
c * 

101 gc**ql 

If (ikk. oq. tr ) 2) 
if (ikk.cq. iri2t) gci*q3 

C if (ikk.cq.ir 122) sa»q2 

do 4k kk a l,nncr2 
x-otsnc (kk+l) 
y*0. 

2 * 0 . 

cc«i. 

*5*0. 

do kk n*l,nb) 
call wxyz 

wxetr (kk) *wxctr (kk) +ux*cc+uy*ss 
wyctr (kk) «*wyctr (kk) +uy*cc-UX*8S 
wsctr (kfc)-wzctr (kk)+uz 
aav-x 


x«3avftc-y*s 

y*3*sav+y*c 

sav**cc 

CC«CCOC”SS»S 


S5"3S<tC+SaV*S 

kk continue 

gal*: av I 
dgm»sav2 

i f (I trcnt.eq.O) return 
102 i f (i trcnt.ea.C) return 



81 


c * 

C rtaftftrtftftftftftftJtftftftai’tftas'lftaaaartftftftftaViftaftAaftrtrt 

C ft 

c * CONTROL: VERIFICATION OF THE INFLUENCE COEFFICIENTS 
c ft 

I f (k.na, l) goto 302 
do 301 kk“l,nncr2 
do 30J Jj«i,2 
x»xnc (kk+l.jj) 
y*ync (kk+1 , j j) 

2*znc (kk+l.jj) 
cc*»l. 

SS°0. 

do 301 n“l,nbl 
call wxyz 

wxcc (kk) «wxcc (kk) + (uxftcc+tjy*ss) ft.5 
wycc (kk) «wycc (kk)-* (uy*cc-uxftss) ft. 5 
wxcc (kk) *w2cc (kk) +U2ft.5 
sav»x 

x«savftc-yft# 
y»#ftsav+yftc 
*av«cc 

cc"cc*c-ss*s 
ss"ssftc+sav*s 

301 continue 
return 

c * 

302 continue 
do 304 Vk»i,nncr2 
x-etanc (kk+1) 
y«0. 

2 - 0 . 
cc-1. 
ss»0. 

do 304 n«l f nbl 
call wxyz 

wxcc (kk) -wxcc (kk) + (uxftcc+uyi's 3) 
wycc (kk) =wycc (kk) + (uyftcc-uxftss) 
wzcc (kk) «wzcc (kk) +uz 
aav-x 

x-3av*c-y«s 
y-sftsav+yftc 
sav»cc 


* « ffl 


CC"CCftC“S*ftS . 

*S°35ftC+30V*S 

304 continue 

c ft 

return 

c ft 
c * 

C ft 

entry vindn2 (kl,k2) 

C ft 

c * VELOCITIES INDUCED ON THE NEAR WAKE CONTINOUS VORTEX SHEET 
c ft 

nbi-n&ldo 
if (i t.eq.3) nbl*l 
do 15 i l"l,nncrl 
do 15 j l=k 1 ,k2 



\ 


15 


e 

e 

c 

c 


3000 


5011 

3001 


5010 

c * 
c * 
c * 


C2 












A 

A 

A 

A 


C 

c 


A 

A 


x°xnc(i I,jt) 
y^ncd l,jl) 
z*»znc(i l,jl) 
ce«*l. 

83*0. 

do 15 n«»l.nbl 
call wxyz 

wxnc (II , j l) ®wxnc (1 1 ,J 1) +ux*cc+uyAsa 
. ' wync (I l,j 1) «=wync (i 1, j lJ+uyAcc-ux^sa 
w2nc(l,l,jl)*wznc(i l,jl)+u2 

8BV°X 

X«savi>c-y*3 
y»SA3SV^7Ac 
S8V—CC 

CO»Ce*C-SB*3 

sf*3sAc+aav*3 

continue 

return 

VELOCITIES INDUCED AT THE WOOES OF THE ANTICIPATED HEAR HAKE ROLLED UP 
VORTICES 

entry v *ndnl2 (kl.k2) 
nbl-nbldsi 
ifOt.aq.3) nbl-i 
lp°0 

do 5010 k3”nrl l,nr!2 

do 5CI0 jl«2,nnva 

x-xnv (k3,jl) 

yynv(k3,j!) 

z»znv{k3,jl) 

cc-i. 

ea-O. 

Ip-0 

If (xi.eq.x.or.x3.eq.x) ip«2 
If <yl.«q.y.or.y3.cq.y) lp«2 
do 5010 n-l,nbl 
l f (I t.cq.3) go to 3000 
If (n.oq. l.and.Ip.eq.2) go to 3001 
cal 1 wxyz 

If (lt.cq.8.or.it.eq.9) go to 5011 
wxnvr (k3, j l) -wxnvr (k3, j O+ux^ccKiyAas 
wynvr (k3 . j 1} -wynvr (U3, j l) -Hjy^cc-uxAaa 
wznvr (k3.j U-wsnvr (k3, j l)+uz 
aav«*x 

x-3av*c-yAa 

y-aAaav+yAc 

sav**cc 

cc-cc*c-sa*s 

ss-asAc+savtfs 

continue 

return 

VELOCITIES INOUCEO AT THE HODES OF INTERMEDIATE WAKE ROLLED UP VORTICES 
entry vindi2(kl,k2) 


I 


nbl-r.blds 

if (i t.eq.3) nbl-l 
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do 5020 k3«irll,irl2 
do 5020 j“2,nlva 
x»xiv(k3,j) 
y»ytv(k 3 .j) 

2=zlv{k3,j) 
cc**l. 
ss“0. 
ip«*Q 

l f (xl.eq.x.or.x3.eq.x) ip*»l 
do 5021 n«*l,nbl 
if(it.cq. 3 ) go to 3002 
tf (n.eq.l.and. ip.eq. 1) go to 3003 
call wxy2 

wxlvr (k3.j)°wxivp (k3. j)+ux*cc+uy*ss 
wy ivp (k3. j) “wy ivr (k3. j)+uy*cc-ux*ss 
wzivr (k3*j) <B wzivr (k3.j)+uz 
sav"x 

x»sav*c-yfts 
y»s*sav+y*c 

SSV-CC 

cc«*ee*e-ss*s 

S£“SB*e+savfts 

continue 
continue 
return 

t ********************************** ft* 

c * 

c * VERIFICATION OF THE INFLUENCE COEFFICIENTS 
c * CALCULATION OF THE VELOCITIES BY THE INFLUENCE COEFFICIENTS 
c * AND HR I NT SUT 
c * 

entry cxecn2 

c * 

do 1»0I* I"l,nncr2 
wxcca (i) B 0. 
wycca (i) E 0. 

W2CC3 (i) “0. 
do 1*01* j®l,nncr2 

wxcca (i) »wxcca (i)+twx(i , j) ftgairenc (j+1) 
wycca (i) *>wycca (i)+twy (i , j) *gammc (j+1) 
wzcca (i) *wzcca (i)t-twz (i , j) *gammc (j+l) 

1*01* continue 

wr i to (iwr , 100) 

IOO format (" EXECN *ftft") 

wr I te (iwr , 1 10) Ixc , (wxcc (ki) ,ki»l,nncr2) 
wr i te ( i vrr , 1 10) lyc , (wycc (ki) ,ki»l,nncr2) 
wr i te (iwr , 1 10) I zc , (w2cc (ki) , kl«l ,nncr2) 
wr i te (iwr , 1 10) I xca, (wxcca (ki) ,ki«l.nncr2) 
write (iwr, 110) lyca, (wycca (ki) ,ki»l ,nncr2) 
wr i te (iwr , 1 10) Izca, (wzcca (ki) ,ki“l ,nncr2) 
return 

c * 

c ****** **** ft** ft** *********** ft *** A Aft* ****** 

C * 

entry inicn2 

c * 

do L05 i ■ 1 ,nncr2 
wxcc (i) *0. 
wycc (i) E 0. 




3002 

3003 

5021 

5020 

c * 



34 


. *y»r •*rv: • 


405 

e ft 
110 
c ft 


r r 


W2CC{!}»0. 

continue 

return 

forisat(ix,a4.1x.6f 10.6,/,6x,6f 10.6,/,6x,6f 10.6) 

end 




C:, y rr '-':o 

*0 




0? 







ii 



c 

c 


c 

c 


e 

c 


e 


c 


c 


c 

c 

c 


* 

* 


* 

ft 


ft 

ft 


* 


ft 


ft 


* 

ft 

ft 
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subroutine indvel 
real*A !0 

common /vel/x,y,z,ux,uy,uz.xl.yl,2l,x2,y2,22,x3.y3,z3,xA.yA,2A. 
figm.gml ,gm2,dgm, it.rho.zpf ,tl ,ep3t,eps2,stro, jgm 
common /farc/eta (AO) ,wx (22, AO) ,wy (22, AO) ,wz (22, AO) 
common/parm/iwr , ird, i trace, lpsem.nmcs, i trcw, i trcl 1 , I trc! 2, i trcg, 
61 trcnt, i trcf , iplot. icgen, iplotv, i trtg, I vergr , i test, isameb, i same t, 
61 ims, I iml, 1 im2,ni ter, 1 trans, ipn, j punch, I plotw, I sw, icont, iplotg, 
£ipr,iplott 
c oxmon/ s em i /d o l 2 


• data pf/ 3 . 1 A 1592653 / 
fpi : value of l./(A.*pi) 

•'••• data fpl, twopi/.079577A72, 6 . 283105306/ 

data aa,bb,cc,dd/5.12160o-2, 11,0501,-6.769760-7. .527592/ 


OR5SSSM. PAGE 13 
0 ? pc-na QUAbT? 


* . 

entry coord 

save 2mod2,xct,yct,2ct,dx,dy,xxl ,yy 1 , 22 I,xx2,yy2, 222 , xx3,yy3. 223 , 
6 rmax 2 ,ds,ds 2 ,ds 2 ,day,dsx, fva,epsl 2 
If (I t.cq. 2 ) go to 50 


surface element 
S®*«5 A (sml+gm2) 
dgffi*.5ft (gm2-gml) 
epsl 2 ". 25 *op 3 lftft 2 


agral-abs (gnl) 
agra2*»abs (gm2) 
xct“ (xl+x2+x3+xA) *.25, 
yct“ (y l+y2+y3+yA) *.25 
ZCt- (2 1+Z2+Z3+ZA) ft . 25 
first unit vector 

xxl-agmlft(x3-xl)+agm2ft (xA-x2) 
yyl«*B5mlft(y3-yl)+agm2ft(yA-y2) 

22 l»agm 1 ft (23-2 1) +agm2ft (zA-z2) 

2 mod l^sqr t (xx l*xx 1+yy l*yy 1 + 221 * 22 1) 

t2l»l./zmodl 

xxl»xxl*tzl 

yyl-yyl*t 2 l 

ZZ 1“22 lfttZ 1 

second unit vector 
xx2«* (x2+xA-xl-x3) 
yy2»(y2+yA-yl-y3) 
zz2-(z2+zA-2t-z3) 

2mod2 o sqrt (xx2*xx2+yy2ftyy2+zz2ftz22) 
find v2 corrected such as vl.v2“0 
we have v»aftvl+bftv2 

/vi/**l t /v2/“l.,/v/«*l, v.v2.gt.0 v.vl-0 
dot*xxI*xx2+yy l ftyy 2+22 1*222 
b"l ,/sqrt (2mod2ft2mod2-dotftdot) 
a—bftdot 
xx2«aftxxl+bftxx2 
yy2 n aftyy l+b»yy2 
Z22-a*Z2l+bft222 
22 3“ xx lftyy2— yyl*xx2 
xx3°yy I>:2z2-22lftyy2 
yy3 K 2zlftxx2-xxl*z22 


c * 

c ft 


dimensions of the elements 
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zmod2 a zxnod2ft.5 
dx o ZR0dl*.5/ Csgn»l*^fi3rn2) 
dy“zmod2 ft .5 

ds-zmod l*2!Bod2/ (asnl+agtn2) 

naax2 a i*Q0.*d3 

return 

c • * 

50 continue 

xct“.5 ft Ul + *3) *. 
yet«.5 < ’ t (yi + y3) ' /».£ 

2 Ct».5 ft (2I+23) fiY 

■■ i«,iSO<*SCt • . ' 

••■■ i.v •(! ;-ii-.r.lpeqrt(*I*xWTrl*yl) 
r3«tqrt (x3 ft «3+y3*y3) 
r0«.5 ft (r l + r-3) 
d8X».5 ft (x3*?£l) 
day*».5*(y3 - Vl) 

t d32«.£*(23*2l) 

’ ,ds2»d3X**Z+d5yft*2+dsz**2 
fvc®cp52*epf2ftds2 
rc , 4J«2“^C0.*ds2 
return 

c * 
c * 

entry wxyz 

c * 
c * 

lpr»6 

IfOt.eq.iO go to 88 
if (I t.oq.3) go to 77 
dxx»x-xct 
dyr»y-yct 
d22«2-2Ct 

r02®dxx**2+dyy* r *2+ , d2z**2 
if (l t.eq.2) go to 2 
if (it.cq.3) go to 2 

flC»cn/2T!ad2 

gv«dgra/ (2sod2*dy) 

xp«xx l*dxx+yy l*dyy+zz l*dz2 

yp=xx2*dxx+yy2*dyy^222*d22 

zp»xx3 A dxx+yy3*dyy+223 ft d22 

if (r02.lt. rnax2) go to J» 

r03“r02*sqrt (r02) 

f«ds**fpi*gc/r03 

wyp«-f*2p 

Wzpaf ftyp 

go to 10 
li contlnuo 

c * 
c * 

e*n correction of zp 

zpp-zp 

tf (sbs(zp) .gt.0.1) go to 63 
if Up) 60,61.61 

£0 zp*-sqrt Up rt *2+epsl2) 

go to 63 

61 zp*sqrt (zp rtft 2+epsl2) 

63 continue 

c * 

xlp*»xp+dx 




of>* ’ 

O* ' 


, '£» 








VI 


.215 


91 

93 


9 4 
96 


226 


xlm*xp-dx 
y'lp*yp+dy 
ylm®yp-dy 
xlp1“xlpftxlp 
xln2 a jlrn*xln 
ylp2-ylpftylp 
ylm2 a y Lipftyjm 
2l2“2p‘*Z{> * 

rl“sqrt(jhp2+ylp2+zi2) 
r2**qrt (xlm2+v\lp2+zl2) 
r3**sqrt (xlp2+ylra2+zl2) 
r4**sqrt (xlm 2 +ylm 2 + 2 l 2 ) 
arg l*»r 2ftzp*x tp*y 1 p-r 1 *zp*x lm>iy 1 p 
arg2«r<4*zpftxlpftylm-r3*zp*xlR!Aylm 
arg3»r l*r2*zl2+xlpi'xlnAylp2 
arg4«r3rtr4 ft 2l2+xlpAxlm*y lra2 
arg7“arg3ftargi» 
slgf«i. 

•. If (arg7. It.O.) slgf»-l. 
arg5®3lgfa(argl*arg4-arg2*arg3) 
arg&«8i,gfa (arg7+orgl*arg2) 
aj l*atan (arg5/arg6) 

If (argl.gt.O. .and.arg3.1 t.O.) aj l»aj 1+pl 
If (argl.lt. 0.. and. arg3.lt. 0.) ajl"ajl-pi 
I f (arg2.gt.O. . and. arg4. It.O.) aj l«*aj 1-pi 
If (arg2. 1 t.O. .and.arg4. 1 t.O.) ajl«*ajl+pi 
I f (argS.gt.O. .and.arg6. 1 t.O.) aj I»aj l+pi 
If (arg5.lt. 0. .and. arg6. It.O.) aj l«*aj 1-pi 
aj 2*»a I og ( (x lp+r 1) * (x ifit+r 4) / ( (x ln+r 2) * (x lp+r 3) ) ) 
kover-0 

If (abs(aj2) .lt. 5 . 03 ) go to 215 

aJ2“sign(5.03,aj2) 

kovar-kover+l 

continue 

tp l -0.00001 

i f (abs (y lp+r 1) . 1 1 . tpl) go to 91 
if (abs(ylm+r3) .lt.tpl) go to 91 
ajpl»xlp*aiog ((ylp+rl)/(ylra+r3) ) 
go to 93 
ajpl»0. 

if (abs(ylp+r2) .lt.tpl) go to 94 
If (abs (ylm+r4) .lt.tpl) go to 94 
aj p2«x lm*a 1 og ( (y 1 p+r 2) / (y ira+r 4) ) 
go to 96 
0 j p2«0 . 
aj3«ajpl-ajp2 

If (kovcr.no. 0) write (6,226) kover, 

Sx, y, z. x 1, y 1, z 1, x2, y 2, 22. x3. y3, z3, x4, y4, z4. xp, yp, zp. r l.r 2 , r 3, r 4, dx, 

6dy,aj2,aj3 

format (//lx,"-'' control point is too close to the edge of the vortex 

Gsheet * kover-", i 2/5x,9f 10. 6/5*. 9f l0.6/5x,6no.6, 2el5.6) 
argl-zp*(xip+y lp+r 1-xlm-y lp-r 2 ) 
arg2-zp>< (xlp+y lm+r3"xlm-y lm-r4) 
arg 3 “ 2 l 2 +(xlp+ylp+r 1) * (xlm+y lp+r 2) 
argi^-z 12+ (x Ip+y lm+r3) * (x lm+y lm+r4) 
arg7 B arg3 ,-; arg4 
sigf-1. 

if (arg7. 1 t.O.) sigf— 1. 
argp-sigf* {arg I>>arg4~arg2*arg3) 
arg6 B s igf * (arg7+arg l*arg2) 




i 

I 





10 , 


c * 
c * 
c * 
c * 
c * 
2 


c * 


337 


33*. 

7 

c * 
c * 
c * 


'* 80 

bj3°atan(arg5/arg6) 

If (argl .gt.O. .and.arg3. I t.O.) bj3"bj3+pi 
If (argl.lt. 0.. end. arg3.lt. 0.) bj3* , bj3*pl 
If (arg2.gt.Q. .and.arg4. 1 t.O.) bj3“bj3~pi 
If (arg2.lt. 0.. and. argli.lt. 0.) Bj3 H bj3+pi 
If (arg5- gt.O. .and. argfi. I t.O.) bj3°bj3+pi 
If (arg5* 1 t.O. .and.argfc. I t.O.) bj3»bj3-pi 
aJ3“aj3+2. ft 2P*bj3 

wyp»-fpift (gc*oj l+gvft(ypftaj l+zpp*aj2)) ^ ( 2 12-cps 12) 12 

•wzp»?.fpl ft (gc ft aj2+gv* (yp*aj2+aj3) ) 

• velocities ’jn tha general system of coordinates 
' ux«wyp’ft>tx2+w2p*xx3 
uy-4yp*yy2+wspftyy3 
uz"wyp*zz2+wzp*zz3 
return 


concetratcd straight vortex element 
continue 

dsnx*ds2*dyy-dsy*d22 
dcny*d3x*dzz"dsz*dxx 
dtraz=dsy*dxx-dsx*dyy 
dcm2°dsnx*dC3x+dsiny*d5Eiy+d3raz*dsmz 
If (dsra2. lo. . lc-20) go to 113 
fvds«fva/dsa2 
I0»l. 

r03«r02*sqrt (r02) 

If (rG2.gt.rnax2) go to 7 
zero Induced velocity on vortex Itself 
• If (r02.lt.. lc-20) go to l 13 

(dxx*d3x+dyy*dsy+d22*dsz) /r02 
If (abs(a-l.) .lo..00i) goto 113 
alpha2®da2/r02 
alpaa"alpha2-a*a 
If (alpha2+(2.*a) . I t.-l.or.alpha2- (2.*a) ,lt.-l.) goto 113 
If (abs(alpaa) . le..le-I5) go to 334 
sqla^sqrt (l.+2.*a+alpha2) 
sq2c a sqrt (l.-2.fta+alpha2) 

If ((sqla. lo. . lc-20) .or. (sq2a. le. . ie-20) ) write (6*337) sqla, 3q2a,x, 
Sy,z.xl,yl,2l.x3.y3.z3 

fornat (///lx,"**overf low*ft",3x,"sqla,sq2a",2ol5.5/3x,"x,y,2‘',3el5 
«.5/3x."xl.yi.zl,x3.y3.23",6oi5.5) 

If (sqla. le..lo -15) 3qla«3qla-«-0.0001 
If (sq2a. lo. . lc-15) sq28 e >3q2&+o.0001 
IO s »((alpha2+a) /sqla+(alpha2-a)/sq2a)/ (2.*alpaa) 
go to 7 

I0«j./((l.-a*a)**2) 

fact°-IOftfpiftgr3*2./r03 

Induced velocity multiplied by 1 . (1.+ (epa2/rho) ft*2) 

■I./(l.+fva/dstn2) 

zero Induced velocity on vortex line Itself 
if (dsm2.lt.. le-20) go to 113 
fvds®fva/dsm2 
fact“fact/ (l.+fvds) 

UX“f 3Ct>'<dSRlX 
uy«f act^dsmy 
uz*fact>'<dsm 2 
return 




113 


89 


c * 
c * 
c * 
c * 
77 
88 

c ft 

ft 


i 


1000 


1001 

1002 


ux«*0. 

uy»0. 

U2**0. 

return 


semi inf inite cyl inder 
continue 

continue 


otar»3qrt (x*ft2+yft*2) 

Zp**2pf-X 
url-0. 

do 1000 k»l,10 
psr».1571+(k-l) *0.311*2 
ul«rho*rho+fitarftetar+2p*2p-2ftrhoftotarftco3 (psr) 
u2®rho*0. 3ll*2* (rho-etar*cos (psr) ) / (etar*otar+rho*rho-2ftrho*etar* 
* &cos(par))*(t.-2p/(sqrt(ul))) 

U2l« (Qm/l2.56S4) *u2+u2 1 
continue 

U2»2ftl.ftuzl/del2 

sqr» (l»*etar*rho) / (2p**2+ (rho+etar) ft*2) 

8qrr°cqrt (sqr) 

sq«l.-sqr 

sqrp-sqrt (sq) 

sqrpp=3qrp**2 

If (sqr.eq.l.) go to 1001 

capk«*loa (I»./sqrp)-i-0.25* (1og(l*./sqrp) -1.) ftsqrpp 

C-1.+.5* (log (Wsqrp) -.5) ftsqrpp 

go to 1002 

capk»10.0 

e-1.0 

uh»- (2 ./sqrr) ftsqrt (rho/etcr) ft (capk* (2.-sqr) -2.*e) 
ps-atan2 (y,x) 

ux» (uhft 1 . ftcos (ps) ) ftgm/ (12. 566l* :ll de 1 2) 

uy- (uhft l . fts i n (pa) ) ftgm/ ( 12 . 5661 *ftde 1 2) 

return 

end 


of Fcoa QUAirrY 



subroutine wnofc2 90 

C * * ft 

C * * PROGRAM FWC * 

c * * * 

C ft * SUBROUTINE WN0FC2 * 

c ft ft ft 

C ft ftftftftftftftAftrtsVrtftft&ftftftftftftft 

c .* 

c >;;• 2/13/^s 

c ft • ' 

C ft EVALUATION OF THE IHOUCEO VELOCITIES AT THE NODES 
c * BY INTERPOLATION OF THE VELOCITIES INOUCEO AT THE CENTERS 

C * AND INTERPOLATION MEW/OLD DISTRIBUTION OF INDUCED VELOCITIES 

C ft 

c ft ENTRY INIWN: EVALUATION OF THE INTERPOLATION COEFFICIENTS 

c * 

ccmmon/lzmb/wls 

cownon/rol I6/ql,q2,q3.qL t lw 

cammon/roll l/harsjt5(26) ,hcmat6(26) 

common/al i/cgl,cg2,cg3,cgl», ir 

cera:on/roU2/wzlvr (8,50) .wjtivr (8,50) .wylvr (3.50) .wzr.vr (28. IS) , 
6wxnvr (28. 18) .wynvr (28. 13) 

ccmmon/parra/iwr . ird, i trace, lpsca.nmes, i trcw, i trcl 1, i trcl2,5 trcg, 

SI trcnt, i trcf , i plot, icgon, iplotv, I trta, ivcrgr , i test, issmeb, i samet, 
61 Iras, I ini, I in2, ni ter, i trans, ipn.j punch, Iplotw, Isw, icont, ip lotg, 
6lpr, Iplott 

common/ gccra/ r.b 1 ds l , nb 1 ds , s I gnia, f mu, otan (25) .knnvr.etai (6) ,knivr, 

61 twi st, thetad (25) , theta (25) ,thctac (21») > thataO, thotQd, alphas, 
6cd0,cdk,dpsind,dpsln,dpsl id.dpsi I ,coeff .ccoff l,c,s, blades, 
Snnvr.nnva.nnvrl.nnval’.nncr.nnca.etanv (25) .etanc(26) , 
6nivr,niva,nivr 1, nival, nicr,nica,etaiv(6) ,ctaie(7) , 
fintva.ntval,ntca,fpsl,fps2.nrl l,nr 1 2. i rl l, ir!2, i model 
common /rcsui/gammc (26) ,twx (2I*,2I*) , twy (21., 21*) . twz (21*,21») . 

6wxc (21*) ,wyc (21*) .wzc (21*) ,wxnc(26. 13) ,wync (26, 19) ,wznc(26. 19) , 

6wxie (S,5l) ,wyic(9,51) ,wzic(9,51) ,wrnv(28.18) ,wtnv(28.18) . 

6w2nv (23.18) ,wriv(8,50) ,wtiv(3,50) ,wziv(3,50) , 

Sxnc (26. 19) *ync(26.!9) ,znc(26.19) ,xic (9,51) ,y Ic (9.5D ,zic (9,51) , 
6xnv (28. 18) ,ynv(23.18) ,2nv(28.l8) ,xlv(8.50) ,ylv(8.50) ,ziv(8,50) , 
6la,ib 

common /cwnof /aan (25) ,aai (25) ,icb{25) ,acb(25) ,jcb(26) ,bcb(26) , 
6lcc(7) ,acc(7) ,jcc(3) ,bcc(3) ,frn,ftn,fzn 
common /wndata/wxnvt (28. 18) .wynvt (28. 18) ,wznvt (28, 18) 

6,wxlvt(8.50) ,wylvt(8.50) ,wzivt(8.50) 
ccman /savtw/twyt (2L.21*) ,twzt(2l*«21*) 

cojKr.on/ti pco/wxntl .wyntl ,wznt I,wxnt2.wynt2,wznt2,wxntp,wyntp,wzntp 
dimension wrnvt (28.1*1*) .wtnvt (28,1*1*) ,Hrivt(8.50) ,wtivt(3,50) 
dimension rs (28.50) 

equivalence (nnvr.nncrl) , (nivr.nicrl) . (nnva.nncal) 
equivalence (niva.nical) , (nnvr l,nncr2) , (ni vr 1 , ni cr2) 
equivalence (wxnvt(l.l) , wrnvt (1, 1)) , (wy-ivt(i.l) ,wtnvt(l, 1)) 
equivalence (wxivt(l.l) .wrivt(l.l)) , (wyivt(l.l) ,wtivt(l,l)5 
data laan, laai , 1 icb, lacb, 1 jcb/" AAN"," AA I " ICB"," ACB"," JCB"/ 
data Ibcb, I icc, lacc, I jcc, Ibcc/" BC3"," ICC"," ACC"," JCC"," DCC"/ 
data 1 xnv, I ynv , I znv , 1 r nv, l tnv/"WXNV" , "WYNV" , "WZNV" , "WRNV" , "WTNV" / 
data Ixiv.lyiv, Iziv, 1 r i v, 1 tiv/"WXI V" ,"WYI V","WZ! V" ,"WRI V'V'WTIV"/ 
data l twy, 1 twz/"TWY ","TVZ "/ 
data dps2,fpi/100., 12.566371/ 

c * 
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goto 70 
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c 

c 

c 

c 

c 

c 

c 

c 

c 


50 

c 


51 

c 

c 

c 

c 


53 

54 

52 

c 

c 

c 

c 


56 

57 

55 

c 

c 

c 

c 

c 


59 

60 

53 


* 

ftftftftftftftftftftftftftftftftftftftrtftftftftftftftftftftftftftfcftftftftftftft 

ft 

entry Sniwn2 

* 

if (isamet.cq. 1) return 

.oittS SflAJ.PAG^i& 

* I HTERPOLATI OH COEFFICIENTS G? F8Q3 QUAUTV 

* / • ‘ 

ft COEFFICIENTS AAN : NODE I VS. CENTERS I AND 1+1 

* FOR THE NEAR WAKg. RADIAL DIRECTION 

do 50 i**l,nnvr 

aan(i)« (etanc (i+1) -etanv (i) ) / (etanc (1+1) -etanc (i)) 
continue 
ft 

do 51 i»l,nivr 

aai (i) « (etaic (i+1) -ctaiv ( i > ) / (etaic (i+1) -etaic (i)) 
continue 

* 

ft COEFFICIENTS ACB ANO THE CORRESPONDING ADRESSES ICB : 

* NODE I OF THE NEAR WAKE VS. THE CENTERS ICB(I) AND ICD(I)+1, 
ft OF THE I NT. WAKE (RADIAL DIRECTION) 

do 52 i a l,nnvr 

et»etanv (?) 

do 53 i l“2,nicr 

if (et. 1 1. etaic li 1)) goto 54 

continue 

icb (i) «i 1-1 

acb (?) "(etaic (i 1) -et) / (etaic {» 1) -etaic (i 1-1)) 
continue 
ft 

ft COEFFICIENTS 3CB ANO THE CORRESPONDING ADRESSES JCS : 

* NODE J OF THE TRANSITION WAKE VS. CENTERS JC3 (J) AND JCB(J)+1 
ft OF THE INTERMEDIATE WAKE, AZIMUTAL DIRECTION 

do 55 j-ntca,nnva 
et"(j-l) *dpsin 
do 56 j 1-2.4 

eti- (ntva-1) ftdpsin+(j 1-1) ftdpsi 1 —dps i i*.5 

if (et. 1 t.eti) goto 57 

continue 

jcb(j)"jl-l 

beb (j) ■ (et?-et) /dpsi i 

continue 

ft 

* 

* COEFFICIENTS BCC AND THE CORRESPONDING ADRESSES JCC : 

* NODES J OF THE I NT. WAKE VS. THE NODES JCC (J) AND JCC(J)+! 
ft OF THE NEAR WAKE (AZIMUTAL DIRECTION ) 

do 58 i«l .nivr 

et-etaiv(i) 

do 59 i l-2,nnvr 

i f (et. 1 1. etanv (i 1) ) go to 60 

continue 

icc (i) a i 1-1 

acc(i) = (etanv (i 1) -et) / (etanv (i 1) -etanv (i 1-1)) 
continue 
ft 

do 61 j“ l, 3 

6t“(ntva-l)*dpsir.+ (j-l) ftdpsi i 


92 


62 

63 

61 

- c 1 
’ *■ ' 

200 


0 s v 


c 

e 

c 

70 

c 

c 

c 

e 

c 


do 62 jl*ntca,nnva 
eti*(j 1-1) *dpsin 

if (j.eq.l. and. jl.eq.ntca. or. at. It. ati) goto 63 

continue 

jcc(j)**j 1-1 

bcc (j) ** (ot i -et) /dps i n 

continue / 

i . * 1 .* 1 * / 

If, (I traco.eq.O) return 
wr i te (iwr ,200) 
format (" *** IHIWNQFC ***") 

> ;wM te (iwr , 101) laan, (aan(i) , i**l,nnvr) 
writo(iwr,10l) laai, (aai (i) ,i*»i,nivr) 
wrl to (iwr , 102) 1 icb, (icb(i) , i»l,nnvr) 
write(iwr,101) lacb, (acb(l) ,i®l r nnvr) 
write (iwr , 102) i jcb, (jcb (i) , i>*ntca,nnva) 
write (iwr, 101) Ibcb, (bcb (i) , i«ntca,nnva) 
write (iwr, 102) 1 icc, (i cc (i) , i**l,nivr) 
wrl to (iwr, 101) lacc, (acc (i) , i«*l,nivr) 
wrl te( iwr, 102) 1 jcc, U« (i) , i**l,3) 
write (iwr, 101) lbcc, (bcc(i) , 1*1,3) 
return 


****************** *ftftft*ftftft****ft**ftft**ftftft*ft*ftftftft**ftft*ft ft* 

* 

continue 






ft 

ft 

ft 

ft 

ft 


179 


INTERPOLATION OF THE VELOCITIES 

.A) INTERPOLATION OF THE VELOCITIES INDUCED AT THE CENTERS 
OF THE NEAR WAKE 1 TO NNCA TO THE NODES 1 TO NNVA 
do 1 i*l,nnvrl 
do 1 j«l,nnva 
a*aan(i) 
ai*l.-a 

wxnvt (i , j) n ( (wxnc (i , j) +wxnc (i , j+l) ) *a+ (wxnc (i+l, j) + 

Cwxnc (i+1, j+l) ) *ai) *.5 

wynvt(i ,j)**((wync(i ,j)+wync(i ,J+1)) *a+(wync(t+l,j) + 
twync (i+1, j+l) ) *al) *.5 

W2nvt(i , j) = ((w2nc (i ,j)+w2nc(i , j+l) ) *a+ (wzne (i+1, j) + 

Awznc (i+l,j+l))*al) *«5 
continue 
a*aan(nnvrl) 
al“l.-a 

wxnvt (nnvr 1 , 1) *wxntp+ (wxnc (nnvr 1 , 2) -wxr>.t 1) *a+ (wxnc (nner 1 , 2) -wxnt2) *al 
wynvt (nnvrl, l)*wyntp+(wync (nnvr 1,2) -wyntl) *a+(wync (nner 1,2) -wynt2) *al 
W2nvt (nnvr 1,1) < c wzntp+ (w2nc (nnvr 1 , 2) -W2nt 1) *a+ (wane (nner 1 , 2) -W2nt2) *a 1 
if (Itrcw.ne.2) goto 10 
write(iwr,179) 

format ("1*** WNOFC *** VELOCITIES AT THE NODES OF THE NEAR WAKE") 
do 11 j»i,nnva 
wr i te (iwr , 103) 


11 

c 

c 

c 

c 


wr i te (iwr , 101) Ixnv, (wxnvt (i , j) , 
wr i to (iwr , 101) lynv, (wynvt (i , j) , 
wr i te (iwr, 101) lznv, (wznvt (i ,j) , 
continue 


•l.nnvr) 
*1 .nnvr) 
■l.nnvr) 


B) INTERPOLATION OF THE VELOCITIES AT THE CENTERS 1-2-3-L 
(AZIKUTAL POSITIONS OF THE I NT . WAKE) AT THE NODES OF THE NEAR WA 
FROM NTCA+1 TO NNCA 





93 


10 


2 

178 



c 

c 


c 


2C0C 

1001 

c 

c 

c 


c 


l t (i trans.eq.O) goto 12 
do 2 i"l,nnvr 
i l«icb (1} 

8«aob (i) 
alM.-a 

do' 2 j»ntca,nnva 
jl“jcb(j) UKiUt.v£^VmaiL lb 

b-bcb(j) OF POOS? QUALITY 

bl«l .-b 

1 f (j .eq.ntca) 0*0.5 
wxnvt {! * j) “wxnvt (i , j)*( ( 
twxic (i 1 , j 1) *&-*-wxi c (1 1 , j 1+1) *bl) rt c+( 

Cwxic (i 1+1 , j 1) ^b+wxic (i l+l , j 1+1) *bl) *a I) 
wynvt (i ,j) -wynvt (i . j) + ( ( 
fcwyic (i 1 , j 1) *b+wy i c (i 1 , j l+l) *bl) *a+( 
twyic (i l+l , j 1) ftb*xy ic ( i l+i , j l+l) *bl) *a l) 
wznvt (i , j) -wznvt (i ,j) + (( 
twzic (1 1, j 1) Ab-rwzic (i 1 , j l+l) *bl) *«:+( 
twzic (i l+l, j 1) *b+wzic (1 1+ 1 , j 1+1) <*bl) *al) 
continue 

if (i trcw.ne.2) goto 12 
write (iwr , 178) 

fornat ("INODES HTCA.NNVA, AFTER ADDITION CF THE VELOCITIES", 

6" INOUCEO OH CENTERS l TO i* OF THE INT.WS") 
do 13 j»ntca,nnva 
wr i ta(lwr , 103) 

wr 1 to (iwr , 101) 1 xnv, (wxr.vt (i , j) , i*l ,nnvr) 
wr i te (iwr , 10!) lynv, (wynvt (i , j) , i“l ,nnvr) 
wr i te (iwr . 101) Iznv. (wznvt (1 , j) , i«l , nr.vr) 
continue 

ideb*5 

TOTAL INDUCED VELOCITY AT INTERRED I ATE WAKE 

if (i trans.eq.O) ideb“l 
nrl21-nr 12-1 
nr 122»nrl2“2 
Irl21-ir 12-1 
irl22«ir!2-2 
do 1001 i*ir 1 1 , ir 12 
do 2000 j-2,niva 

ra (i , j) -sqrt (xiv (i , j) **2+y iv (i , j) **2) 
if (1 .cq. i r 1 1) gni*ql 
if (i .eq.ir 122) gmi“o2 
if (i . cq.irl21) gmi-q3 
i f (i .eq. i r 1 2) gmi “qL 

wzivtti ,j)-wzivr (i , j)+(gmi*wi s) / (i» *3 • U» 159 Ar s (i , j) ) 

wxivt (i , j) -Wxiv (i . j) 

wy ivt (i , j) -wy ivr (i , j) 

conti nue 

ccnti nue 

TOTAL INDUCED VELOCITY AT HEAR WAKE 

do 5016 i -nr 1 1 , nr 1 2 
ao 5016 j“2, n:ivn 

rt ( i , j) -nqr t (xr.v ( i . j) >';^2+ynv ( i , j) **2) 
i f (i .cq.nr 1 1) gni 1 
i f (i .eq.nr 1 22) gn i »q2 



?4 


5016 
c * 


111 

222 


If (1 .eq.nr 121) gni-qj 
If (1 .aq.nrU) gmi-ql. 

wznvt (i , j) -w 2 nvr {! . j) { 5 m! Awl s) / (I,*'). Ul53*«“3 (i . j) ) 
wxnvt (i , j) -wxnvr (i , j} 
wynvt (i ,j) --ynvr (1 , j) 
continue 






O* 


da 111 I«2,nncri 
I f (otanc (!) .gt.cgl) go Co 222 
ibol«i 
pol 1-cgl-etanc ( i 00 1) 
lboll-ibol-1 

pol2-etanc (ibol+l) -atone (I bo 1) 

wznvt (nr! 1 . 1) -wzc (I bo l l) + (wxc (I bo l HI) -wzc (I bo 11)} a (pot l/pol2) 
wynvt (nr 11 « 1) -wyc (I bo 1 1) + (wyc (ibol HI) -wye (1 boll) ) a (pol l /pol 2) 
wxnvt (nr 1 1 , l) -wxc (ibol l) + (wxc (Ibol HI) -wxc (Ibol 1 }) a (pol t/po 1 2) 


e 

c 

c 

e 333 

e LU 

c 

c 

e 

e 

c 


do 333 l» 2 ,nncrl 
if (atonc(i) .gt.cg 2 ) go to Ukk 
lbo 2 «l 

pol 3 "(cg 2 *-«rtanc (ibo 2 ) ) / (otanc (lbo 2 -M) -otanc (ibo 2 ) ) 
lbo 2 l«ibo 2 -l 

wznvt (nr 122 , 1 ) -wzc (lbo 21 ) + (wxc (ibo 2 Ht) -wzc (lbo 2 l) ) Apot 3 
wynvt (nr 1 22 . 1 ) -wyc (ibo 2 1 ) -Mwyc (i bo 2 Hl) -wye (ibo 2 i) ) Apol 3 
wxnvt (nr 1 22. 1) -wxc (i bo2 1) ♦ (wxc (i bo2 HI) -wxc (ibo2 1 } ) Apol 3 

do 555 i« 2 .nncrt 
if (etanc(i) .gt.c;3) go to 666 
555 lbo 3 «i 

666 poU»(cg 3 "otanc(ibc 3 ))/ (otanc (5bo3+l) -etanc (ibo 3 )) 

ibo3l"ibo3-l 

wznvt (nr 1 2 1 , l) -wzc (i fco 31 ) + (wzc ( i bo 3 t-*-X J -wzc (ibo 3 1 ) ) adoII* 
wynvt (nr 1 2 1, 1) -wyc (i fco31) *(wyc (ibo3Hl) -wyc (i bo3 1) ) Apol It 
wxnvt (nr 121, 1) «*wxc (ibo3l)*(wxc (ibo3Hl) -wxc (ibo30 ) Apol l, 

do 777 i « 2 .nncrt 
If (etanc (i) .gt.cgl,) go to 863 
777 lboi»-i 

883 pol5® (cgi,-etanc (iboi,) ) / (etanc (IboL+ 1 ) -etanc (ibol,) ) 

JboLl-lbo!»-l 

w 2 nvt (nr 12 . 1 ) -w*c (ibol,l) + (wzc (ibol, 1 + 1 ) -wzc 0 bol»l) ) Apol 5 
wynvt (nr 12 . i) -wyc (ibol, l) + (wyc (<bol*l-*-l) -wyc (ibol, l) ) Ap 0 l 5 
wxnvt (nr 1 2. l) -wxc (ibol*l) + (wxc (ibol, HI) -wxc (ibol, 1) ) Apol 5 

ft 

* TRANSFORMATION J RECTANGULAR TO POLAR COORDINATES 
A OF THE VELOCITIES INDUCED AT THE NODES 


c 

c 

c 

c 

c 

16 

L60 


contlnuo 

i f (i trcw.cq . 2 ) wr i to { iwr , L60) 

format (" 1 *** WNOFC Ana. VELOCITIES 

do 1,5 j-l.nnva 

do 80 i-nr 1 I , nr 1 2 

wxx-wxnvt (i , j) 

wyy-wynvt (i . j) 

xx-xnv (i , j) 

yy»ynv (i . j) 

r-sqrt (xxAx**vyrtyy) 

wrnvt (I , j) - (xx-'wxxf yyAwyy) /r 


AT THE NODES, POLAR C00R."5 



95 


wtnvt (i . j) ■ (“w\x*yyt-nyy ft *x) /r 
SO conti nuo 

I f (i trcw.ne. 2 ) goto t»$ 
wr I to (iwr , 1 03 ) 

wr I te (iwr , 10 1) 1 rnv, (wrnvt ( ! , j) , i« 2 ,nnvr) 
wr i tc (iwr , 10 1 ) 1 tnv, (wtnvt ( 1 . j) , i« l .nnvr) 

1,5 continue URiGlW/U. PAGE (3 

do 1,6 j- 2 .niva OF POOH QUALflY 

00 8 1 i " i r > l , i r 1 2 

• • 

. m WXX-WxiVt (i , j) 
wyy-wylvt (i , j) 
xx«xiv (i , j) 
yy-yiv(i.j) 
r-tqrt ( *xftxx+yyftyy) 
wr ivl (i , j) - (xx*wxx*yyftwyy) /r 
Wt ivt (i , j) * (-wxxftyy+wyyftxx ) /r 
cont i nua 

if (i trcw.no. 2 ) goto 1,6 
wr I to (iwr , 103 ) 

wr I te (iwr , 101 ) Ir Iv, (wr Ivt (i , J) , i«l r 1 1 , ir 12 ) 
wr i te (iwr , 10 l) l t i v, (wt 1 vt (i , j) , i ■! r t 1 , i r 1 2 ) 
continue 

wr ivt (i r 1 1 , l) -wrnvt (nr 1 1 ,nnva) 
wt i vt (i r 1 1 , 1 ) -wtnvt (nr 1 l ,nnva) 
wr ivt (i r 1 2 , l) -wrnvt (nr 1 2 ,nnva) 
wt ivt (i r 1 2 , l) -wtnvt (nr 1 2 ,nnva) 
wzlvt (ir I 2 , 1 ) -wznvt (nr 1 .'.nnv;) 
wzivt (ir 1 1 . 1 ) -wznvt (nr 1 1 .nnva) 
wr ivt (i r I 2 1 , !) -wrnvt (nr 1 21 ,nnva) 
wt i vt (i r l 2 1 , 1 ) -wtnvt (nr I 21 , nnva) 
wzivt (1 r 12 1 . 1 ) -wznvt (nr 12 1 ,nnva) 
wr Ivt (Ir 122, w -wrr.vt (nr 122. nnva) 
wtlvt(irl 22 .!)-wtnvt (nr 1 22 , nnva) 
wzi vt (ir I 22 , l) -wznvt (nr 122 . nnva) 

c A 

C * INTERPOLATION OF THE 'OLD' AND THE 'NEW' DISTRIBUTIONS 
c * OF VELOCITIES 
c ft 

1 f (i trcw.no. 0 ) wri to (iwr, 1,62) 

1,62 format ("iftftft WNOFC ftftft : VELOCITIES AT THE NODES ") 

c ft 

C ft WEIGHTING FACTORS FOR THE 'HEW' AND THE 'OLD' VELOC'TIES 
f 1 . t, 

f tn-.7 
ftn-.J 

if (ni tor. 30 . )) goto US 
frn-.j 
f tr-.S 
fzn-,3 

i f (ni tor . eq . 2 ) goto L 9 

frn-o 

f tn-.? 
fzn *. 3 

US contini -3 

if ( i t t cw . r e . 0 ) wr i t e ( i wr , nr. j) frn,f;n,fzn 

f orr-.a • ; 1 WtlGHT.SG UCTOHS: FRN.FTN.F2N ^ f 6 . 2 ) 
f rv*- 1 . - 1 1 n 



951 



ftv*l.*f tn 
f zv«l ,-fzn 


96 


c 


648 

1.3 


k 7 


649 

44 


48 

c 

c 

e 

c 

e 


A 

do 1*7 j“ 1* nova 
do 43 « “or I l,nrl2 

wrnv (i . J) "wrnvt d . j) »'f rrvH*rnv (i , J) <*f rv 
, wtnvd , j) -wtnvt d , j) *f tn+utnv (i , j) *f tv 

w 2 nv(l . j) -wznvt (I ,J) *f 2 n+wznv (i , j) *fzv 
go to 43 

wrnv (I , j) «*wrnvt (J . j) 
wtnv(l , j)-wtnvt (I ,J) 
wznvd , j) “wznvt (i ,j) 
conti nua 

If (itrcw.eq.O) goto 47 
wr I to (i wr , 103) 

wrlte(iwr, 101) Irnv, (wrnv(i, j) ,l*»l,nnvr) 
wrl te (Iwr, 101) 1 tnv, (wtnv(l . J) , l»I ,nnvr) 
write (Iwr, 101) Iznv, (wznv (I ,j) , i»l,nnvr) 
continue 
do 48 j“l,nlva 
do 44 IHrl l,lr!2 

wr Iv (I ,J)*wr Ivt (1 ,j) ftf rrt+wr I v (1, j) *frv 
wtlv(i , j) "wr ivt (l ,J) *f tn+wtiv (I , j) *f tv 
wzlv(l , j) -wzlvtd , j) efzn+vziv (i , j) ftfzv 
go to 44 

wr iv (I , j) «wr Ivt (l , j) 
wtlvd ,j) «wtlvt(i ,J) 

W 2 iv(I ,j)*WZlVt(l.J) 
continue 

If. (itrcw.eq.G) goto 48 
writo(Iwr, 103 ) * 

wrl to (Iwr, 101) Irlv, (wr iv(! ,j) , t*»Irl 1, Ir 12) 
wrl to (Iwr , 10 1) 1 tlv, (wtl v (I , j) , i«*I r 1 1, 1 r 1 2) 
wrl to (iwr , 101) Iziv, (wzlv{l.j) .!«Irll,lrl2) 
continue 

* 

AAAAAAAAAAAAAAAAANAAAAiRAAAAAAAAAAAAAAAAftAAAAAAAAAAftA 

A 

* INTERPOLATION OF THE INFLUENCE COEFFICIENTS 

A 



if (niter. eg. 1) return 

fzn*».5 

fzvM .-tzn 

I f (I trcw.eq.0) goto 90 
wr I te (iwr , 500 ) 

500 formatC *** WNOFC: INTERPOLATION OF INF. COEFF’",/, 

S" INF. COEFF FROM L00P2S") 
do 91 j-l,nncr2 

write (iwr, 101) Itwy, (twy(l.j) , I»l,nncr2) 
wr I te (iwr, 101) I twz, (twz (i , j) , !■! ,nncr2) 
wr i te (iwr . 103) 

91 continue 

wrl te (iwr, 501 ) 

501 formatC 0L0 INF. COEFF. ") 

do 92 j"l,nncr2 

wr i te (iwr , 10 1 ) I twy , (twy t (i , j) , I- 1 ,nncr2) 
wr i te (iwr , 101) I twz, (twzt (i • j) • i*l , nncr2) 
wr I te (iwr, 103 ) 

92 continue 

90 continue 
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do 82 i“1,nncr2 
do 82 j«l ,nncr2 

twy (i , j) »f zn^twy (i , j)+fzv*twyt (I ,j) 
tWZ (i ,j) -fzn^twz (i , j)+f2V*tW2t (i ,j) 

continue 

I f (I trcw.eq .0) return 
wr I to (iwr , 502 ) f zn, f zv 

. format (" I NT . INF. COEFF.: F2N-",e20.5," F2V-", 
do 95 J“l,nncr2 

wr i to (iwr , 101) I twy, (twy (i ,j) , 1-1 ,nncr2) 
wr 1 to (iwr , 101) ltwz, (twz(I ,j) ,i-l,nncr2) 
write (iwr, 103) 
conti nuo 
return 

formst(lx,ai,.lx,5M3.9./.6*.5fl3.9./.6x,5M3.3) 

format (lx,al«, Ix,9i 1 3./ . 6x . 9 i 13) 
format (lx) 


ORIGINAL PAGti lb 
OF POOR Qt)AL(T\ 

o20.5) 



-■1 


subroutine wnofc3 

c * 

C * fteartAytnrtrtArtrtftftrtrtrtiirtitfte 

c * ft * 

c ft ft PROGRAM FWC ft 

c ft ft ft 

c ft ft SUBROUTINE WN0FC3 ft 

c * ft ft 

' c ft eftOftfteftenftftAftftaftrtftAftftft 

e * 

c ft 2/13/78 

C ft 

c * EVALUATION OF THE INDUCED VELOCITIES AT THE NODES 
C * BY INTERPOLATION OF THE VELOCITIES INDUCED AT THE CENTERS 

C * AND INTERPOLATION NEW/OLO DISTRIBUTION OF INDUCED VELOCITIES 

c * 

e ft ENTRY INIVN: EVALUATION OF THE INTERPOLATION COEFFICIENTS 

c * 

ccmraon/lamb/wis 
common/al i/cgl,cg2,eg3,cgL, Ir 
ccmmon/rol l6/ql,q2,q3.qli. iw 
eemmon/rol I l/hamrat5 (26) ,hammt6 (26) 

eemmon/rol 12/wzivr (3.50) .wxivr (8.50) .wylvr (8.50) .wznvr (28. 18) . 
twxnvr (28. 18) .wynvr (28. 18) 

correaon/parm/iwr , Ird, i trace, I psora, nines, i trcw, i trcl l, i trcl2, 1 trcg, 
ti trcnt, i trcf , iplot, iegen, iplotv, I trtg, ivergr , I test, i samab, i siir.et, 
tl ins, I iral. I Ira2.nl tar, I tron3, Ipn.j punch, ipiotw, isw, Icont, iplotg, 
tipr, Iplott 

ccnrnon/gectn/nblds 1 .nblds. s I gma.fnu.etan (25) .knnvr.otai (&) .Itnivr, 
Altwi*t.thotad(25) .thata(25) ,thctac(2i.) .thetaO.thetOd, alphas, 
tcdO.cdk.dpsind.dosin.dosi Id, dps! i ,coef f ,coef f l ,c, s ,bl ados, 
6nnvr,nnva,nnvrl,nnval,nncr,nnca,atanv(25) ,atanc(26) . 
tnlvr ,nlva, nivrl.nl val.nicr.ni ca, eta iv (b) .ctaic (7) , 
tntva.ntvai.ntca. fpsl. fps2,nr 1 l.nr 12, ir 1 1, Ir 12. i model 
common /resul /gamine (26) , twx (21,, 24) , twy (2L.2I.) , twx (2l»,2i.) , 

Cwxc (21*) ,wyc (2L) ,wzc (2L) ,wxnc (26. 19) ,*ync (26. 19) ,wznc (26. 19) , 
twx ic (9. 51) .wy ic (9.5U .wzic (9,51) ,wrnv (28. 18) .wtnv (23. 18) . 
lwxnv(28.18) .wriv(8.50) ,wtiv(8,5C) ,wziv(8,50) , 
txnc (26. 19) ,ync (26. 19) .znc(26.l9) ,x» c (9.51) ,yic(9,51) ,zlc(9.51) . 
txnv (28. 18) ,ynv (28, 18) ,znv(28.18) ,xlv(8,50) ,yiv(8.50) ,zlv(3,50) , 
tia, lb 

common /cwnof /aan (25) ,aai (25) . i cb (25) ,acb(25) ,Jcb(26) ,bcb(26) . 
ticc(7) ,acc(7) ,jcc(3) ,bcc(3) .frn.ftn.fzn 
common /wndata/wxnvt (28.18) .wynvt (28. 18) .wznvt (28. 10) 
t.wxl vt (8.50) .wy ivt (8.50) ,wzi vt (8.50) 
common /sovtw/twy t (21*. 21*) , twzt (21* , 2I«) 

convnon/tipco/wxntl.wyntl.wznt l ,wxnt2,wynt2,wznt2,wxntp,wyntp,wzntp 
dimension wrnvt (28.LL) .wtnvt (28,1*1*) ,wrivt( 8 . 50 ) .wtivt ( 8 . 50 ) 
dimension rs (23 . 50) 

equivalence (nnvp , nner l) , (ni vr ,nl cr 1) . (nnvs, nnca 1) 
equivalence (niva.nical) , (nnvr 1 ,nncr2) , (nivr l ,nicr2) 
equivalence (wxnvt(i.l) .wrnvt(l.l)) , (wynvt(l.l) .wtnvt(l.l)) 
equivalence (wxivt(l.l) ,wr ivt (1 , 1) ) , (wyivt(l.l) .wtivt(l.D) 
data laan, laai , I icb. Ucb. 1 jcb/" AAII 1 ', " AAI"," IC3"," ACB" JCS"/ 
data Ibcb, I Icc. lacc. I jee, ibcc/" 3CB"," ICC"," ACC"," JCC" ," 3CC"/ 
data lxnv, lynv, Iznv, Irnv. I tnv/"WXNV","V/YNV","WZMV","VRNV" ,"WTNV"/ 
data Ixi v, ly iv, Iziv, Ir iv, I tiv/"WXI V","WYIV","WZIV","WRIV","WTI V"/ 
data Itwy, I twz/"TWY ","TWZ "/ 
data dps2,fpi/100.. 12.566371/ 

c * 


9C 


.A 1 * 1 

J o^° 



e 

e 

c 


goto 70 


99 


c 

c 

c 

c 

c 

c 


50 
c 

51 

c 

c 

c 

c 


53 

54 

52 

c 

c 

c 

c 


56 

57 

55 

c 

c 

c 

c 

c 


59 

60 

58 


ft 

ftrtftftftftftftrtftrtitftftftftftftftrtftftftftrtftftftftftrtftnrtftftftftrtrtrt 


ft 

A 

A 

A 

A 

ft 

A 


A 


ft 

ft 

ft 

A 


A 

A 

A 

A 


entry iniwn3 
I f (isamet.eq. 1) return 
INTERPOLATION COEFFICIENTS 


ORIGINAL PASS {3 
OF POOR QUAUTV 


COEFFICIENTS AAN : NOOE I VS. CENTERS I AND 1+1 
FOR THE NEAR WAKE, RADIAL DIRECTION 
do 50 i"l,nnvr 

aan (I) - (etanc (l + l) -etanv (i) ) / (etanc (i+i) -otanc (I) ) 
continue 

do 51 1 - 1 .nlvr 

aai (I) - (etaic (i+1) -etaiv (i) ) / (etaic (i+l) - etalc (1) ) 
continue 

COEFFICIENTS ACB AND THE CORRESPONDING AORESSES ICB t 
NOOE I OF THE NEAR WAKE VS. THE CENTERS I CD (I) AND ICQ(l)+l, 
OF THE INT. WAKE (RADIAL DIRECTION) 
do 52 i-l,nnvr 
ct-etanv (i) 
do 53 I l-2,nicr 
If (et. I t.etaic (i 1) ) goto 54 
continue 
icb (i) -i l-l 

acb (i) - (eta ic (1 1) -et) / (eta I c (1 1) -etaic (I l- 1)) 
ccnr. inua 


COEFFICIENTS SC9 AMD THE CORRESPONDING AORESSES JC3 : 

NODE J OF THE TRANSITION WAKE VS. CENTERS JC3 (J) AND JC3(J) + t 
OF THE INTERMEDIATE WAKE, AZIMUTAL OIRECTICN 
do 55 j-ntca.nnva 
et* (j-1) *dpsin 
do 56 j 1-2.4 

etl-(ntva-l) *dpsln+(j 1-1) *dpsi i-dpsi I*. 5 

1 f (et. I t.eti) goto 57 

continue 

jcb (j) - j l-l 

beb (j) - (eti -ct) /dpsi i 

continue 


COEFFICIENTS BCC ANO THE CORRESPONDING AORESSES JCC : 

NOOES J OF THE INT. WAKE VS. THE NODES JCC (J) AND JCC(J) + 1 
OF THE NEAR WAKE (AZIMUTAL DIRECTION ) 
do 53 i-l ,nivr 
et-utaiv (i) 
do 59 i !"2 , nnvr 
i f (ct . 1 1 .etanv (i I) ) go to 60 
cont i nue 
ice (i) -i 1-1 

acc (i)- (etanv (i t) -et) / (etanv (i !)-etsnv(i 1-1)) 
continue 


ft 


do 61 j*l .3 

ct - (i-tva- 1) fteps i nr (j - 1) «dps i i 



100 


62 

63 

61 

c ' 


200 


c 

e 

c 

70 

c 

c 

e 

c 

e 


do 62 j l“ntca,nnva 
etl"(j l-l) <»dpsin 

if (j.eq.l. and. ji.eq. ntca.or.et.lt. eti) goto 63 
continue • t . , a, 

Jcc(j)-j l-l ** . 

bcc (j)- (etl-et)/dpsin 
continue 

i f (i trace. oq.C) return 
wrl te( Iwr, 200) 
formatC* *** INIWNOFC **0") 
writ* (lwr, 101) loan, (aan(i) , l-l,nnvr) 
write (Iwr, 101) laoi , (aai (i) , i«l,nlvr) 
wr I te (Iwr , 102) l icb, (icb (I) . i"l,nnvr) 
wri te (Iwr, 101) lacb, (acb (i) , !■! ,nnvr) 
wr I te (Iwr , 102) I jcb, (job (l) , i»ntca,nnva) 
wr I te (iwr, 101) tbcb, (bcb (I) , i«ntca,nnva) 
write (Iwr, 102) lice, (lcc(i) ,1-l.nivr) 
write (Iwr, 101) lacc, (see (I) , l«l,nivr) 
wrl to (Iwr, 102) ijee, (jee (1) , I *1,3) 
wrl to. (iwr, 101) I bcc, (bcc (i) • I “l, 3) 
return 

* 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAftAftAAftAAAAAA 






of 




continue 

INTERPOLATION OF THE VELOCITIES 


A) 

OF 


OF THE VELOCITIES INCUCE9 AT THE CENTERS 
TO NNCA TO THE NODES 1 TO MNVA 


179 


11 

c 

c 

c 

c 


INTERPOLATION 
THE NEAR WAKE I 
do 1 l-l, nnvr l 
do 1 j»l,nnva 
a-aan(i) 
al»I.-a 

wxnvt (i , j)“ ( (wxnc (i , j)+wxnc (i , j+1) ) *a+(wxnc (i+1, j) + 
fiwxnc (S + l , j+1) ) *al) *.5 

wynvt (I , j) « ( (wync (i . j) +wync (I , j+1) ) *a+(wync (l + l, j) + 

Swync (1 + 1 , j+l) ) >'<al) *.5 

wznvt (I , j)*» ( (wznc (i , j)+wznc (i ,j+l) ) na+(wznc (1+1, j) + 

&wznc(i + l,j+l))*al)*.5 
continue 
a«aan(nnvrl) 
al»l.-a 

wxnvt (nnvr 1 , 1) -wxntp+ (wxnc (nnvr 1 ,2) -wxnt 1) *a+ (wxnc (nner l , 2) -wxnt2) 1 
wynvt (nnvr 1,1) =wyntp+ (wync (nnvr 1,2) -wyntl) *a+ (wync (nner 1,2) -wynt2) >'al 
wznvt (nnvr 1 , 1) -wzntp+ (wznc (nnvr 1 , 2) -wznt 1) >'a+ (wznc (nner 1 , 2) -wznt2) *a 1 
I f (i trcw.no. 2) goto 10 
wri to(lwr, 179 ) 

format ("1*** WNOFC *** VELOCITIES AT THE NODES OF THE NEAR WAKE") 
do 11 j-l,nnva 
wr 1 te (iwr , 103 ) 

wr 1 te (iwr , 101) Ixnv, (wxnvt (i , j) , !■! .nnvr) 
wri te (iwr, 101) 1 ynv, (wynvt (i ,j) , i»l,nnvr) 
wr i te (iwr , 10 1) 1 2 nv, (wznvt (i , j) , i *•! .nnvr) 
conti nue 


A 

A 

A 

A 


B) INTERPOLATION OF THE VELOCITIES AT 
(AZIrtUTAL POSITIONS OF THE I NT. WAKE) 
FROM NTCA+1 TO NNCA 


THE CENTERS J-2-3-L 

AT THE MOOES OF THE NEAR WA 



101 


10 


2 

178 


13 

c * 
12 
c 
c 
c 


c 

c 


c 


2000 

1001 

c 

c 

c 


c 


If (i trans.eq.O) goto 12 
do 2 i«l,nnvr 
1 l*»lcb ( • ) 
a«acb (i) 
al-I .-a 

do 2 j*ntca,nnva 
j l"jcb (j) 
b*bcb (j' 
bl**l .-b 

I f (j .eq.ntca) b«0.5 
wxnvt (i , j) -wxnvt (i , j) + (( 

Swxic (1 1 ,j 1) *b+wx i c (i 1 , j l+l)ft'oI) *a+( 

Swx ic (1 1+1, j 1) ftb+wxic (i l+i,ji+l)ftbl)*al) 
wynvt (1 , j) «wynvt (1 , j) + { ( 

6wyic(i l,jl)*b+wyic(l 1 , j l+l) ftb 1) *a+ ( 
Swyic (i l+t,j l) ftb+wyic (i l+l,ji+l) *b 1) ftal) 
wznvt (I , j) ” vznvt (i , j) + ( ( 

Cwzic (i 1, j 1) ftb+wzic (i 1, j 1+1) ftbl) *a+( 
fiwz i c ( i 1+1, j l) ftb+wzic {i 1+1 , j 1+1) *b 1) *a 1) 
continue 

if (i trcw.ne.2) goto 12 
wri te (iwr , 178) 
format ("1N0QES 
6" INDUCED ON 
do 13 j*ntca,nnva 
wri te(iwr, 103 ) 


urtiiilttAl. PAGE lb 
OF POtft QUALITY 


NTCA.NNVA, AFTER ADDITION OF THE VELOCITIES", 
CENTERS 1 TO l* OF THE INT.WS") 


write(iwr,10!) Ixnv, (wxnvt(i ,j) ,i»l,nnvr) 
wri te (iwr, 101) lynv, (wynvt (i , j) , i»l,nnvr) 
wr i te (iwr , 101) Iznv, (wznvt (i ,j) , i-l,rnvr) 
continue 


ldeb*>5 

TOTAL INDUCED VELOCITY AT INTERMEDIATE WAKE 

i f (i trans.eq.O) ideb*l 
nrl21“nr 12-1 
nr!22«nr 12-2 
lrl2l-i r 12-1 
i r 1 22-i r l 2-2 
do 1001 i-trU.lrl2 
do 2000 j«2,niva 

rs (i , j) -sqrt (xiv (i , j) ftft2+yiv(l ,j) ftft2) 
if (i .eq. irl 1) gmi^ql 
i f (i .c-q . i r 1 22) gmi-q2 
I f (i .eq. ir 121) gmi“q3 
if (i .eq. ir 12) gmi-ql* 

wzivt (i . j) -wzivr (i , j) + (gmi *wi s) / (4*3 . U 159 rt rs (I , j)) 

wxivt (i , j) -wxivr (i , j) 

wy I v t ( i . j ) »wv i vr ( I , j ) 

cont i nue 

conti nue 


TOTAL INDUCED VELOCITY AT NEAR WAKE 


do 5016 i*nr I 1 , nr I 2 
do 5016 j“2,nnva 

rs (i , j) “sqr t Ur.v (i , j) A*2+ynv(i , j) e*2) 
i f (i .cq .nr 1 1) gni^ql 
i f (i .cq.nr *22) g:n : »q 2 
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5016 
c * 


111 

222 


c 

c 

e 333 

c UUU 

e 

e 

e 

c 

e 


555 

666 


c 


if (I .eq. nr 121) smi<*q3 
If (I .eq.nrll) gmi-qi, 

wznvt (i ,j) «w 2 nvr (i ,j) x (gni*wis) / (i , j ) ) 

wxnvt (I , j) -wxnvr (i , j) 
wynvt (i , j) -wynvr (i , j) 

COntinue v rAV' « 

or 

do 111 i-2.nncrl ' 

'i f (etanc (i) .gt.cgl) go to 222 






i bo 1-1 . «y.« 

pol^Xcgt-otanc (ibol) ) / (otonc (ibol+l) -etanc (ibol) ) 

. Itoo H- i bo 1 - 1 

'wtnvtfnr 1 1. 1) -wzc (lboll) + (w2c (i bo 11+1) -W2C ( « bo 11)} *pol l 
wynvt (nr 1 1 , 1) ®wyc (ibol 1) + (wyc (ibol 1+1) -wyc ( I bo 11)) *po1 1 
wxnvt (nr 1 1 , l)**wxc (i bo 11) + (wxc (ibol 1+1) -wxc (ibol l) ) *pol 1 


do 333 i«2,nncrl 

if (atone (i) .gt.cg2) go to t,i,l« 

ibo2-i 

poI3« (cg2-8tonc (ibo2) ) / (otonc (ibo2+l) -otonc (ibo2) ) 
ibo21-ibo2-l 

wxnvt (nr 122, l) -wzc (Ibo21) + (w 2 c (ibo21+l) -wzc(ibo21)) *pol 3 
wynvt (nr 122, 1) -wyc (ibo2l) + (wyc (ibo21+l) -wyc (i bo2 1) ) *pol 3 
wxnvt (nr 122, 1) -wxc (ibo21) + (wxc (ibo2i+l) -wxc (ibo2l) ) *pol 3 


do 555 i-2.nncrl 

if (atone (i) .gt.cg3) go to 666 

ibo3"i 

pol U m (cg3"Otanc (i bo3) ) / (otanc (i bo3+l) -etanc (ibo3) ) 
ibo31"ibo3-l 

wxnvt (nr 12 1 , 1) -wzc (ibo31) + (wzc (ibo3l+l) -wxc (ibo 3 1 ) ) *poll» 
wynvt (nr 121 , 1) -wyc (ibo3l) + (wyc (ibo 31 +l) -wyc (ibo 3 1 ) ) *poll» 
wxnvt (nr 121 , 1) °wxc (ibo31) + (wxc (1 bo31+l) -wxc ( i bo3 1) ) ft poll, 


do 777 i-2,nncrl 
If (etanc (i) .gt.cgl,) go to 883 
777 I bou™ i 

888 pol5" (cgli-otanc (ibol,) ) / (otonc (ibol»+l) -etanc (ibol,) ) 
iboi,l-ibol,-l 

wxnvt (nr 12, 1) -w 2 c (iboU) + (wzc (i boLl+1) -wzc (ibol»l) ) *pol5 
wynvt (nr 12.1) -wyc (ibol*l) + (wyc (ibol, 1+1) -wyc (ibol,!) ) «pol5 
wxnvt (nr 1 2 , 1) -wxc ( I boM) + (wzc ( i bol, 1+1) -wzc ( i bol,l) ) *po 1 5 

C * 

c * TRANSFORMATION : RECTANGULAR TO POLAR COORDINATES 

c * OF THE VELOCITIES INDUCED AT THE N00ES 

c * 


16 continue 

if (itrcw.cq.2) wrl to (iwr ,1,60) 

A60 format ("1*** WNOFC VELOCITIES AT THE NODES. POLAR COQR.") 

do US j-l,nnva 

1 1- nr 1 1 

1 2- nr 1 2 


7000 do 30 i-i 1 , 12 

if (i .eq.nr 1 2) i 3*0 
if (i .ne.nr 12) i3*l 
wxx-wxnvt ( i , j) 
wyy-wv.nvt (i , j) 
xx-xnv (i , j) 
yy-ynv (i . j) 



r*sqr t (xxAxx+yyAyy) 
wrnvt (I . j ) •« (xx'-'-wxx-t-yyrtwyy) /r 
wtnvt (i , j) - (-wxxrtyy+wyyrtxx) /r 
80 ^ continue 

if (itVcw.nc.2) goto 7001 
wr i tc (iwr , 103) 

wri te (iwr, 101) Irnv, (wrnvt (i , j) , 1*1 ,nnvr) 
wr i te (iwr , 101V1 tnv, (wrnvt (i , j) , !“1 ,nnvr) 
7001 • 1-1 

I 2-onvr 

if(i3.eq.0) go to 7000 
1*5 continue 

do 1*6 j“2,niva 
do 81 i"irll,irl2 


103 


ORIGIN 

of poos QUALtfV 


81 


1*6 


c 

c 

c 

c 

c 

c 

c 

c 

c 

162 

c 

c 


wxx-wxivt (i , j) 

wyy«wyivt (i ,j) 

xx-x i v ( i , j ) 

yy-yiv(i ,j) 

r-sqrt (xx*xx+yyAyy) 

wrivt (i . j) «* (xx*wxx+yy*wyy) /r 

wtivt (i , j) - (-wxx*yy+wyy rt xx) /r 

continue 

if (i trcw.ne.2) goto U6 
write(iwr, 103 ) 

wr i te (iwr , 101) Iriv, (wrivt (i ,j) , i-irl 1, ir 12) 
write (iwr, 101) 1 tiv, (wtivt (i ,j) , l-trll, i r 1 2) 
continue 

wr i vt (i r 11 , 1) -wrnvt (nr 1 1 ,nnva) 
wtivt (i r 1 1 , 1) -wtnvt (nr 1 1 .nnva) 
wr ivt (ir *2, 1) -wrnvt (nr 12. rnv3) 

• wtl vt (i r 12, 1) -wtnvt (nr 1 2, nnva) 

wzivt (i r 1 2, 1) "wznv; (nr 1 2, nnva) 
w2i vt (i r 1 1 , 1) -wrnvt (nr 1 1 ,nnva) 
wrivt (ir 1 21, 1) -wrnvt (nr 121, nnva) 
wt ivt (ir 121 , l) -wtnvt (nr I 21 , nnva) 
wzivt(irl21,i) -wznvt (nr I 2 1 .nnva) 
wr ivt (i r 1 22, 1) -wrnvt (nr 1 22, nnva) 
wti vt (i r 1 22, 1) -wtnvt (nr 1 22, nnva) 
wzivt (ir 122, 1) -wznvt (nr 1 22 .nnva) 

* 

ftftftftftftftftftftftftftftftftftftftftftftftfcftftftftftftftftftftftftftftftftftftftftftftftftft 

A 

* INTERPOLATION 0? THE 'OLD' AND THE 'NEW' DISTRIBUTIONS 

* OF VELOCITIES 

ft 

if (i trcw.no. 0) wr i te (iwr ,1*62) 

format ("lftftft WNOFC ftftft : VELOCITIES AT THE NODES ") 

* 

* WEIGHTING FACTORS FOR THE 'NEW' AND THE 'OLD' VELOCITIES 
f rn«.li 
f tn-. 7 
fzn-.3 

if (ni ter.ge.3) goto 1*9 
frn-.5 
f tn-.8 
f zr.-.5 

if (niter. cq. 2) goto 1*9 

f rn». 5 

ftn=.9 
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f*n».5 

continue 

i f (i trcw.ne.O) wr i te (iwr ,950 f rn, f tn, fzn 

formatC *** WEIGHTING FACTORS: FRN.FTN.F2N 3f 6-2) 

frv«l .-frn 

ftv»l .-ftn 

f2V-l.-f2n 


*\ 

v. w ' 




do 47 j • 1 . nnva 

1 1- nr I ! oO 0 ^ 

12- nrl2 0* 

7003 do 43 I "i 1 » • 2 

If (i .eq.nr!2) i 3*0 
if (i .ne.nr!2) i3“l 

wrnv(i , j) “wrnvt (i , j) *frn+wrnv (i , j) *frv 
wtnv(i , j) “wtnvt (i , j) *ftrt+wtnv(i , j) *ftv 
wznv(i , j) «wznvt (i , j) *f 2 n+wznv (J , j) *fzv 
go to 43 

648 wrnv(i , j) -wrnvt (i , j) 
wtnv ( i , j) -wtnvt (i , j) 

W2nv(l ,j)«x2nvt(i ,j) 

43 continue 
If (I trcw.eq.Q) goto 7002 
wri te (Iwr , 103) 

write (iwr, 101) Irnv, (wrnv(i ,J) , i»l,nnvr) 
wr I te (Iwr , 101) 1 tnv, (wtnv (i , j) , i«l ,nnvr) 
wri te(iwr, 101) Iznv, (wznv(i , j) , i»l,nnvr) 

7002 i x-i 

12*nnvr 1 

if(l3-eq.0) go to 7003 

47 continue 
do 48 j“l,niva 
do 44 i-irll,irJ2 

wr i v (1 , j) **wr i vt (i , j) *f rrt+wr I v (i . j) *f rv 
wtiv(i , j) ®wtivt (i , j) *f tn+wtiv (i ,j) *f tv 
wriv(i t j) **W2 i vt (i ,j) *f2n+w2iv(i ,j) *fzv 
go to 44 

649 wriv(i ,j)«wrivt(i ,j) 
wtiv (i , j) "wt i vt (i , j) 
wziv (i , j) -wzivt (i , j) 

44 continue 
if (itrcw.eq.0) goto 48 
wr i te (iwr, 103 ) 

write (iwr, 101) Iriv, (wriv(i.J) , i»ir 1 1, 1 r 12) 
write (iwr, 101) I tiv, (wtiv(i , j) , i-ir 1 1 , ir 1 2) 
wri te (iwr, 101) Iziv, (wziv (i , j) , i»ir 1 1, ir 12) 

48 continue 
c * 

c ***************************************************** 

C * 

c * INTERPOLATION OF THE INFLUENCE COEFFICIENTS 
c * 

if (ni ter .cq. 1) return 
fzn«.5 
f zv**l .-fzn 

i f (i trcw.eq.0) goto 90 
wr i te (iwr , 500 ) 

500 formatC *** WN0FC: INTERPOLATION OF INF. COEFF 

S" INF. COEFF FROM L00P2&") 
do 91 j*l ,nncr2 
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wrl te (iwr , 101) 1 twy, (twy (i , j) , i-1 . nncr2) 
wr i te (iwr , 101) 1 twz, (twz (i , j) , i-1 , nncr2) 
wr i tc (iwr, 103 ) 
continue 
wr i te (iwr , 501 ) 
format (" OLD INF. COEFF. ") 
do 92 j-l ,nncr 2 OKiGINAL. PAGs 15 

wr i te (iwr, 101) Itwy, (twyt(i.j) , i«l,nncr2) qF POOR Qil/tJ-tW 
wr i te (iwr , 10 1) 1 twz, (twzt (i , j) , i-1 ,nncr2) 
wr i te (iwr , 103 ) 
continue 
continue 

do 32 i«l,nncr2 
do 82 j-l,nncr2 

twy (i , j) -fzn*twy (i , j)+fzv*twy t (i , j) 
twz (i , j) -f zn*twz (i , j)+fzv*twzt (i , j) 
continue 

if (i trcw.eq.0) return 
write (iwr, 502 ) f 2 n,fzv 

format(" IN7. INF. COEFF.: FZN-",e20.5,‘‘ FZV-",e20.5) 
do 93 j-l,nncr2 

wr i te (iwr , 101) i twy, (twy (i , j) , i-1 ,nncr2) 

wr i te (iwr , 101) I twz, (twz (i , j) , i-1 ,nncr2) 

wr i te (iwr, 103 ) 

continue 

return 

format(lx.ai»,lx,5fl3.9.'/.6x.5fl3.9./.6x,5fl3.9) 
format (lx, aL. lx. 9 • 13./.6x,9i 13) 
format (lx) 



subroutine i nt;r2 
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» f. 

c 

e 

c 

c 

c 

c 

e 

c 

c 

c 

c 

c 

c 


a 

n 

ft 

ft 

« 

ft 

ft 

ft 

ft 

* 

* 

A 

ft 

ft 


13 

c 

c 

c 

c 


92 

91 

c 

c 

c 

c 

c 


ftftftftftftftftftftftftftrtftrtftftftftft* 
ft ft 

ft PROGRAM fWC * 

ft ft 

ft SUBROUTINE 1NTGR2 ft 
ft ' • ft 

ftftftftftftftftftftftftftftftftftftftftftft 


p|\G£ 




03/29/78 

INTEGRATION OF THE VELOCITIES AT THE NODES , 

COMPUTED BY WN0FC2, TO GET THE NEW POSITIONS OF THE NODES 

cc-wvon/al I/cgt,cg2,cgJ,cgi., Ir 

consnon/parm/'lwr , Ird, 1 trace, Ipscst.nffles, I trcw, I trcl 1 , 1 trcll, I trcg, 

& I trcnt , i trcf , i plot, I cgsn, i pi otv, i trtg, Ivergr , i test, I saoeb, i samet, 
(I ins, I ini. I ira2.nl ter, i trans, 1 pn.j punch, i plotw, i sw, Scant, iplotg, 
tlpr, iplott 

coffswn/gaoa/nbldsl.nblds, a i gs*a. fru.et an (25) .knnvr.etai (6) .knivr, 

H twist. thetad(25) ,thcta(25) ,thetac(20 , thetaO, thelOd, alphas, 
ScdO.cdk.dpsInd.dpsin.dPsI id.dosl i.cosf f .coeff l,c,s, blades, 
innvr.nnva.nnvr l.nnval.nncr .nnca.etanv (25) .etanc (26) , 

An Ivr .niva.nl vr l.nival.nlcr.nica.etaiv (6) , etaic (7) , 
tntva.ntva l , ntea, fpt 1 , f pt2 , nr 1 l.nrl2, Ir 1 1, i r 1 2. Iftcdel 
cession /resul/gaisrac (26) ,twx (2b, 2b) . twy (2020 , i* 2 .(lU, 20 , 

Swxc (20 ,wyc (20 ,wjc (21.) ,wxnc (26, 19) .wync i26. I2f .wsnc (26. 19) • 
lwxic(9.5U ,wylc(9.50 ,waic(5,5l) ,wrnv(28. IS? ,wtnv(2S. 18) . 

Sw 2 nv (25 , 13) ,wr iv(8.«0) ,wtiv(3.50) ,w 2 lv (E.=0) , 

txnc (26. IS) ,ync(26.i3) ,xnc(26.iS) ,xic(9.5i) ,yic(9.50 .rJc(9.5S). 

6xnv(23.18) .yny(23.IS) ,2nv(23,13) ,xiv(8.50 ,ylv(3.S0) ,2iv(3o0). 

& i a , i b 

cosmon /wdata/xt (2S) .yt(23) ,:t (28) .r In (23) ,vrln(28) ,vt!n(28) , 

Av 2 l (28) ,vx2(23) ,vy2 (23) , va2 (2S) 
equivalence (nnvr ,nncr 1) .(nivr.nicrl) , (nnva, nnca 1) 
equivalence (niva.nicel) , (nnvr l ,nncr2) , (nicr l,nicr2) 

ft 

etam-.jft (etanv (nnvr) +etanv (l) ) 

ft 

xep»0. 

do 13 xpt«l,nncr2 
xep-xepfwac (kpt) 
xcp«*xep/f I oat (nncr2) 

ftftftftftftftAftoftftxitftftfteeefteftftftM 

ft 

ft VERIFICATION OF THE INTEGRATION 


lsav»i trtg 

1 f (ivergr, ne. 1) go to gi 

1 trtg-l 

goto 91 

Conti nue 

ivergr-0 

cont i nue 

ft 

Attftrtftftftftftftftftftftrtftftftftftftftrtftftft 

ft 

* LOOP ON THE TWO WAKES 

ft KW-1: INTEGRATION ?0R THE NEAR WAKE 
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c * 
c « 


C ft 
C ft 

c * 

c 

c 


c 


U2 


c * 

C ft 

60 


c 


52 

c * 


C ft 

c ft 

61 
100 


KW-2: INTEGRATION FOR THE INTERMEDIATE WAKE 

do 90 Aw* l, 2 
If (Lw.eq.2) goto 60 

initialisation FOR THE NEAR WAKE 

nr 1 2 l*>nr 12*1 
nr 1 22**nr 1 2-2 
• Ir I21*i r 1 2- I 
I r 122-lr 1 2-2 
do 1*2 i«nrl l, nr 12 
If (I .eq.nr!2) r 1 n ( I ) "cgi* 

If (I .eq.nrl21) rln{l)*cg3 
I f (I .eq.nrl 1) rln(i)«cgl 
if (i .eq.nr 122) rln(i)*cg2 
vr In (i) «wrnv (i . 1) 
vt In (i ) «wtnv ti , l) *r In (i) 

V2 1 (I ) »wrnv (i . 1) ♦fmu 
xnv(i , l) «r in (I) 
ynv (1 . 1) *0. 

2 nv (I , 1) «0. 
xt (i) «r In (I) 
yt (i) *0. 

2 1 (i) -0. 
contl nuo 
dpsi»dpsin*.5 

dsr*3 . ftobs ( (xop-rf.-nu) *dp* i n) 

nva*nnva 

nvr*nr 12 

goto 61 

initialisation tor the intermediate wake 

continue 

do 52 !■* r 1 1 . Ir 1 2 
I f (I .eq. ir 1 1) i l»nr I 1 
1 f (i . eq. i r 1 2 1) i J*nr 1 2 1 
if (i .eq. i r 1 22) Il«nrl22 
if (I .eq. irl2) I l*nr 12 
xlv (• , t) *xnv (i 1 ,ntva) 
y iv (I , l) *ynv (i l , ntva) 

2 1 v (I , l) «2nv (i l.ntva) 

xt (I) *xiv (1 , 1) 
yt (I) “yiv(» , 1) 

2t (I) * 2 lv (i . 1) 

r ln(i) -sqrt (xiv (i . 1) ft>'2+y i v (I , 1) AA2) 
vr In ( i ) ~wr i v (i , 1) 
vtln(i)«wtiv(i,l)+rln(i) 

V2 t (• ) >w: iv (i , 1) + .'mo 
cont i nue 
dpji®apii i ft. 5 

dir«*3. l (xcdi fmu) ftdpsi i) 

nvo a niv,i 

nvr» i r I 2 


cont i nue 

fornot ("Iftftft IN7CR ftftft; KW«",il," NVA- , ',i2," S VR*" 


UKIGINAL PAGE 

of poor? QUAUT? 


ik«", i2. 



If 

S" OPSI f 3 - 2 ) 

c * 

e * THE RAOUL AND TANGENTIAL VELOCITIES ARE WEIGHTED BY 
C * RNEW/ROLD . 

C * ALL THE JOINTS ARE CONSTRAINED TO 0 E IN THE PLANE (P) 

t ft WHICH CONTAINS THE NEW °0 I HT OF THE REFERENCE STREAML 

C ft ANO THE PROJECTION OF THIS POINT ON THE Z AXIS 

€••'*' IN ORDER TO LIMIT THE OlSTORSICNS OF THE ELEMENTS 

e * 

c * LOOP ON THE AZ I MITTAL POSITIONS 

c * 

do 81 J" 2 ,nva 

c * 

c * UNCORRECTED NEW POINTS (FIRST LOOP ON THE RAOUL POSITIONS) 
if(kw.eq.J) lj-nrll 
lf(kw.eq.2) I j«* I p 1 1 

C * 

do 82 l-lj,nvr ’ 

I f (kw.eq.2) goto 711 
C * OLO RADIUS OF THE POINT 

r2o**sqrt (xnv(l , j) **2-*-ynv(l , j) **2) 

C * RADIAL VELOCITY AT THE POINT (UNCQRRECTEO) 
vr2o«"wrnv(i , j) 

C * TANGENTIAL VELOCITY AT THE POINT (UNCORRECTEO) 

vt2o®wtnv(i ,j)+r2o 
C * AXIAL VELOCITY 

vz2 (1) -wznv (i , j)+fnu 
go to 712 

711 r 2 o»sqrt(xlv{i ,j)** 2 +ylv(l ,J) ** 2 ) 
vr2o-vr iv(i , j) 
vt2e-wtiv(i ,j)+r2o 
V22 (I)-wziv (I , j)* 1 frsu 
712 continue 

c ft 

c A NEW RADIUS OF THE POINT: 

C * MEAN RADIAL VELOCITY:- .5* (VR 1 H-*VR 2 N) 

C * VR 1 M: VELOCITY AT THE PRECEDENT POINT (CORRECTED) 

C * VR 2 N: VELOCITY AT THE CURRENT POINT (COREECTEO) 

C * VR 2 M:«VR 20 *R 2 M/R 20 

r 2 r»- (r In (I) ■•■dps i *vr In (i) ) / (I .-dps I <*vr 2 o/r 2 o) 

C * CORRECTIONS ON THE RADIAL ANO TANGENTIAL COMPONENTS 

C * OF THE VELOCITIES 

ar-r2n/r2o 

! f (at.gt. 1 .05) r 2 n-r 2 o* 1.05 
If (at.gt.1.05) at- 1.05 
vr2n-vr2o*at 
vt2n«vt2o*at 

C * NEW AZIMUTH (TEMPORARY) 

th 2 «atan 2 (yt(i) ,xt (i))+dpsU 2 .A(vtln (i) +vt 2 n) / (rln (i)+r 2 n) 
cc-cos (th2) 

*S“* i n (th 2 ) 

C ft RECTANGULAR COORDINATES Of THE UNCCRRECTED NEW POINT 

xt (i ) -r2n><cc 
yt (i) “r2n*js 
Ztd-dpsi A (vz 2 ( » ) -*-V2 1 {;)) 
i f (abs (ztd) .gt.dzr) ztd-sign (dzr , ztd) 
zt(I)-zt(i)+zto 

C * VELOCITY AT THIS POINT 

vx2 (i) -vr2n»cc-vt2n*t i 
vy2 (i) *vr2n*'i+vt2n’>cc 
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c * 
c * 


82 ' 
c * 


721 


83 

81 

90 

c * 


c * 


c * 


SAVE RADIAL AND TANGENTIAL VELOCITIES. AHO RADI 
WILL BE USED FOR the NEXT azih.uthal POSITION 
vtln(i)-vt2n 
vrtn(i)-vr2n 
vzl (i)«vj2 (i) 
rln(i)»r2n 
continue 

do 83 j .nvr 

If (kw.eq.2) goto 721 

xnv(i ,j)-xt (i) 

ynv(i,j)»yt(!) 

znv (i .j)-zt(i) 

goto 83 

cont i nue 

xiv(i , j) «xt (i) 

ylv(i ,j)»y: (i) 

zlv(i ,j)- 2 t (i) 

continue 

continue 

continue 

i f (ivorgr .eq. 1) goto 92 
I trtg"i oav 

return 

end 




AL COORDINATE 


UHIGINAL PAGE 13 
OF POOF? QUALITY 


A 


c 

c 

c 

e 

e 

c 

e 

c 

c 

c 

c 

e 

e 

c 


« 


c 

c 


13 

c 

c 

c 

c 


92 

91 

c 

c 

c 

c 

c 


subroutina intgr3 

* 

* ftftnftftftftrtftftftrtftftrtftftftftftftft 

A A A 

* A PROGRAM FWC * 

ft * * 

ft ft SUBROUTINE IN73R3 * 
ft ft ft 

ft ftftftftftftftftftftrtftftftftftftftftftftft 

* 

* ' 03 / 29/78 

ft 

ft INTEGRATION OF THE VELOCITIES AT THE NODES , 
ft COMPUTED BY VN0FC3. TO GET THE NEW POSITIONS OF THE NOOES 

ft 

eatmon/al I /eg 1 .CQ2.cg3.cg!*. I r 

n l * ' eoraon/parn/lwr, i rd, I trace, Ipsen, renss, I trew, itrcl l, I trc12, I treg, 
t 61 trent . I tref . I plot. iegen, iplotv, i trtg, ivergr , I test* i saraeb, i samet, 
&l Iras, 1 ini. 1 in2.nl ter, i trans, I pn, jpunch, iplotw, i tw, Iconc, iplotg, 
&ipr , Iplott 

ccJKnon/geesn/nbldsl.nblds, sigma, fmu.ctan (25) .knnvr.otal (6) .knivr, 

&l twi it, rhotad (25) , thota (25) , thotac (2U) , thctaO, thatOd, alphas, 
ScdO.cdb.dpsInd.dpsin.dpsi id, dps 1 1 .cocff ,coof f 1 ,c,s, blades, 
*nnvr,nnva,r,nvrl,nnval,nncr,nnca,etanv(25) ,etanc(26) , 
tnlvr.niva.nlvr l, nival, nlcr.nl ca, eta iv (6) , etaic (7) , 
Sntva,ntval,ntca,fpsl,fp»2,nr I l.nrl2, ir 1 1, ir)2, incdel 
ce.T,non /resul/garvr.c (26) ,twx(2b,2b) . twy (2b, 2b) , twz(2b,2b) , 

6wxc (2b) ,wyc (2b) ,wzc (2b) .wane (26, 19) .wync (26. 19) .wane (26. IS) » 

Cwxic (9.51) .wy!c(9.51) ,w=ic{5.51) ,wrnv(23,l8) ,wtnv(28.18) . 
&wznv(23, 18) ,wrlv(3,50) ,wtlv(8,50) ,«*xiv (8.50) , 

6xnc(26.1S) ,ync (26.19) .me (26. 19) ,xic (9.5 1) .y ic (9.515 ,zic(9.5l) . 
*xnv (23, 18) .ynv (23.13) ,2nv(28.l8) .xlv(8,50) ,yiv(8.50) .ziv(3.50) , 
Sta.ib 

ccsnon /v,data/xt (23) . y t (28) ,zt(23) ,rln(23) ,vrin(2E) ,vtln(28) , 
tvzl (28) .vxi (23) ,vy2 (28) .vz2 (23) 
cquivalonca (nnvr.nncrl) , (nivr.nicrl) , (nrtva.nncal) 
equlvalcncft (niva.nical) , (nnvrl,nncr2) , (nicr l,nicr2) 

* 

#tora».5ft (etanv (nnvr 1) +etanv (l) ) 

* 

xcp*0. 

do 13 kpt«l,nncr2 
xop“xep*-wzc (kpt) 
xep«xep/f Icat (nncr2) 

ftftftftftftftftftftftftftftftnftftftftftftftftftft 

ft 

ft VERIFICATION OF THE INTEGRATION 

ft 

ltov-1 trtg 

If (ivergr.ne.l) go to 91 

I trtg»l 

goto 91 

cont i nue 

ivergr»0 

continue 

ft 

A 

* LOOP ON THE TWO WAXES 

ft KW-1: INTEGRATION FOR THE NEAR WAKE 
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, 



■***1 


C ft KW-2: INTEGRATION FOR THE INTERMEDIATE WAKE 
c ft 


c 

c 

c . 


<42 

43 

41 


50 


c 

c 


60 


do 90 kv,-1.2 
.1 f (kw.eq.2) goto 60 
ft 

ft INITIALISATION FOR THE NEAR WAKE 
ft 

IrI2l-irl2-l 

nr!21-nrl2-l 

* do 42 » -nr 1 l, nr 12 

rt If (1 .eq.nr 12) rln(i)«cg4 
If (I .cq.nrl l) rln(i)-cgl 
if (I .eq.nr 121) rln(i)«cg3 
vrln(i)-wrnv(i ,1) 
vr In ( i) -wtnv (i , i)+r In (i) 
vzl (i) «wznv(I , l)+fmu 
xnv (1 , 1) »r In (i) 
ynv(i,l)»0. 

2 nv(l , 1) -0. 

xt(i)-r ln(i) 

yt(i) -0. 

zt(l) “0. 

continue 

do 43 l“l,nnvrl 

I f (etanv (i) .gt.etam) go to 41 

I k-t 

I f (etanv (ik+l) -etam. I t.etam-etanv (ik)) 

do 50 i*l,nnvr 1 

rln(i)-etanv(i) 

vrlntD-wrnvU , I) 

vtln (i) *wt.iv (i , U+etanv (i) 

V2l (i)-wznv(i , lf+fmu 

xnv(i , 1) -etanv (i) 

ynv(i , 1) “0. 

znv (1 , 1) «=0. 

xt(i)-etanv(i) 

yt(i)-0. 

zt (i) *0. 

continue 

dpsi-dpsinft.5 

dzr-3.ftabs ( (xep+fmu) *dpsin) 
nva«nnva 
nvr-r.r 1 2 
goto 61 
ft 

* INITIALISATION FOR THE INTERMEDIATE WAKE 

continue 

do 52 i-irl 1, I r 1 2 
I f (1 . eq . i r 11) i l«nr 1 1 
i f (i .eq. ir 1 21) i 1-nr 121 
1 f (i .eq . ir 1 2) i 1-nr 1 2 
xlv(i , 1) »xnv(i l.ntva) 
y i v (i , 1) -ynv (i 1 , ntva) 
ziv (i , 1) -znv (i l.ntva) 


ik-ik+l 
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xt (i) «xiv(i , 1) 
yt(i) a yiv(i , l) 

Z t (i ) -z i v (i . I) 

rln(i)»5ort(xiv(i, 1) ft*2+yiv(i , 1) >i*2) 
vr In (!) »wr i v (i , 1) 
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-"T" 


/ 


52 

c * 


c * 
c ft 
61 
100 

c ft 
c * 

C ft 
C ft 
C ft 
C ft 
C ft 

c * 

C ft 

c * 


C ft 

C ft 


C ft 
1000 


c * 
C ft 

c * 

C ft 


712 
c * 
C ft 

c * 
C ft 

c * 

C ft 

C ft 
C ft 


v t In (l > »wt 1 v (1 , l)+r|n(l) 

V2 1 (I ) -wzi v (i , l)+fmu 
continue 
dpij-dp*i ift.5 

d2f"3**a&* *dp*i i) 

nva-niva j 
• ;nvr»lrl2 




,g& 




o? 


?G0^ 


continue 

format ("Iftftft INTGR ftftft: KW-",!1," NVA-",I2," NVR«",I2," IK»",I2, 
A" DPSI-",f3.2) 

THE RADIAL ANO TANCEMTIAL VELOCITIES ARE WEIGHTED QY 
RNEW/ROLD. 

ALL THE POINTS ARE CONSTRAINED TO BE IN THE PLANE (P) 

WHICH CONTAINS THE NEW POINT OF THE REFERENCE STREAMLINE (I K) 

ANO THE PROJECTION OF THIS POINT ON THE Z AXIS 
IN OROER TO LIMIT THE OISTORSlONS OF THE ELEMENTS 

LOOP ON THE AZIMUTAL POSITIONS 

do 81 j"2,nva 
I 2-nvr 
I3«0. 

UNCORRECTEO NEW POINTS (FIRST LOOP ON THE RADIAL POSITIONS) 
If(kw.eq.l) ll«nrll 
if(kw.eq.2) il-irll 

do 82 1-11,12 
If (l2.eo.nnvri) 13-1 
lf(kw.eq.2) goto 7U 
OLO RAO I US OF THE POINT 

r2o-sqrt (xnv(l , j) ftft2+ynv (1 ,j) **2) 

RADIAL VELOCITY AT THE POINT (UNCORRECTED) 
vr2o-wrnv (I ,J) 

TANGENTIAL VELOCITY AT THE POINT (UNCORRECTEO) 
vt2o-wtnv(l ,j)+r2o 
AXIAL VELOCITY 

vz2 (l) -wznv (I , j)+fmu 
Co to 712 

711 r2o»sqrt (xiv(I , j) ft*2+y Iv (I ,J) ft*2) 

vr2o-wr iv (i , j) 

vt2o-wttv (I , j)+r2o 

V22 (I) -wzlv (i , j)+fmu 

continue 

NEW RADIUS OF THE POINT: 

MEAN RADIAL VELOCITY:- -5* (VR1N+VR2N) 

VR1N: VELOCITY AT THE PRECEDENT POINT (CORRECTED) 

VR2N: VELOCITY AT THE CURRENT POINT (COREECTEO) 

VR2N;-VR2Q rt R2N/R20 

r 2n- (r in (i)+dpsiftvr In (i) ) / (l .-dps i ftvr2o/r2o) 

CORRECTIONS ON THE RADIAL ANO TANGENTIAL COMPONENTS 
OF THE VELOCITIES 
at-r2n/r2o 

I f (at.gt. i .05) r2n-r2o* l .05 
If (at.gt.l.05)at-l.C5 




vr2n«vr2oAat 

vt?n=vt2o*at 


C * NEW AZIMUTH (TEMPORARY) 

th2-atan2 {yt (i) ,xt (i) )+dpsi*2.A (vtln (i)+vt2n) / (rln (i) +r2n) 
cc*cos(th2) '■» • 

ss*sin(th2) * ' *'* 

c * RECTANGULAR COORDINATES Of. THE UNCORRECTED NEW POINT 

xt (i) -r2n*cc „ 

yt(l)-r2n*ss OHIGINAV PAG , E JS 

2 td*dps i * (vz2 (i)+vzl (I)) OF PtJOH QUALITY 

If (abs (z td) .gt.dzr) ztd-sign(dzr.ztd) 
zt (I) *zt (i ) +ztd 

C * VELOCITY AT THIS POINT 

vx2 (i) "vr2nftcc-vt2rv>55 
vy2 (i) *vr2n*ss+vt2n*cc 

C A SAVE RADIAL AND TANGENTIAL VELOCITIES, AND RADIAL COORDINATE 

c A WILL BE USED FOR THE NEXT AZIMUTHAL POSITION 
vt In (!) -.vtZn 
vrln(i)*vr2n 
vz I (I)*vz2 (i) 
r In (i ) —r2n 

82 continue 

if (I2.eq.nnvrl) go to 1001 
go to 1003 

1001 xO-xt(ik) 

yO^yt(ik) 
zO*zt(ik) 
vx*vx2 (ik) 
vy*vy2 (ik) 
vz*V22 (i k) 

alpha*- (x0*vx+y0*vy) / (xOAxO+yOAyO) 

. vx^/x+alpnartxO 

vy»vyj-alpha*yO 

do 85 i*i 1, i2 

beta* ( (xO-xt (i) ) *vx+(y0-yt (') ) ! '‘vy+ (zO-zt (i) ) nvz) / 

& (vx2 (i) *vx+vy2 (i) *vy+vz2 (i) >'<vz) 
xt (i) *xt (i)+beta*vx2 (I) 
yt (1) *yt (i)+beta*vy2 (1) 
zt (i) *zt (i)+beta>' [ vz2 (i) 
r In (i) -sqrt (xt (i) *>'2+yt (i) **2) 
vzl (i) *V 22 (J) 
xnv (i ,j) -xt (i) 
ynv (i , j) *yt (1) 
znv(i ,j)* 2 t(i) 

85 continue 

go to 81 

1003 do 83 i-11,12 

If (kw.eq.2) goto 721 
xnv (i , j) *xt (i) 
ynv (! , j) *y t (i) 
znv (i , j) *zt (i) 
goto 33 

721 continue 

xi v (i , j) -xt (i) 
yiv(i .j)-yt(i) 

Ziv (i , j) -zt (i) 

83 continue 
i 1-1 

i2*nnvrl 

if (i3.nq.0.and.kw.eq. 1) gc to 1000 



conti nue 
continue 
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If (ivergr.eq.l) goto 92 
ltrtg°isav 

return 

f . 

end '• ' 






subroutine loop!2 
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ftftftftftftftftrtftftftftftftftftftftftftft 


ft ft 

* PROGRAM FWC A 

ft ft 

* SUBROUTINE LOOP12 ft 

ft ft 


ftftftftftftftftftftftftftftftftftftftftftft 


CRISlf&U PAG* ?3 
OF POOR QUflJXTV 


07/U/80 

EVALUATION OF THE INDUCED VELOCITIES AND OF THE DISTRIBUTION 
OF CIRCULATION ALONG THE BLADES, 

FROM THE INFLUENCE COEFFICIENTS RETURNED BY L00P2 


THE SUBROUTINE RETURNS: 

ITEST-0 IF CONVERGENCE IS REACHED 

ITEST-1 IF CONVERGENCE IS REACHED FOR THE LOOP ON THE CIRCULATION 
ITEST-2 IF CONVERGENCE IS NOT REACHED WITHIN LIM1 ITERATIONS 
ON THE CIRCULATION 

TWX(I.J) CONTAINS THE DERIVATIVE : 

0 (WX(l))/0 (GAMMA (J+l)) 

TWY: SAME FOR THE Y COMPONENT 
TWZ: SAME FOR THE Z COMPONENT 

NOTE : THE SUBSCRIPTS I AND J ARE REFERED TO THE INNER 

CENTERS OF THE ELEMENTS 

common/rol lp/wxctr (2=) ,wyctr (25) .wzctr (25) 

common/rql I i0/wr.ctr2 (25) ,wyctr2 (25) ,wzctr2 (25) 

common/parm/iwr , i rd, i trace, Ipsem.r.mes, I trcw, I trcl 1, 1 trcl2, 1 trcg, 

& i trcnt, i trcf , iplot, icgen, iplotv, i trtg, ivergr , i test, isomeb, I same t 
61 ims, I • *n 1 , 1 1m2,nl ter , 1 trans, i pn, j punch, i plotw, I sw, 1 cent, I plotg, 
6ipr, Iplott 

common/geom/nbldsl.nblds, sigma, fmu.etan (25) .knnvr.etai (6) ,knivr, 

SI twist, the tad (25) .theta (25) ,thetac(21*) , thetaO. thetOd. alphas, 
6cdO.cdlc.dps i nd.dps i n.dps i id, dps i I ,coef f , cocf f l ,c, s, blades, 
6nnvr,nnva,nnvrl,nnval,nncr,nnca.etanv(25) .etanc (26) , 

6nivr, niva.nivrl, nival, nicr,nica,etaiv(6) ,etaic(7) , 

Sntva.ntval ,ntca, f psl , f ps2 ,nr 1 1 ,nr 12, 1 r ! 1, 1 r 1 2 , i model 
common /resul/gammc (26) ,twx(2£»,2it) ,twy (21*, 21*) , twz (2I», 24) , 

Swxc ( 2 U) ,wyc (21*) ,wzc (21*) ,wxnc (26. 19) .wync (26, 19) ,wznc (26. 19) . 
Cwxic (9,5») ,wy ic (9,51) ,w2ic (9,51) ,wrnv (28. 18) ,wtnv (28, 18) , 

Swznv (23, 13) ,wriv(8,50) ,wtiv(8,50) ,wziv(8,50) , 

Sxnc (26, IS) ,ync (26, 19) .znc(26, 19) ,xle(9.50 ,V* c (3.51) ,zic(9.5D . 
Sxnv (25, 18) ,ynv(28, 185 ,znv(28, 18) ,xlv(8,50) ,yiv (8,50) ,ziv(8.50) , 

6 i a , i b 

dimension wyct (21») ,w3ct(21*) 
dimension w/ett (2U) .wzett (21*) 
dimension gttmp(26) 
dimension truest (21*) 

cqui va I ence (nnvr , n.ncr 1) , (ni vr , n i cr 1) , (nnva , nnca 1) 
equivalence (niva.nical) , (nnvr 1 , nncr2) , (ni vr l , ni cr2) 
data pi 2/9.369601*1*01,/ 
data pi/3. W*15S2o53/ 

data lx. iy , !z/"WXC ","WYC ","VZC "/ 
data lg. Jyt. !2t/"GMC ","WYCT", M WZCT"/ 
data lts/"TE3V7 



c 


ft 


I test-0 
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do 80 > l-l,nner 
gtcmp(h+l)“ga'twnc (i + 1) 
wyctt (i)-wyc (i) 

W2Ctt (i) -W2C (i) 

continue . 
gtcmp (npcr'O “0.0 

* •* 

do 90 i ter»l, I iml 
ft 




* EVALUATE THE DISTRIBUTION OF CIRCULATION AND EPS (TEST VALUE) 


eps-O. 

nncr3-nncr2-l 


TIP LOSS FACTOR 


flammc (nncr 0-0.0 
do l I«»l,fincr2 
cps-eps+wseO) 

tltm-(w 2 c(i)+fmu)/'(etanc(l+l)+wyc(l)) 

1 1 amda-atan (t I am) 

If (flarada.lt. flanda— .15 
u- (etane ( i + l) +wyc ( I ) ) /cos (f I amda) 
alp-f larada-thetac (i) 

I f (a&s (alp) .gt.alphas.and.alp.gt.O.) alp-alphas 

If (sbs (alp) .gt. alphas. and .alp. 1 1 .0 .) alp— -a I phas 

If (i .eq.nncr2) go to l 

gansac ( i+1) -p i *p i ft l . ft S I gma*u*a I p/b I ades 

continue 


c TIP LOSS FACTOR 

c 

c * TEST VALUE: .01 TIMES THE MEAN Z INDUCED VELOCITY 
eps-.02*abs (aps) /nncr2 
epp-l .S^eps 

epp2-3.fteps 

c ft 


ktest-l 
nncr3-nncr2-l 
do 2 i-I,nncr2 
wyct (i) -0. 
wact (i) -0. 
wactr (I) -0. 
wyctr (i) -0. 
do 3 J-l,nncr2 

wyct (i) -wyct (i)+twy (i . j) ftganrac (j+1) 
wzct (i) -wzct (i)+tw 2 (i , j) agamic (j+1) 

3 continue 

2 continue 

call 1 oopO 

do 1000 i-l,nncr2 

wyct (i) -wyct (i)+wyctr (i) 

1000 wact (i) -w2Ct (i)-rw2ctr (i) 

do 1100 i-l,nncr2 

jf (abs (wyct(i)-wyc(i)) .gt.epp) ktest-0 
if (abs (wzct(i)-wzc(i)) .gt.epp) ktest*0 
conti nue 



1100 
c ft 


if (itrcl l.eq.O) goto 60 A1/ 

wr 1 tc (iwr, 100} I ter 

100 format {/" *** L00P1 ***: NITER-", i 2) 

wr i te {iwr, 101) ly , (wyc (kk) ,kk-l ,nncr2) 
wr I tc (iwr , 101) 12 , (wzc (kk) ,kk«l ,nncr2) 
wri tetiwr, 101) Ig . (gammc (kk) ,kk-2.nncr l) 
wr i te (iwr , 101) iyt, (wyct (kk) ,kk-l,nncr2) 

write (iwr, 10!) Izt. (wzct(kk) ,kk-l,nncr2) ORIGINAL PAGE IS 

write (iwr, 150) eps.epp OF POOR QUALITY 

150 format ("eps-",f 10.5, 10x,"epp- ",f 10.5) 

101 format (lx, ai,, lx,6f 10.6,/. 6x,6f 10.6) 
c * 

60 continue 

c * 

c * CONVERGENCE IS REACHED: GOTO 50 
if (ktest.eq. 1) goto 50 

c * 

c * NEW APPROXIMATION OF THE INDUCED VELOCITIES 

c * 

facn-.2 
facv-l .-faen 

c * 

do 6 i-l,nner2 . 

wyc (i) -wyc (i) *facv+wyct (i) *facn 
W 2 C (i) -wzc (i) *facv+wzct (i) *facn 
6 continue 

90 continue 

c * 

c * NO CONVERGENCE WITHIN LIM1 ITERATIONS 
c * 

i test- 2 

• if (1 iml.eq. 1) itest-1 

goto 25 

c * 

C * 

50 continue 

if (1 iml.eq. 1) go to 25 

c * 

c * TEST OF CONVERGENCE FOR THE OUTER LOOP: 

C * THE CONVERGED VALUE OF THE INDUCED VELOCITIES 

C * MUST BE WITHIN (EPS) OF THE ORIGINAL INDUCED VELOCITIES 
c * 

ktest-l 

if (Itrrl 1. cq . 1 ) wr i to (iwr , 150) eps 
do 51 i**l,nncr2 
mtest (i) -0 

if (abs (wyc (i) -wyett (i) ) . 1 1 .epp2) go to 52 

mtest(i)-l 

ktest-0 

52 continue 

i f (abs (wzc (i) -wzett (i) ) . 1 1 .epp2) go to 5‘ 

mtest (i ) -mtest (i) +2 

ktcst-0 

51 continue 

i f (i trcll .ne.O) wr i te (iwr , 1 10) 1 ts , (mtest (kk) , kk”l ,nncr2) 

1 10 format ( !x , ak , lx. 3 i 10 ,/, 6x, 9 i 10) 

i f (ktest.eq. 1) goto 25 
i test- 1 


C 
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c ft MEW APPROXIMATION OF THE INDUCED VELOCITIES 
c ft ANO OF THE DISTRIBUTION OF CIRCULATION 
C ft 

facn*. 5 

If (niter. le. 2) facn*.l 

facv-l .-facn 

if (I trcl 1 .ne.O) write (iwr. 200) facn 
200 format (" ft*ft L00P1 ft** OLO/NEW FACT: FACN-", f 10.5) 

i\* do 91 i"2,nncrl 

1 gammc {i)*facn*gammc (i)+facv*gtemp(i) 

91 *'* * continue 
c TIP LOSC> FACTOR 



garonc(nncrl)*0.0 

c 

If (i trcl l.eq. 1) wr i te (iwr, 101) 1g, (gammc (kk) ,kk-2,nncr 1) 

do 7 i“l,nncr2 

wyc(i)-Q. 

wzc (i) "0. 

wxctr (i) *0. 

wyctr (i) -0* 

wzctr (i) *0. 

do 8 j*l,nncr2 

wyc (i) *wyc (i)+twy (i , j) ftgsmmc (j+1) 
wzc (i) *wzc (i)+twz (i ,j) *gammc (j+1) 

8 continue 

7 continue 

call loopO 
do 11 i*l,nncr2 
wyc(i)*wyc(i)+wyctr (i) 
wzc(i)-wzc(i)+wzctr (i) 

11 continue 
c ft 

25 continue 

c * 

C * X COMPONENTS OF THE INOUCED VELOCITIES 
c * 

do IU i"l,nncr2 
wxc (i) -0. 
wxctr (i) -0. 
do 10 j*l,nncr2 

wxc (1) “wxc (i)+twx (i ,j)*gammcG+l) 

10 continue 

U* continue 

call IcopO 
do 12 i*l,nncr2 
wxc(i)*wxc(i)+wxctr (i) 

12 • continue 

I f (i trcl 1 .eq.O) goto 9 
wr i te liwr , 1 1 1) iter 

111 format {/,"*** L0QP1 *** LAST ITERATI ON:", i3) 

wr i te (iwr ,101) lx , (wxc (kk) ,kk*»l,nncr2) 
wr i te (iwr , 101) ly , (wyc (kk) ,kk*l,nncr2) 
write(iwr, 101) lz , (wzc (kk). ,kk*l ,nncr2) 
write (iwr, 120) 

write (iwr, 101) lyt, (wyctt(kk) ,kk-l ,nncr2) 
wr i te (iwr , 101) lzt, (wzctt(kk) ,kk-l ,nncr2) 
wr i te (iwr , 120) 

wr i te (iwr , 10 1) 1 g , (gammc (kk) , kk«2 ,nncr 2) 
wr i te (iwr, 121) epp 


121 

9 

c ft 


format (/," E?P-",f 10.6,/} 
cont i nue 


119 


i f (i trace.cq. 1 .or . i trcl 1 .eo. 1 ) wr I te (iwr , 103 ) i test 
103 format(" ftftft L00P1 *** : RETURN CODE : ",i3) 

i f (i test.eq. 2 ) call out 2 {l) 
if (itest.aq. 2 ) wr i te (iwr , 101 ,) liml 

101, ’ • format {" *** L00P1 *** ; NO CONVERGENCE WITHIN ",i3," IT r RATI ONS") 

120 format (lx) 

write ( 6 . 7000 ) (gamme (i) , i» 2 ,nncr 1 ) 

7000 format (5x."GAMW>",20f 10 . 5 ) 

write ( 6 , 5000 ) 2 nv{ 13 , 13 ) ,xnv(l3, 1 3 ) ,ynv( 13 , 13 ) 

5000 format (/ ,5x, 'ZNV- ',fl0.5, 'XNV- \f 10. 5, 'YNV- • ,f 10. 5) ORIGINAL PASS IS ' 

re * urn OF POOR QUALITY 


120 


<* 

c* 

e* 


e* 


». 

subroutine loopO 


cownon/roJ 15/wxctr (25) .wyctr (25) .wsct r (25) 
conwon/rot I !C/wxctr2 (25) ,wyct.-2 (25) ,* 2 ctr 2 (2J) 
ccwmon/scfsi/delr 

ccnmon/parm/iwr , i rd, 1 traca, Ipaen.nnaa. i trcw, I trell, 5 trc!2, 1 trcg, 
61 trcnt. i trcf , i plot. icgen, iplotv, i tr tg, I vcrgr , i teai. i a awsft, ! aanet , 
Slims, I ini, I im2.ni tar , 1 trans, ipn.j punch, iplotw, law, i cont, iplotg, 

SIpr, iplott 

ecnmon/gecffl/nb!dsl,nbld8,8lgma,fmu,Qtan(25) .knnvr.etal (6) .knivr, 
61 twist, the t»d (25) . theta (25) , thetac (20 , thctaO. the tOd, alpha*. 
6cdO,cdk.dpsind.dptin,dpsi id.dpsi I ,cocf f , coof f l.c.s.blades.nnvr, 
Snnva.nnvr l,nnval,r.ncr,nnca.etanv(25) .etanc (26) .nivr.niva.nivr I, nival, 
Snicr.nica.etaiv (6) .etaic (7) .ntva.ntvel.ntca, f pa i , fps2, nr 1 1 , nr 1 2, 1 r 1 1, 
Sir 12, imodel 

ccrsnon/reaul/gairenc (26) .twx (21*. 20 , twy (2L.2L) , twt (26, 26) , 

Swxc (20 ,wyc (21*) ,w 2 c (21*) .wxnc (26, 19) .wync (26. 12) .wane (26,19) , 

Swxlc (9.51) .wyic (9.51) .W 2 le (9.50 .wrnv (23. 18) ,wtnv(28, 18) . 

Swznv (28. 18) .wr Iv (8.50 .wtiv(8.50) ,wslv(8.50) . 

Sxnc(26.19) .ync(26.19) ,2ne(26.l9) .xic(9.50 .ylc(9.50 .tic (9.50 , 
Sxnv(28.18).ynv(28.18) .tnv(28.l8) ,xiv(8.50) .yiv(S.SO) , 2 iv (S.50) . ia. lb 
cosaon/vel /x.y , 2 ,ux,uy .ur.xl ,y l ,z l .x2.y2.22,xJ.y3,l3.x 1 »,yU,tl», 
Sgn.qml .gmD.dgn, i t, rho.spf . 1 1 ,epal , epa2 , atra, jgti 
data twopi/6. 283185306/ 




c* THIS SUBROUTINE COMPUTES THE INDUCED VELOCITIES CM THE BLADE 

C* DUE TO INTERMEDIATE AMO FAR WAXE ROLLED UP VORTICES 

c* 

nncr2"nncr-2 
nncr5*nncr2+3 
do 10CC j»3,nncr2 
k«nncr5“j 

I f (aos (gaawc (O ) , I t.abt (ga.nmc (k-l) )) go to 1000 
go to 2CC0 
1000 continue 

2000 qk^-ganrac (k) 

q3«garmc (k) -ganmc (k-3) 
ql-gannc (k-3) -gamme (2) 
lt-2 

Irl2l-lrl2-l 
c ir!22-irl2-2 

nrl2t«nrl2-l 
e nr!22"nrl2-2 

do 3000 i«irl 1, ir 12 

1 f (i .eq . i r 1 1) gn»ql 

C l f (1 .cq. 1 r I 22) gn»q2 

I f (I .eq. i r 1 21) gn*q3 

l f (i .eq. ir I 2) gm»qk 
do 3500 j-l .nival 
dpa7»eps2 
xl-xiv (i ,j) 
yl-yiv(i ,j) 

2 1*2 i V ( i . j) 
x3<*xiv(i ,j + l) 
y3-yiv(; .j-n) 

23 ": iv (i , j + l) 

pnib»twooi/bl adei 
kaoa* : f i x (phi b/aps i i ) 
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3100 

3000 

c* 

c* 

c* 


c 


3200 


nvsp*niva-knca 

phi !*d?si n* (r.nva-U+dpsi i* (j-1) 
if (phi l.gt.pnib) ep*2“.l 
call coord 

If (I trans.eq.O.ar.j .gt.2) call vindb2(L,i) 

eps2°dps2 

continue 

continue 

FOR 2 SPIRALS : 720 degrees 


on«SH» Pf02O 

op peon QO^ Lir 


do 3200 i ■ i r 1 1 , i r 1 2 
if (i ,eq. i r I J) gm-qlM. 
if (i .eq. ir 122) gra»q2*i. 
tf (i .eq. i r 1 2 1) gra«q3*l. 
i f {i .eq. i r 12) gn»ql|ftl. 

2 pf “ 2 i v ( 1 ,nivo) * 2 - 2 iv (I ,nvsp) 
del 2 * 2 iv (! .nivj) “siw (i ,nvsp) 

rho*.5 rt sqrt ( (xi v (i ,niva) +xiv (I .nival) } **2+(yiv (i .niva) + 
6y!v(i .nival)) **2) 

it-3 

call vindb2 (U, i) 

continue 

return 

end 



subroutine out2(‘cut) 
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e ft 

5 * AAAAArtAAArtftAAAAAAAAft 

e *» ft*. ft 

C * * PROGRAM FWC * 

C A ft * 

c * » SUBROUTINE OUT ft 

e ft ft ft 

C ft ftAAAAftAAAAAAAAAAAAAA 

C ft 

c ft 02 / 14/78 

c * 

jsV * OUTPUT OF THE RESULTS 
,e . * 

c <» ENTRY OUTINT: OUTPUT FOR INTERMEDIATE RESULTS 

c • 

real tip.lt. tp, Ipi 

conairon/parn/iwr. I rd, i traca, Ipscm.nme*. I trew, l trel 1 , I trcl 2 . i treg, 
ti trent. i tref , ip lot, iegan, i plotv, i trtg, i vergr , i tost. leemcD, i semat, 
41 ins, 1 1 ml, 1 ioil.ni ter. i trans, ipn, jpuneh, iplotw, i tw, Icont, iplotg, 
tlpr , iplott 
cemnon/hparn/i fmu 

ccnjBon/geoni/nbldst.nbtd-. sigma, fmu.etan ( 25 ) .knnvr.etai ( 6 ) .knivr, 

41 twi at, thatad ( 25 ) .theta ( 25 ) ,thctac{ 24 ) , thataO.thatOd. alphas, 

CcdO, edk, dssind.eps in, dps i id.dpsl i .ccaff .coaf ft, c.s. blades, 
4nnvr,nnva,nnvrt,nnvat,nncr,nnca,at8nv{25) ,etanc(2S) , 
4nivr,niva.nivrl,nival.nicr.nica.etsiv(6) .ataic (7) , 
intva.ntvat.ntca.fpsl.fps2.nr It, nr 12, ir 1 1, i r 12. i medal 
ccrnian /resul/gsirsnc (26) , tv<x (2L, 2L) , twy (24.20 . tvx (24. 2L) , 

Awac ( 24 ) .wyc ( 24 ) ,wzc ( 24 ) ,wxnc ( 26 . 15 ) >wvnc ( 26 . 59 ) .wane ( 26 . 19 ) . 

Swxlc ( 3 . 51 ) .wyic ( 9 , 51 ) .wsic ( 9 . 51 ) .wrnv ( 23 , 13 ) ,wtnv( 23 . 18 ) . 

Swznv ( 23 . IS) .wr i v ( 8 . 50 ) .wt I v ( 8 . 50 ) ,wz I v ( 8 . ; 0 ) . 

Cane ( 26 , 19 ) .ync( 26 .l 9 ) .rnc ( 26 . 15 ) .xic( 9 , 5 U . y ic ( 5 , 5 D .sic ( 9 . 51 ). 
Cxnv ( 28 . 18 / ,ynv( 28 . 1 S) ,xnv< 28 . 13 ) .xiv( 8 . 50 ) .yiv( 8 . 50 ) . 2 iv( 8 . 50 ) , 
ila, ib 

common /otdata/ps i ( 28 ) .r ( 28 ) ,gmc( 24 ) , 

Iwxt ( 25 ) .wyt ( 25 ) . f lamda ( 24 ) ,u( 24 ) ,al pha ( 24 ) , lip ( 24 ) ,flip( 24 ) , 

4 dp ( 24 ) ,fdp ( 24 ) . tdp ( 24 ) ,tp( 24 ) ,fp( 24 ) ,tlip( 2 k) ,xt( 28 ) ,yt( 28 ) . 

4 wt ( 28 ) ,wr ( 28 ) 

CCWHion/terap/l tl 

equivalence (nnvr.nncrl) , (nivr.nicrl) , (nnva.nncal) 
equivalence (niva.nical; , (nnvr l ,nncr 2 ) , (nivr 1 ,nicr 2 ) 
data pi 2 / 9 , 8656044 / 
data pl/ 3 * 14159 * 653 / 

data leta, Iwxc, iwyc, lwic/"ETA ","WXC ","WYC "»"WZC "/ 

data I u , Igmc, I I cm, lnlp/"U "."CMC ","LAM ", ,f ALP "/ 

data Ithe, II ip, Itlp, If !p/"THE ","LIP ","TLIP*'."FUP"/ 

data Idp , I tdp, 1 fdp, I tp /"OP "."TOP ","F 0 P "."TP "/ 

data lotn.lxn ,lyn ,lzn /"ETAN","XNV ","YNV "."ZNV "/ 

data leti.lxi ,lyi ,lzi /"ETAI'VXI V M ."Y1V "."ZIV ’■/ 

data Iwrn, Iwtn, lw:n/"V;RNV", ,, VfTNV","W 2 NV , V 

data Iwri , twti , :w 2 i/"WRIV", ,, WTIV ,, ,"WZIV"/ 

data Ir ,lpsi/"R "."PSt "/ 

data I f p/"FP "/ 

data 1 twy, I twz/"TWY ","TVZ "/ 

c a 

1 10 format (lx.ai, lx,9f 10,6./. 6x.9f 10.6) 

111 format (lx. 34 . Ix.of JO. 3 ,/, 6 x. 9 f 10 . 3 ) 

120 format(U) 

140 format (i 4) 


O* 



c * 
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wr I te ({wr , 100) 

100 format (10*. //."fRfE WAKE ANALYSIS") 
wr i to (iwr , 196) cdO.cdL 

l$6 format!//, lC»,"(C00-",f5. 3.". CDK-", f 5. 3»") ") 

goto uo 

c ^ 

C ftftrtftftftftftftftftrtnnrtrtftrtrtrtrtftftftftrtftftftftrtftftrtftftrtftftftftftftrtftrtftrtftftftft 
C * 

entry outint2 

c ft 

lout-1 

wr I te (Iwr , 101) ni tar 

101 format (lhl,15x,"l-REE WAKE ANALYSIS ", 

t" (INTERMEDIATE RESULTS, ITERATION NUMBER:", 1 3, lh) ) 

c * 

UO continue 

c ft 
C * 

1 1-0. 

Ip-O. 
lpl -0. 

C ft 

C TIP LOSS FACTOR 

c 

gomrac (nncr 1) -0.0 

c 

do 1 l«l,nncr2 

tlarr- (wjc ( i ) mY.ru) / (etanc (i + l) *wyc (I)) 
f lamda (I) -stan u 'am) 
u (1) - (a tone C I -*• 1) -w-yc < « ) ) /co* (f 1 amda (I) ) 
a I pha (i) «f 1 amda (i ) -thatac (i) 
alpha i-al pna (i i 

if (aes (alpha ( : ) ) .gt. alphas .and . alpha (I) .gt.O.) alohal-alpha* 

I f Lbs (a I pna (i) ) . gt .a I phas .and . a I pra ( i ) . 1 t.O.) a I pha 1 --alphas 
I f (i .oq.nncrl) go to 1000 
gmc (I) -pi rt p i fts I gmaft l . ftu (I ) "a I pha 1 /blades 
1000 cd-coO+cdWfta'pha (I) ftft2 

I f (abs (alpha (!) ) .gt. alphas) cc-2.ftcd0+cdkft3lpha (I) ft*2 

1 1 p (i ) -u (i ) ! <gmc t i ) 

tl Ip (I) -I *p (i ) *co» (f lamda (i ) ) 

f 1 1 p (i) — I ip (i) nsin (f lamda (i) ) 

dp (i) -u (I) ftu (i) ft (pi As! gma/b l ados) fted* . 5 

f dp (!) -dp (i ) ftcos (f I amda (i) ) 

tdp(i) -dp (i) <>s i n (f lamda (i) ) 

tp (i) - tl i p (I) +tdp (i) 

fp(i)-f I ip(i) ♦fdp(i) 

1 t-l t+tp (i) ft {etanv (i + j) -otanv (I) ) 
lp-lp+f p (i 5 ft (etanv ( i + 1) fti>2 -o t anv (i ) ft ft 2) 
lpi-lp;+f 1 ip (i) ft (etanv (i + 1) ftftj-etanv (i) ft* 2) 

1 continue 

I t-l t ft b I ades/p i 
lp-lpi' blaaos/ (o •' "2.) 

Ipi - i pi ftbl aces/ (pi ;*2.) 
tap- .0 1 ft 1 1 
tapp-1 1 1- 1 1 

if (abs (taco) . 1 1 .abs .tab) ) wr i te (iwr , 105 ) 
format (//," CONVERGENCE f OR IT REACMCD", /) 

1 t 1-1 t 

if (i imj.cq.O) go to 5 


OHlGlCiAL IS 

OF POOR QUALitV 


105 



124 


*• 


5 


83 


89 

109 
c * 

c * 

115 


12 

II 

10 

c * 


' • *- • • *. 
cp-lp/fmu*P3 
Ct-I t/fnu rtft 2 
ctcp«ct/cp 
conti nuc 
wr I to (iwr , 
wri to (iwr , 
wr!to(iwr, 
wri to (iwr, 
wri te (iwr, 
wri to (iwr, 
wr I to (iwr , 
wri to (iwr, 
wri to (iwr , 
wri to (iwr , 
wri to (iwr, 
writo(iwr, 
wri ta(!wr, 
wri te (iwr, 
wr t to (iwr , 
wri to (iwr, 
wr I to (iwr , 
wr ito(iwr, 
wri to (iwr, 
wr i to (iwr , 
wri to (Iwr, 
wr i te (iwr , 


• « 


10) Iota, (atanc(i) , i-2,nncr l) 

:o) 

10) iwxc, (wxc (i) , l«l,nner2) 

10) iwyc, (wyc(i) , i-l,nncr2) 

10) I w 2 c , (wic(i) , l-l,nncr2) 

20) 

10 ) lu, (u (I) ,I-I,nner2) 

10) lgmc, (gmc (i) , l-l,nncr2) 

20 ) 

10) 1 lam, (f lcmda(l) , I-l,nncr2) 
10) I a I p , (alpha (H , i-l,nncr2) 
10) l tho, (thotcc (l)t I-l,nncr2) 
20 ) 

10) II ip, (llp(i) , i-l.nncri) 

10) l tip, (tl ip (i) , i-l ,nncr2) 
10) 1 f Ip. (fllp(i) ,1-I,nncr2) 
10) 1 tap , (dp ( l ) , i«» 1 . nncr 2) 
10)ltdp , (tdp(lj , i-l,nncr2) 
10) I f dp , (fdp(lj ,I-l,nncr2) 

, (tp (lj ,i»i,nncr2i 
t (fp (i) , l-l,nncr2) 


1C) I tp 
10) lfp 
20 ) 






if (lfmu.no. 0) wrlto(lwr,88) 1 1, Ip.ct.cp.ctcp 

If (Ifmu.cq.C) wr i te (iwr ,83) It, Ip 

format (/// , 10x."LT -",el2.5/ , .Ox."tP « ,, .el2.5/10x."CT 

6/ lOx . "CPS-" . c 1 2 .5//9*"CT/CP-" , o 12 .5//) 

wri to (iwr ,89) ) p i 

format (///, l0x,"LPI -",012.5) 

wri to (iwr . 109) 

format (lh 1) 


•". 012.5 


if (lout. eq. 2) return 


continue 

lw»l 

nvr-nr 1 2 

write (Iwr, 110) letn, (otanv(l) , l-l,nnvr) 
do 10 J-l,nnva 
wr 1 te (iwr , 11,0) j 

wri to (iwr , lie) I xn , (xnv(l , j) , 1-nr 1 1 , nr 1 2) 
wr I to (iwr , 110) lyn , (ynv(l , J ) , I -nr 1 1 , nr 1 2) 
wri to (iwr, 110) Ian , (anw(i.j) ,1-nrl l . nr 1 2) 
wr l tc (Iwr , 120) 

wri to (Iwr , 1 10) lwrn, (wrnv(i , j) , I -nr ) I, nr 1 2) 

wri to (iwr , 1 10) Iwtn, (wtnv (i , j) , i-nr I !,nrl2) 

wr i to (iwr . 1 10) lwan, (wanv (i , j) , 1 -nr 1 1 , nr 1 2) 

do 12 i-nrll,nrl2 

wr (i ) — wrnv (» . j) 

wt (i) -wtnv (1 , j) 

xt (I) -xnv (i , j) 

yt (i) -ynv (i . j) 

conti nue 

goto 30 

cont i nue 

conti nue 
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lw»2 

nvr-ir!2 

wr ! te (Iwr , 110) let i , (etniv (I) , !»1 ,ni vr) 

* - do 20 j«l,niva 
wr i tc (iwr , 11,0) j 

wr i te H wr , 1 10) ix i , (xiv(i ,j) , i-lr 1 1, ir 12) 
wr i te (iwr , 1 10) 1 y i , (yi v (i , j) , i-i r 1 1 , i r 12) 
wr I te (iwr , 1 10; lzi , (ziv(i ,j) , i-irl 1, Ir 12) 
wr i te (iwr, ? 20) 

wr 1 tc (iwr , 1 10) Iwr i , (wr i v (i , j) , i«*i r 1 1 , Ir 12) 
wr i te (iwr , ! 10) Iwt I , (wt iv (i , j) , i-lr 1 1 , Ir ! 2) 
write (iwr; 1 10) Iwzi . (wzl v (1 , j) , i-lr 1 1, ir 12) 
do 22 1— i r 1 1 . ir 12 
wr (!) -wr iv(i , j) 
wt (i) -wt iv (i , j) 

X t ( i ) “X I V ( i , j) 

jrt(i)-yiv(i,j) 

22 contir.ue 

goto 30 

21 continue 

20 continue 

1.1*3 format (8f 10. 5) 

UUU format (2f8*l*) 

1,2,5 format(6M0.5) 

C * 

wr i te (iwr , 203 ) 
do 93 I-l ,nncr2 

wr I te (iwr .201) I twy, (twy (i ,k) ,k-l ,nncr2) 
wr i te 'iwr ,201) 1 twz. (tws (i ,k) ,k»l,nncr2) 
write (iwr, 202) 

93 continue 

201 format (lx. au, lx,9f 1 3 - 9/ox . Sf 1 3 • 9> 

202 format (//) 

203 format (" INFLUENCE COEFF. TW ( I , K) A5AM/1A (K) -Wl ( I ) “/) 
c * 

i f (i model .ne.2) go to 201. 
do 55 j-l.nnva 
write (iwr, 1 1,0) j 

wr i te (iwr , 110) ixn , (xnv(i , j) , 1-1 ,nnvr) 
wr i te (iwr , 1 10) lyn . (ynv(i ,j) , I-l, nnvr) 
wr i to (iwr, 1 10) l 2 n , (znv (i , j) , i-l ,nnvr) 
wr i te (iwr, 120) 

wr i te (i wr , 1 10) Iwrn, (wrnv (i , j) , i -l , nnvr) 
wr i te (i wr , 1 10) iwtn, (wtnv (i , j) , i-l , nnvr) 
wr i te (iwr , 1 10) Iwzn, (w 2 nv (i , j) , i -I ,nnvr) 
do 56 i-l, nnvr 
wr (i) -wrnv (3 , j) 
wt (i)-wtnv(i , j) 
xt (i) -xnv (i , j) 
yt (i) -ynv (i , j) 

r (i) -sqrt (xt (i) **2+yt (i) **2) 
cc-xt (i) /r (i) 
ss«yt (i)/r (i) 
phi-atan (yi (i) /xt (i) ) 
i f (xt (i) . I t.0.) phi-phi+pi 

if ((xt (i) ,gc.0.) .and. (yt (i) . 1 1 . 0 .)) phi-phi+pi+pi 

psi (i) "57.2?6"'phi 
conti nuc 
wr i to (iwr , 120) 

wr i te ( i wr , 2 1 1) I psi , (ps i (i ) , i - 1 , ,mvr) 


OKiGSKAL PAGS U* 

o? poos QUAJ.rr/ 


56 
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wrlte(!wr,UO) lr . (r (!) , I-l.nnvr) 
wrt to (iwr , 120) 
continue 
return 

continue 

.If(iw.eq.l) Ij-nrll 
. i f (iw.eq . 2) lj-Irll 
do 31 I-Ij.nvr 
r (!) -sqrt (xt (1) "*2+yt (!) **2) 

CC«xt(D/r (!) 

•«»yt (!) /r (!) 
phi"atsn(yt(I)/xt{!)) 

If (xr (!) . 1 t.O.) phi-phi+p! 

If ((xt(!) .ge.O.) .and. (yt(!) .It.O.)) phl-ph!+p.+pi 

ps! (i)-57.296*phi 

continue 

wr ! to (iwr , 120) 

wr l to (iwr , 1 1 1) Ipsi • (pal (!) • 1*1 J «nvr) 

write (iwr, 110) lr , (r (!) , l-IJ .nvr) 

write(iwr,l20) 

if(iw.eq.l) goto 11 

goto 21 



end 
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#• 


Of 


!>ubscrinced Variables 


TABLE 2 


of peon 


etanv 


Radial positions of Che blade stations and Che near wake 
vortices before they deform 


etanc 


Radiol positions of the blade centers 


gammc/Gnu- 


The bound circulation distribution 


gummtl 

garont2 


grain 

gra2n 


hammtl 

hammtl 


hmln 

hra2n 


Strength of the root and tip vortex in the near wake, due to 
a unit circulation distribution 

First and second lumped strengths of the vortex sheet elements 
in Che near wake 

Coefficients for ganatl, gamnc2 respectively 

Note: by definition: gammtl(i) ■ 2 hamcitl (i,j)* 

gamme (j+1) 

Coefficients for grain, grain respectively. 

Note: by definition: graln(i) <* 2 hmln (i,j) * gararac(J+l) 


twx 

twy 

twz 


X,v,z components of the influent coefficient. They arc the 
induced velocities at the blade centers due to a unit bound 
circulation 
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vrriv 
. >' . 

vtlv ■ 
vziv 


Radial, tangential and axial components of the induced velocity 


1 1 ... V 


:it. at node (i.J) of the intermediate waka vortices 


wmv 

wtnv 

vznv 


Sane as wriv, vtiv, wziv, but for intermediate wake vortices 


vzc 

vyc 

wzc 


X,y,z components of the total Induced velocity at the blade 
centers. For example: wxc(i) ■ wxctr(i) + 2 twx(i,j) * 

gammc(J+l) 


wxetr 

wyctr 

wzetr 


X,y,z components of the induced velocity at the blade cancers 
due to the intermediate and far wake vortices 


wicivr 

vyivr 

wzivr 


X,y,z components of the induced velocity at rolled up vortices 
in Che intermediate wake 


wxnc 

wync 

wane 


X,y,z components of the induced velocity at centers of the 
near wake (at coordinates xnc.ync, znc) 


vxnvr 

wynvr 

wznvr 


X,y,z components of the induced velocity at the anticipated 
rolled up vorcices in the near wake 
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xiv 

yiv 



xnc 

ync 

znc 

xnv 

ynv 

znv 


■* t» 


OF P0Q3 (JUAi-VTi 


X»y,z cartesian coordinates of the node of intermediate wake 
vortex 


X,y,z cartesian coordinates of the center of near wake vortex 


Same as xiv, yiv, ziv but now these are for the near wake. 

Xnv is the x coordinate of the node, located on the stream lines 
originating from the blade stations (etanv) 



r 


mr < 
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variables 

alphas 

capk 


'* p . G c \S 

stall angle /In radians) OR\G^^ r\\\P&X^ 

• ■ • ' of p° oR u 


elliptical integral of the first kind 


fj 

cdo, cdk drag coefficients (cd » cdo + cT cdk) 


cgl 

eg 2 

cg3 

cg4 


spanswise location of the rolled up root vortex 

spanovise location of Che rolled up 2nd middle vortex 

spanswise location of the rolled up 1st middle vortex 

spanswise location of the rolled up tip vortex 


at t> » 0 


coeffl/coeff parameter of root and tip vortices 

(parameter k) * 


dels 


Vertical spacing between the la3t node of the intermediate 
wake and the starting location of the far wake. 


dpsiid 


azimuthal size of each vortex element in the intermediate 
wake 


dpsind 


sane as dpsiid but for the near wake vortex element 


e 


elliptical integral of the second kind 
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epsl 

eps2 

etan 

frau 

fpsl 

fps2 

imodel 


irll 

irl2 

irl21 

ir!22 

it 


itest 




• 2 second vortex sheet model 
■ 3 ; third vortex sheet model 
** A ; fourth vortex sheet model 


rolled up root vortex 
rolled up tip vortex 
rolled up first middle vortex 
rolled up second middle vortex 

m 1 . the vorcex is a sheet/rectangular element 

“ 2 ; the vorcex is a filament/segment element 

- 3 : the vorcex is a slmi infinite cylinder 

■O' convergence is reached 

** 1 . convergence i3 not reached, but the iteration 
on the bound circulation has converged 


all in. the 

intermediate 

wake 



2 , the iteration on the bound circulation does not 
converge within the liml limit. 


kunvr 


liin 1 


lira 2 


Ltwist 


nblds 


niva 


nicer 


nnca 


nncal 


nncr 


nncrl 


nr.cr2 


■I l ' -t ' ■ 1 wyyiT 1 i wr? i n L ‘ » m.i • v.’H ’ .u ■y n w-’i u ' . i n r i ■ " » ' ff 

• « tl f» ' 

«»»*«£ •• *.»• 

of f 00 

sane as nnvr. 

limit of the number of iterations on the bound 
circulation 

limit of the number of iterations on the wake 
geometry 

■1 I twisted blade 
■0 * untwisted blade 

■ number of rotor blades 

m number of axinuthal nodes in the intermediate wake 

■ current loop index for the iteration on the wake 
geometry 

■ it is the number of centers in the azimuthal direction 
of the near wake (equal to nnva +1) 

equals to nnca -1, equivalent to nnv 2 

it is the nuober of centers in the radial direction 
of the blade and near wake. 

equals to nncr -1, equivalent to nnvr 

equals to nncr -2, equivalent to nnvr 1 
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nnva 

nnval 

nnvr 

nnvr 1 

nrll 
nr 12 
nrL’.i 
nr 122 
Note: 

ntva 

qi 

q2 

q3 

rho 
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number of nodes in the azimuthal direction of the 
near wake. 


equals to nnva -1 


number of stations/nodes in the radial direction 
of the blade and the near wake. 


equals to nnvr -1. It is the number of vortex elements 
of the near wake (due to tip loss) 

anticipated rolled up root vortex \ 

anticipated rolled up tip vortex I in the 

anticipated rolled up first middle vortex l near wake 

anticipated rolled up second middle vortex J 

nrll must be higher than nnvr. If there are 

only 3 rolled up vortices, then the nrl22 is not used. 


same as nnva 

strength of the rolled up root vortex 
strength of the rolled up second middle vortex 
strength of the rolled up first middle vortex 
strength of the rolled up tip vortex 


radius of a vortex cylinder in the far wake 
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a 1 Rail 


thecad 


ux 

uy 

us ] 
wt 


wll 


wls 


x.y.s 


* 2 '^ 2 |S 2 

x 3* y 3 ,s 3 


x, ,y. z. 


:p: 


solidity of rotor 

tha pitch angle at each spansvlsa station (In degrees) 
For untwisted blade, thotdd equals to the pitch 
angle (In degrees) suet be specified. 

x.y.s components of the Induced velocity 
due to one vortex eleaont. 

vortical velocity Induced at a segment elesent due to 
a vortex ring (excluding the self-contribution). 

lamb vortex ring self-induced velocity. 

equals to vtl-vl. Is the self-eleaent vertical induced 
velocity. 

point at which the Induced velocities are computed 
coordinates of the first point on thu rectangular 
on the segment element 

coordinates of the second point on the rectangular 
element 

coordinates of the third point on the rectangular 
element or the second point on the segment element 
coordinates of the fourth point on the rectangular 
element 

starting }\iai location o: the vertex tv Under 


of vo° R w 
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Variables 


ORIGINAL PAGE IS 
OF POOR QUALfTV 


Table 3 


ALP 

DP 

ETA 

FdP 

FP 

CAMM 

LAM 

LIP 

LP 

LT 


angle of attack 

drag distribution 

location of blado centers 

In plane component of the drag force 

thrust distribution 

bound circulation distribution at the blade centers 
In flow awlc 

in plane cotaponent of the lift 

p 

power coefficient ■ v .t — 

thrust coefficient • — A ■■ ■ ■ - 

PTli- R 


?SI 


R 

THET 

TUP 

T? 


axinuthai position (in degrees) 
note: the range is frota 0 to 360 degrees 

radial position of the node in the particular PSI 
pitch angle 

thrust component of the lift 

In plane component of the force distribution 


TVX 


TVY 

T.sZ 


Same as in Ttble 2 


U 


resultant of cne induced velocity relative to the 


blado 



WRIV 


v 


VTIV 

vziv 


« 

Sane as In Table 2 



UK,*V 

WTKV 

WZJJV 


» * 

San* as in Tabic 2 



UX 

WY 

VZ 


x »7** component of the Induced velocic 7 at the 
blade centers 


XIV 

TIV 

ZIV 


Sane as in Table 2 


XNV 

YNV 

ZNV 


Sane as in Table 2 



Table 4. TYPICAl. OUTPUT 


\ 
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